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CAPÍTULO 1 


CAUSAS DO ALARGAMENTO 
E DO DESLOCAMENTO DAS 
LINHAS ESPECTRAIS 


1.1 Introdução 


Várias são as causas do alargamento e do deslocamento 
das linhas espectrais. Citaremos algumas delas e apresentaremos 
apenas as expressões finais. As avaliações numéricas poderão 
ser vistas em Cattani (1968). 

Podemos citar, por exemplo, as seguintes causas 
[Townes e Schawlow (1955) e Cowley (1970)]: 

a) Alargamento e deslocamento naturais e induzidos. 

b) Efeito Doppler térmico. 

c) Colisões moleculares. 

d) Colisões com as paredes do recipiente. 

e) Efeito de saturação. 

f) Ffeitos Stark e Zeeman. 


a) Alargamento e deslocamento naturais e induzidos 


O alargamento natural e induzido deve-se à interação do 
sistema emissor com o seu próprio campo e com o campo 
externo, respectivamente. 

No caso de uma transição de freqiiência V, de um estado 
excitado para um estado fundamental, a semilargura da linha 
pode ser avaliada por [Heitler (1954)]: 


py já | [+ N(v)]. (1.1.1) 


em que d é o momento de dipolo do sistema, V a frequência da 
transição e N(v) o número de fótons do campo externo com 
frequência V. No caso de interação do sistema com seu próprio 
campo (alargamento natural), N(v) = O. Para termos apenas uma 
ordem de grandeza do efeito, vamos calcular Av usando a 
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expressão (1.1.1). Assim, para comprimentos de onda 4-1 cm 
(microondas) temos Av - 107 Hz e para transições na região de 
luz visível, ou seja, À - 10 cm, temos Av - 10º Hz. 

O alargamento e o deslocamento naturais são da mesma 
ordem de grandeza. 


b) Efeito Doppler térmico 


Consideremos uma molécula emissora que se move com 
uma velocidade v e que em um dado instante emite um fóton. 
Por efeito Doppler a frequência (V) do fóton emitido por um 
átomo com velocidade v ao longo da linha de observação é dada 
por: 


v=vd+), (1.1.2) 


em que Vo é a frequência ressonante sem o efeito Doppler. 

Se o sistema emissor é um gás ideal a uma temperatura T 
constituído de moléculas que têm massa M, as velocidades v 
dessas moléculas obedecem a uma distribuição de Maxwell P(v) 
dada por: 


1/2 2 
ro -(5] x exp[ | 
27KT PET) q13 


Neste caso, devido ao movimento térmico, o perfil da 
linha T(v) gerado pelo efeito Doppler é dado por: 


2 
Mc? [v—v 
19=com FEL 5 (1.1.4) 


que é uma linha gaussiana simétrica. 
E fácil ver, usando a expressão (1.1.4), que a 
semilargura da linha é dada por: 


1/2 
Apt [27 ema) | (1.15) 
Cc 


Pode-se mostrar também que devido ao efeito Doppler a 
linha emitida centra-se numa fregiência Av, —-R, em que 


R=(hv,) M2Mc”) e c é a velocidade da luz no vácuo. Ou seja, 


ela vai sofrer um deslocamento de frequência dado por R. No 
caso de 4-1 cm, T = 300 K eM - 50 obtemos Av-25 


kHz e R/h-3,6x10” Hz. 


c) Colisões moleculares 


Suponhamos que temos um gás numa dada pressão e 
temperatura. Se de alguma maneira excitarmos os estados 
internos moleculares de um gás, quanto tempo esses estados 
excitados demorarão para retornar ao estado fundamental? Se 
excitarmos níveis ópticos, as colisões moleculares serão pouco 
eficientes na restauração do equilíbrio térmico que será em geral 
atingido após um tempo igual ao da vida média natural. 

Se, no entanto, forem excitados níveis rotacionais (por 
exemplo, níveis na região das microondas), as colisões 
moleculares são capazes de restaurar o equilíbrio térmico. No 
capítulo seguinte estudaremos as condições que devem ocorrer 
para que o equilíbrio térmico seja restaurado. 


d) Colisões com as paredes do recipiente 


As colisões das moléculas com as paredes do recipiente 
onde está contido provocam transições que alteram as vidas 
médias dos estados excitados. Tratando as colisões com as 
paredes como se fossem colisões intermoleculares, pode-se 
mostrar que a semilargura da linha é dada por: 


1/2 
Av = af E (1.1.6) 
v ist M 


em que R é a constante dos gases ideais, A é a área das paredes 
do recipiente e M a massa molecular do gás. Conforme veremos 
a seguir, as colisões intermoleculares geram uma semilargura 
dada por Av=1/(277), em que Té o tempo livre médio entre 
duas colisões. 

Para T = 300 K,M - 50 e estando o gás contido numa 
guia de onda com uma parede de 3 mm x 5 mm, teremos 
Av - 1,5 kHz. 


e) Efeito de Saturação 


Campos eletromagnéticos muito intensos tendem a 
igualar as populações dos níveis entre os quais ocorre transição, 
enquanto as colisões moleculares tendem a restabelecer o 
equilíbrio térmico. Esses dois efeitos são competitivos. 

Sendo I a intensidade do campo incidente, a semilargura 
é dada por: 


Av 


24412 
: [+ 1ómIdP2nti | cita 


“ DRT 3ch? 


/ 


em que 7 é o tempo livre médio de colisão, d é o momento de 
dipolo da transição. 

Este efeito foi estudado por Karplus e Schwinger (1948) 
em que se vê que é a própria radiação incidente que provoca a 
desexcitação das moléculas. 


f) Efeitos Stark e Zeeman 


Quando não existem campos estáticos aplicados (elétrico 
ou magnético), os níveis de energia de um sistema emissor são 
degenerados. A remoção dessa degenerescência ocorre quando 
há um campo externo estático — efeito Stark — se o campo for 
elétrico, — efeito Zeeman — se o campo for magnético. Os níveis 
de energia são separados e, se esta separação for da ordem da 
largura original da linha degenerada, observar-se-á um 
alargamento da largura original da linha. 


Referências — C1 


Cattani, M. 1968. Tese de Doutoramento, Universidade de 
São Paulo (não publicada). 


Cowley, C. H. 1970. The Theory of Stellar Spectra, 
Gordon and Breach Science Publishers. 


Heitler, W. 1954. Quantum Theory of Radiation, Oxford 
University Press. 


Karplus, R. and Schwinger, J. 1948. Physical Review 73, p. 
1020. 


Townes, C. H. and Schawlow, A. L. 1955. Microwave 
Spectroscopy, McGraw-Hill Book Company, Inc. 


CAPÍTULO 2 


TEORIAS CLÁSSICA E SEMICLÁSSICA 
DO ALARGAMENTO E DO 
DESLOCAMENTO PROVOCADOS 
POR COLISÕES MOLECULARES 
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2.1 Introdução 


As teorias desenvolvidas para explicar o alargamento e 
deslocamento de linhas espectrais provocados por colisão são 
divididas em dois grupos: teoria do impacto e teoria estática ou 
quase estática (impropriamente denominada de estatística). Na 
primeira, considera-se que as moléculas estejam a maior parte 
livre e que as interações entre elas ocorram em um intervalo de 
tempo muito menor do que o tempo livre médio de colisões (7); 
na segunda, as moléculas estão paradas e interagindo entre si 
continuamente. 

Para um melhor entendimento da teoria do impacto, 
vamos dividí-la em duas partes: Teoria das Colisões Fortes e 
Teoria da Integral de Fourier. 


2.2 Teoria das Colisões Fortes 


Os trabalhos desenvolvidos para explicar o alargamento 
de linhas espectrais tiveram início com Lorentz (1906), 
seguindo-se os de Debye (1945), Van Vleck e Weisskopf (1945) 
e Karplus e Schwinger (1948). Esses trabalhos consideravam 
fundamentalmente a idéia de que as colisões moleculares eram 
capazes de restaurar o equilíbrio térmico perturbado pela 
radiação incidente. Essa hipótese é satisfatória apenas na região 
de microondas [Cattani (1968)]. 


No modelo desenvolvido por Debye (1945), o 
Hamiltoniano de interação entre a molécula emissora e o campo 
extemo É é dado por H = d.É,em que déo dipolo elétrico 
da molécula. Em seu trabalho, Van Vleck e Weisskopf (1945) 


12 


obtêm o resultado de Debye, mostrando que o coeficiente de 
absorção (ot) da radiação incidente por unidade de comprimento 
é dado por: 


O 4mN gp q (2.2.1) 


em que w é a fregiiência (angular) da radiação, 7 é o tempo 
livre médio entre colisões e N é o número de moléculas por 
unidade de volume. 

Pela expressão (2.2.1) vê-se que a teoria de Debye é não- 
ressonante pois ela não depende das frequências próprias do 
sistema. Fisicamente significa que não há transições internas, 
uma vez que o campo externo faz apenas girar a molécula 
mantendo |J | (módulo do momento angular) sempre constante. 
É fácil ver que a linha de absorção de Debye é centrada na 
origem (mw = 0). 


Em seu trabalho, Lorentz (1906) substituiu a molécula 
emissora por um oscilador harmônico de frequência própria 0 e 
assumiu que o Hamiltoniano de interação entre a molécula 
emissora e o campo externo (radiação incidente) é dado por: 


2 


2 
H = > 5 m (mx) — e x E cos(0t) . (220) 


Para resolver a equação do oscilador harmônico, Lorentz 
assumiu a hipótese de que, logo após a colisão, valores positivos 
e negativos do valor médio da posição e da velocidade são 
igualmente prováveis, isto é <x>=< x >=0, o que equivale 
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a admitir que a molécula após cada colisão é não-polarizada. 
Como as colisões não ocorrem ao mesmo tempo, temos: 


«>= [E expf- 8) do. (2.23) 


Por sua vez, Van Vleck e Weisskopf (1945) assumiram 
que após o choque o oscilador entra em equilíbrio térmico e 
calcularam os valores médios < x > e < x > usando a 
distribuição de Maxwell-Boltzmann. O coeficiente de absorção 
obtido por eles é dado por: 


a = 28N 62 (2) mo lx 


meo) (oo) +17) (oro) +17) | 
(2.2.4) 


[Observe-se que Van Vleck e Weisskopf obtêm o mesmo 
resultado de Lorentz. Porém, quando wo = 0, não conseguiram 
obter o resultado de Debye]. 


A expressão (2.2.4) é a forma clássica da linha de 
absorção. Para estender essas expressões ao caso semiclássico, 
bastará fazer a seguinte correspondência [Sakurai (1967)]: 


e? 87N -2 
5 >v, 
m 3h * 


ES MR? 


if = 2 


E 


O, D27V, =0, = (Wj-W)/h. 
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Essas correspondências equivalem a identificar cada 
transição entre estados estacionários (i — j) como um oscilador 
“virtual”. Nessas expressões Idil é o elemento de matriz do 
momento de dipolo (elétrico ou magnético) relacionado com o 
processo da transição (1 > j); w; é a fregiiência de Bohr 
entre os estados (1) e ()) e Wi e Wj; suas respectivas 
energias. 

Se as moléculas entre os vários estados estacionários são 
distribuídas de acordo com a lei de Maxwell-Boltzmann e 
somando sobre as possíveis transições, teremos: 


E Id;É f(00,0) exp [-W; /(KT)] 


3 hc L exp -W, / (KT) | 


= 01;;(09), 
ij 
(2.2.5) 


em que o fator de forma f(w;;.0) é dado por: 


[0) Aq 
0; (co; -o) +(A9) (co; +09) +(A9) 


(2.2.6) 


com Aq = 1/7. 

A expressão (2.2.5) dá o coeficiente de absorção para o 
caso de um estado, ou seja, aquele em que apenas um estado é 
perturbado pelas colisões. No caso de dois estados, sendo a 
probabilidade de absorção a mesma para cada caso, a absorção 
líquida será: 


15 
Oj(60) = o Id, (00,0) (N;-N,) . (2.2.7) 


Se o feixe de radiação é fraco, o equilíbrio térmico é 
restaurado pelas colisões [Ch'en e Takeo (1957)]: 


Tia Po 2.2.8 
N, “PU KT (2.2.8) 


e, portanto: 


“Mov GP ho; 
Cj(0) = SEO N, ldyj f (00,0) ' E es] e) (029) 


A expressão (2.2.9) é a forma quântica ou semi- 
clássica da linha de absorção de Van Vleck-Weisskopf. 


As fórmulas de Debye e de Lorentz são obtidas do perfil 
de Van Vleck e Weisskopf da seguinte maneira: 


a) Perfil de Debye 


Se, ho, <<KT, chamada região de microondas, é fácil 


ver que [Townes e Schawlow (1955)], fazendo di =de 


O = 0 6 Ao >> q; 


;j» Obtemos o perfil de Debye: 


40N po qr 
o = Idl . (2.2.10) 
3KTe 14 (om) 
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b) Perfil de Lorentz 


Se hw,;>>KT, chamada região óptica, é fácil ver que 
sendo 


[Oj— OI<< O, 
então: 


Otorentz = Olvan Vleck e Weisskopf » 


na região de ressonância óptica. 

A explicação física para este fato é a seguinte: como a 
frequência de oscilação é extremamente alta, a fase oscila 
rapidamente durante uma colisão fazendo com que o equilíbrio 
térmico não seja restaurado após a colisão. Neste caso a hipótese 


da não-polarização de Lorentz é válida. 


Em seu trabalho, Karplus e Schwinger (1948) obtêm o 
perfil de Van Vleck-Weisskopf por intermédio da matriz 
densidade, considerando: 


P(D = poll + DCD, (2.2.11) 


Pol = expf- gi / [tr [exp(- E) (2.212) 


ig 


indica o equilíbrio térmico no instante t, H o Hamiltoniano do 
sistema e D(t) representa o desvio da posição instantânea de 
equilíbrio, que contribui para o coeficiente de absorção, e 
resolvendo a equação de evolução temporal da matriz densidade: 


2 p(0 = [Ep] / (ih). (22.13) 


As teorias vistas acima, chamadas Teorias das 
Colisões Fortes apresentavam alguns resultados experimentais 
positivos razoáveis e certas limitações que passaremos a analisar 
rapidamente: 


a) Resultados positivos [Townes e Schawlow (1955)] 


D O perfil teórico de Van Vleck-Weisskopf ajusta muito 
bem a forma das curvas experimentais a pressões 
suficientemente baixas; 

ID) Para um grande “range” de pressões baixas, a semi- 
largura Av depende linearmente da pressão: 


fia E SM NR E 


TT o : (2.2.14) 


em que T é o tempo livre médio de colisão, <V > a velocidade 
relativa média e O é a secção de choque; 


II) A frequência aparente de ressonância é constante 
para um grande “range” de pressões baixas; 


IV) Em altas frequências, a absorção é constante. 
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b) Limitações da teoria 


D Os resultados experimentais para a semilargura AV 
são maiores que 1/(2717), portanto AV é um parâmetro empírico 
que é ajustado para concordar com os dados experimentais; 

I) A hipótese da interação forte não permite introduzir 
os potenciais de interação entre as moléculas colidentes que, 
como veremos, são responsáveis pelo alargamento e 
deslocamento da linha; 

NI) Não é previsto nenhum deslocamento das 
frequências de ressonância que já haviam sido detectadas 
[Bleaney e Loubser (1948)]; 

IV) As hipóteses formuladas só são válidas para a região 
de microondas e infravermelho longínquo, onde as energias 
cinéticas dos átomos perturbadores são muito maiores que as 
diferenças de energia entre os estados que contribuem para as 
transições [Cattani (1968)]. 


Pelas razões expostas acima, teorias mais completas 
deveriam ser desenvolvidas. A primeira que procurou corrigir 


essas limitações foi a chamada Teoria da Integral de Fourier, 
que passaremos a estudar. 


2.3 Teoria da Integral de Fourier 


2.3.1 Teoria clássica 


Para o desenvolvimento desta teoria vários trabalhos 
foram publicados, entre os quais destacaremos os de Weisskopf 
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(1933), Foley (1946), Lindholm (1942), Lenz (1933), 
Mizushima (1946) e Reinsberg (1938). 


Em seu artigo, Weisskopf (1933) formulou um modelo 
para a molécula emissora no qual esta era considerada como um 
oscilador clássico cuja freqiiência variava no tempo de acordo 
com as perturbações provocadas pelas colisões que ele admitia 
serem fracas (adiabáticas). Classicamente essa colisão significa 
uma mudança de fase (“phase shift”) na molécula emissora 
depois da colisão e, quanticamente, uma mudança nos níveis de 
energia da molécula emissora após a colisão. 

A mudança de fase do oscilador é dada por: 


8 = | Ab dt, (2.3.1.1) 

AT 
em que AT é o tempo de duração da colisão e Aot) é a variação 
da fregiiência durante e entre as colisões. Segundo Weisskopf, a 
soma de pequenas fases daria uma mudança de fase — 0 — grande 
que afetaria o alargamento e deslocamento das linhas espectrais. 


Portanto, sendo a lei de movimento do oscilador dada 
por: 


X(D =x es + fot) | (2.3:1.2) 
e sendo [Heitler (1954)]: 


2 
K(0) = e ESCOTA (2.3.1.3) 
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uma transformada de Fourier de x(t), virá [Ch'en e Takeo 
(1957)): 


2 
2.4 

e O 2 

(co) = RCE Xo 
31c 


—oo 


oo t 
fat es + ot) dt' + it 


(2.3.1.4) 
Podemos agora formular as seguintes hipóteses: 


a) Se a oscilação é completamente interrompida a cada 
colisão, ou seja, o trem de onda emitido é finito, o que equivale 
a fazer a integração det = O a t=% e como as colisões não 
ocorrem ao mesmo tempo, uma média de I(w) sobre a 
distribuição de x deverá ser feita, então, a forma do perfil de 
Van Vleck-Weisskopf (1945) é obtida; 

b) A mudança de fase — O — poderá ser considerada como 
um interruptor temporário do processo de radiação, sempre que 
for superior a um valor crítico mínimo arbitrário. Com esse 
valor crítico assumido, o diâmetro de colisão e a semilargura 
podem ser obtidos em função das interações moleculares. 
Weisskopf (1933), Lenz (1933) e outros obtiveram resultados 
semelhantes para O = 1,71/2 e 1.1 radianos, respectivamente. 


[Note-se que Reinsberg (1938), Lindholm (1942) e 
Foley (1946) removeram a arbitrariedade do valor crítico 
mínimo de O incluindo todos os valores de deslocamento de 
fase]. 


Foley (1946), seguindo as mesmas hipóteses de 
Weisskopf (1933), constatam [vide o excelente “review” de 
Ch'en e Takeo (1957)] que a intensidade da radiação I(w) pode 
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ser obtida através da transformada de Fourier da função de 
correlação d(T) dada por: 


AT) = Í dt, esaf f sata (2.3.1.5) 
Portanto: 
I(w) = 2º o* x; Re Í exp (1Awr) 67) dr. (2.3.1.6) 


Para resolver esta integral Foley considera a 
aproximação por diferença de fase na qual as durações das 
colisões são extremamente curtas, binárias e independentes. 
Assim: 


exp | Hat? dt'| = =epli>do (2.3.1.7) 
ih 5 


1=0 


para N colisões. Se p(0)dO/t, é o número de colisões que 
provocam um deslocamento de fase entre O e O + dO, em que To 
é o tempo livre médio entre N colisões no intervalo de tempo 7 
dada pela distribuição de Poisson, então: 


N 
[5] E expf- a), (2.3.1.8) 
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Desse modo, Foley (1946) demonstrou que: 


K(0) = 20: o (2.3.1.9) 
Tic [(1-4)/7.| +[B/t,— Ad) 
onde: 
A = [ p(8) do cos 6, (2.3.1.9) 
B= Í p(9) dO sen 0, (2.3.1.90) 


—oco 


e, portanto, a semilargura e o deslocamento são dados por: 


Aq = 2 [ p(8) (I-cos 6) dg/r 


—oo 


(2.3.1.10) 


o? 


oo 


Am, = | p(9) sen 6 de/7, . (2.3.1.11) 


—oco 


[Note-se que a expressão explícita de p(0)d9/7T, é vista 
no artigo de Ch'en e Takeo (1957)]. 

A mudança de fase 0 pode ser calculada a partir da lei de 
interação molecular. Em muitos casos o potencial de interação 
pode ser tomado como: 
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ER Ea (2.3.1.12) 
Rº 
em que Y é uma função do estado na qual se encontra a molécula 


emissora e R(t) a distância, num instante t, à molécula 
perturbadora, então: 


r Ay 
; 2.3.1.13 
= Rs ( ) 


Na expressão acima, o parâmetro p toma os valores 


p=4,5, ... no caso de a interação entre as moléculas ser dipolo- 
dipolo, dipolo-quadrupolo, ... . [Ch'en e Takeo (1957)]. 


Limitações da teoria clássica da integral de Fourier 


D Não inclui as colisões capazes de provocar transições 
internas, o que classicamente equivale a dizer que a amplitude 
do oscilador varia, colisões essas chamadas de fortes [Anderson 
(1949)]; 


IN) Para colisões fortes (0 > 1) o deslocamento não é 
previsto. Portanto, teremos: (AWm = 0); 


II Ela é válida apenas para baixa densidade e altas 
temperaturas [Foley (1946)]; 


IV) A razão entre a semilargura e o deslocamento 
provocados por moléculas diferentes das emissoras não depende 
da pressão [Ch'en e Takeo (1957)]. 


Desse modo, verificou-se que a Teoria Clássica da 
Integral de Fourier aplicava-se razoavelmente bem na região 
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óptica e na infravermelha, não dando entretanto bons resultados 
na região de microondas [Townes e Schawlow (1955)]. 


2.3.2 Teoria semiclássica 


Esta teoria foi desenvolvida nos trabalhos de Foley 
(1946), Bloom e Margenau (1953), Mizushima (1946) e 
generalizada por Anderson (1949). 

Nesta teoria é usada, fundamentalmente, a classical path 
approximation, isto é, os graus de liberdade translacionais das 
moléculas emissoras e perturbadoras são tratadas classicamente. 
Assim, as coordenadas das moléculas colidentes são 
parametrizadas no tempo e a interação molecular pode ser 
tratada pela teoria das perturbações dependentes do tempo. Os 
graus internos moleculares são tratados quanticamente. 

Para o caso de gases pesados, como OCS, N5, 0,5, NHs, 
etc, a temperaturas altas, o tratamento clássico do espalhamento 
é razoável, pois o “range” da força molecular, tão importante no 
alargamento da linha no limite de validade da Teoria da Integral 
de Fourier, é muitas vezes maior do que o pacote de onda 
associado à molécula colidente, isto é, quando bmin >> À, em 
que bmin é a distância mínima de aproximação da colisão e À = 
h/p é o comprimento de onda da molécula colidente [Massey e 
Burhop (1952)]. 

Os trabalhos de Foley (1946), Bloom e Margenau 
(1953), Mizushima (1946) e Anderson (1949) foram 
apresentados de uma maneira muito clara por Tsao e Curnutte 
(1962), de que brevemente faremos um retrospecto. 

Com efeito, Tsao e Curnutte (1962), supondo trajetória 
clássica para as moléculas colidentes, assumem o seguinte 
Hamiltoniano: 

H= H, + Ho, +He(t) + HR + Hip, (2.3.2.1) 
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em que Hi; e Ha são os Hamiltonianos livres das moléculas 1 
(emissora) e 2 (perturbadora), He(t) é o Hamiltoniano da 
interação molecular colisional, Hr é o Hamiltoniano livre do 
campo de radiação incidente e Hir é o Hamiltoniano da 
interação entre a molécula emissora e o fóton incidente. 
Calculando quanticamente a probabilidade de transição na 
aproximação de dipolo e usando o formalismo da segunda 
quantização [Sakurai (1967)] obtém-se para o coeficiente de 
absorção — potência absorvida — a expressão: 


a = 20 Re | dt exp(-im,) Tr [p, d, Um (1) d, Un(0)), 
0 


(23,22) 
em que: 

« = frequência do fóton absorvido. 

Po = matriz densidade de estados não perturbados. 

d, = componente na direção do campo de radiação do 
momento de dípolo. 

Um = operador de evolução temporal associado ao 
Hamiltoniano: 


H=H,+H + Hd(t) 


UA nai 


TS mUm> U, (0) =1 


e que satisfaz a equação. 


A, = indica uma média sobre as várias trajetórias das 
moléculas perturbadoras. 
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Anderson (1949), em seu trabalho, considera as colisões 
adiabáticas — aquelas que provocam apenas mudança de fase — e 
as diabáticas ou fortes — aquelas que são capazes de provocar 
transições internas. O autor propõe um modelo para a colisão 
diabática (forte): esta ocorre quando a colisão se efetua abaixo 
de um certo parâmetro de impacto crítico - bo [Cattani (1968)]. 
Tsao e Curnutte (1962), ao levarem em consideração tanto as 
colisões ressonantes quanto as não-ressonantes, generalizaram o 
trabalho de Anderson (1949), pois este não havia considerado as 
não-ressonantes. Assim, segundo Tsao e Curnutte (1962), o 
alargamento e o deslocamento de uma linha espectral são a parte 
real e a imaginária da seção de choque: 


6 = | 21b S(b) db , (2.3.2.4) 
0 


em que b é o parâmetro de impacto e S(b) é uma soma sobre os 
números quânticos das moléculas emissora e perturbadora para 
cada b. A forma de S(b) pode ser vista, por exemplo, em Cattani 
(1968). 


Anderson (1949) utiliza a aproximação de impacto, 
segundo a qual o tempo de colisão é muito menor do que o 
intervalo entre colisões, que é considerado sempre entre duas 
moléculas; além disso, entre duas colisões sucessivas a molécula 
colidente é considerada livre. Ele ainda considera trajetórias 
retilíneas para as moléculas colidentes, justificando estas 
trajetórias por considerar fracas as colisões que se processam 


com parâmetro de impacto b > bo. 
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CAPÍTULO 3 


TEORIAS SEMICLÁSSICA E 
QUÂNTICA DO ALARGAMENTO E 
DO DESLOCAMENTO PROVOCADOS 
POR ELÉTRONS E ÍONS 
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3.1 Teoria Geral da Forma de Linhas Devido a 
Colisões 


Segundo a Eletrodinâmica Quântica [Tsao e Curnutte 
(1962)], um sistema quântico emite ou absorve radiação 
segundo a lei: 


P(0) = de n,(09) F(0), (3.1.19) 
Cc 
com: 
T 2 
Fo) = Hr | dt exp Go Gym Idi y;(0)| (3.1.1b) 
0 
em que: 


P(o) do é a potência emitida ou absorvida com fregiiência 
entre me o + do; 


<yit) | e <yi(t) | são os estados final e inicial do 
sistema (moléculas emissoras, colidentes e o fóton emitido); 


d é o momento de dipolo do sistema; 
1+ N(0), 


n,(0) = 
No) 


JA 


sendo N(w) o número de fótons com fregiiência entre o e q + 
da da radiação incidente. A primeira expressão ocorre para o 
caso de emissão e a segunda, para o caso de absorção. 

Ao tratar de um sistema com muitas partículas a uma 
dada temperatura T, a forma da linha F(w) é mais complicada 
que a expressão (3.1.1b), pois temos de somar sobre todos os 
estados finais, e fazer uma média sobre os estados iniciais das 
moléculas emissoras, fazer uma média sobre os estados internos 
e sobre todos os possíveis movimentos das moléculas 
perturbadoras (estas duas últimas médias indicaremos por Av.). 

Quando a média Av. indicada acima é calculada 
classicamente (os graus de liberdade translacionais das moléculas 
emissoras e perturbadoras são tratados classicamente, isto é, elas 
descrevem trajetórias retilíneas), obtemos o resultado da Teoria 
Semiclássica da Integral de Fourier (TSCIF), desenvolvida 
fundamentalmente nos trabalhos de Lindholm (1942), Foley 
(1946) e generalizada por Anderson (1949), ao incluir as 
colisões fortes (diabáticas). Aliás, um excelente artigo de revisão 
sobre a TSCIF foi apresentado por Ch'en e Takeo (1957). 

Como eram insatisfatórias as previsões de deslocamento 
de linhas, feitas na teoria de Anderson (1949), um tratamento 
completamente quântico do problema se tornava necessário, o 
que foi feito por Baranger (1958), Kolb e Griem (1958), Fano 
(1963) e Di Giacomo (1964), calculando quanticamente a média 
Av. Porém, nenhum deles considera explicitamente o tipo de 
interação (dipolo-dipolo, dipolo-quadrupolo, etc) entre as 
moléculas colidentes e emissoras, como também não são 
consideradas as colisões fortes. 

Em vista do exposto acima, Cattani (1968a; 1970a-d) 
usou um formalismo quântico no qual considerou explicitamente 
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as interações moleculares, expandiu a matriz S de espalhamento 
até segunda ordem de Born e tratou adequadamente as colisões 
“fortes” e “fracas”. Desse modo, mostrou como poderiam ser 
obtidos o deslocamento e o alargamento de linhas espectrais na 
região de microondas, infravermelho longínquo, infravermelho 
próximo, e óptica. 


Agora passamos a expor a Teoria Quântica baseados em 
Baranger (1958) e Cattani (1968a, 1970a-d). Vamos partir da 
expressão (3.1.1b) e reescrevê-la na seguinte forma: 


El E 
F(o) = Tr aj dt exp [io(t-t) | ó(t.t) . (3.1.2a) 


em que: 
ó(tt) = » pa (il) ldhye(t)) (us(gtánys(0)| » 6.1.2) 


pois d é um operador Hermitiano. A função d(t,t) é chamada 
Time Correlaction Function (“Função de Correlação 
Temporal”). 


Para uma linha isolada, bastará considerar os estados 
inicial e final apropriados sem fazer a soma em if. No caso de 
linhas que se superpõem, fazemos a soma em if sobre estes 
estados. 

Definindo s = t — t? e levando em conta que os tempos 
de interação intermoleculares são muito curtos comparados com 
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os tempos de observação (T), as expressões (3.1.2a,b) ficarão 
(Baranger, 1958): 


oo 


F(0) = mo exp (ias) O(s) ds (3.1.3a) 
oii 
o(s) = E pal (ys(0) tdi (0). (yes (Idys(s)) | (313) 
Introduzindo agora o operador evolução T(s) definido 
DOE 


Te) y(0)> = | y(a)>, (3.1.4) 


e considerando que [Cattani (1968b)]: 


> [y:(0) (y;(O) [=1, PRE dp: 
f 
então: 


0 (s) = Tr já T*(s) d T(s) p| (3.1.5) 


Av 
em que (+) indica Hermitiano conjugado. 


Como F(m) é real, então d(-s) = 0*(s), em que (*) indica 
complexo conjugado, e portanto: 
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Fo) = a Re (es (icos) O(s) os (3.1.6) 
0 


Vê-se pela expressão (3.1.6) que a forma da linha 
espectral é obtida através da transformada de Fourier da 
Função de Correlação Dipolo-Dipolo. 


Para calcularmos a função de correlação O(s), vamos 
considerar o caso de emissão, ou seja, aquele em que o estado de 
energia mais baixo é o estado final. Como em geral o estado de 
menor energia é menos polarizável (menos perturbável), iremos 
tratar apenas com um estado perturbável, o inicial, ficando o 
estado final com energia bem definida (estacionária). Esta 
situação é denominada por Baranger (1958) de Caso de Um 
Estado (“One-State Case”). Ele trata também, em seu trabalho, 
o caso em que os dois estados são perturbados (“Two-State 
Case”). 

Assim, calculando explicitamente o traço indicado na 
expressão (3.1.5), resultará: 


álvoty dr ovstyldly,) tw |T(s) | vo)x 


ó69)=5 5 |(y, 
x(Vo olv, 1. (3.1.7) 


Como o estado final é bem definido, isto é: 
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T(9) | W;) = exp [-À E, 5) yo), (3.1.8) 
então: 
9(s) = Tr [DT(S) pJ,,- (3.1.9a) 
em que: 
De E Ly) exp [= E, ] (wi | | (3.1.9) 


Se o estado final for degenerado, e se a energia do estado 
final — Es — é tomada como origem das energias, virá: 


D=5 div). (Wild. (3.1.10) 


Usando a aproximação de impacto em que a matriz 
densidade é independente do tempo [Baranger (1958), Cattani, 
(1968b)], teremos: 

4 (s) = Tr[DT(S) p). (3.1.11) 

Calculemos, portanto, Tav(s). Inicialmente, tomemos a 


equação de Schrôdinger: 


ER by(s) = Hhy(s) = (E, + Vr(9) hyts)) , (3.1.12) 


al 


em que Vr é o Hamiltoniano de interação molécula emissora- 
perturbadora e H, o da molécula não-perturbadora. 
Sendo: 


T() hy(O) =) > ih ER T(s) = [H,+Vr(9)] T(S) . 


(3.1.13a,b) 
Usando a representação de interação, ou seja: 
T(s) = exp [- iH, 5] Us) , (3.1.14) 
teremos: 
h EN U(s) = Vy(s) U(S) , (3.1.15a) 
em que: 
Vi(s) = exp LH, s] Vr(s) exp [- LH, 5) (3.1.15b) 


Formalmente, a solução desta equação diferencial 
[Bloom e Margenau (1953), Merzbacher (1961), Davydov 
(1968)] é dada por: 


U(s) = — à [vio U(t) dt=3 exp - n Í NATO a 
0 0 


(3.1.16) 
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em que 3 é o operador de Dyson ou operador de 
ordenação temporal [Dyson (1949)]. A interação Vir (t) é uma 
soma de interações Vi(t), Vy(t). ... .VN(t) devido às N 
moléculas perturbadoras. 

Portanto: 


U(s) = 3 ef ú Í [Vi(D SAO vio] , 
0 
(31.17) 


então: 


N 
Uav(s) = 3 deseiivo a) | (3.1.18) 
0 Av 


A expressão (3.1.18) é extremamente complicada, pois 
se trata de calcular a potência N-ésima de um operador, tendo 
ainda um “entanglement” [Feynman (1951)] produzido pelo 
operador de Dyson. Portanto, para realizar essa potência, 
Baranger (1958) usa a aproximação de impacto, que veremos 
agora com detalhes. 


a) o tempo livre médio entre colisões — s — é muito maior 
do que o tempo de colisão — Tc — isto é: s >> T;; 


b) a molécula emissora colide com as perturbadoras 
obedecendo a uma sequência temporal e com uma de cada vez, 
isto é, através de colisões binárias. 
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c) na maior parte do tempo a molécula emissora é 
tomada como livre, ou seja, sem interagir com as outras. 
Durante a colisão sua função de onda é perturbada pela colisão 
e é determinada pela matriz — S — de colisão. Desse modo: 


yi)=S, by); ny?) =8, ny) = (SS) Iv), 


e assim sucessivamente, em que hy;) é a função de onda da 


molécula i-ésima antes da primeira colisão e Iyi) sua função de 


onda após a primeira colisão, etc.. 


d) a mudança na função de onda ocorre de um modo 
substancial após um número muito grande de colisões. 


e) a energia cinética das moléculas perturbadoras deve 
ser muito maior do que as perturbações causadas pelas colisões. 


Pode-se mostrar que [Baranger (1958), Cattani (1968a)] 
para estar satisfeita a aproximação de impacto devemos ter: 


o? <<1/n2 ou ovT<I1/n?, 


em que O é a secção de choque de colisão, n(2) a densidade de 
moléculas perturbadoras, v a velocidade relativa média entre as 
moléculas colidentes e Tc o tempo de colisão. 


Usando-se a aproximação de impacto referida acima, a 
expressão (3.1.18) torna-se: 


40 


N 
Uav(s) = E esc iivo a) | (3.1.19) 
0 Av 


Considerando: 


S;0) = ef 7 Í V;(t) a (3.1.20) 
0 


e a matriz de colisão — matriz S — para a j-ésima colisão 
ocorrendo no instante t = 0, então: 


Users Sto tas) sa S(to) S(tj)|,,» (3.1.21) 


em que ti =0O,..,tn=s e N é o número de colisões 
ocorrendo entre O e s. 


Se a colisão no intervalo O — s ocorrer no instante t, 
teremos, usando a representação de Heisenberg: 


SO = es] E Í H, [5/0 ef =)H, ] (ÃO 
0 


O —, 0 


comO<t<s. 


Como na aproximação de impacto a grande maioria das 


4Z 


colisões é “fraca”, podemos escrever que: 


S(D=1+f (3.1.23a) 
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em que fi, que é o termo devido à colisão, é muito menor que 1 
(f, <<1). Assim: 


= exp [é ÍH, ] S,(0) exp - à im, dh (3123) 
0 0 


Classicamente, a probabilidade de uma colisão com um 
único perturbador no intervalo representado por te t + dt é: 


Ni dtdv, 


em que N é o número de moléculas perturbadoras do sistema e 
V indica todos os possíveis tipos de colisões, isto é, todos os 
possíveis valores para parâmetro de impacto, energias 
translacionais das moléculas colidentes e estados internos das 
perturbadoras. 

Quanticamente, levaríamos em conta todos os possíveis 
estados translacionais das moléculas colidentes e todos os 
possíveis estados internos das perturbadoras. 


A matriz S para todos os tipos de colisão e para todos os 
instantes, no intervalo O > s, será dada por: 


exp - i Í V'(t) a = 
0 


Av 
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s 


E ! E If NT dt a exp =H, ) S(0) exp - ! H, dj 


| Ni fd exp (im t)[[(S- av] eo (im) 


Chamando: 
E ad 

H= sa] (1-S)dv, (3.1.24a) 
Hi s-e [ÁH ] Hip Bm, io (3.1.24b) 
teremos: 

1 DAVA) a cp dE —1 f ' 
Eos prai] = iN H(9 a) (3.1.25) 
e, portanto: 


o 


N 
U)=9[-ân ] H(y) a (3.1.26) 
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N 
Quando N >> 1, exp(x) = [1 + a) , então a expressão 


(3.1.26), torna-se: 


U,(9)= 3 o é Í Ho a) , (3.1.27) 


Tav(s) = 3 exp 5 [ (Ho + H(0) a) - (3128) 


o —, 


Se: 


Ts) = exp [ii + H)s| => 


jh a Tuw(S) = (H, +H) Tu(S) (3.1.29a,b) 


cuja solução é [Merzbacher (1961), Davydov 1968)]: 


Tav(s) = 3 exp pá 


(o) 


[Ho + H'(0] a) » (3.130) 


portanto: 


Tav(S)= exp [- z (Ho + H) | (3.1.30b) 
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Obtido Tav(s) , vamos agora calcular o perfil da linha 
para o caso de um estado. Segundo as expressões (3.1.5) e 
(3.1.6), temos: 


F(0) = À Re [exp (1005) Tr[D Ta (9) p] ds. (3.131) 
0 
Como consegiiência da aproximação de impacto 


[Cattani (1968b)]: 


p(s) = Po 


então: 


tT po! F(wo) = Re [ expúos) Tr [D exp E E (Ho + H) y| ds. 
0 


(3.1.32) 


Lembrando que D, Ho e H não dependem de s, a 
integração indicada acima resultará: 


aTp;'F(w)= Re| TrfiD(o- H,lh-H ny = 


=-ImTr| D(0-H,/h-H/ |. 


(3.1.33) 
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Introduzindo agora os auto-estados do operador Ho + 


H,ou seja: 
1 : 
(Ho + B) ld) = (mi wi) ld), 


com wi; > O . Como H não é Hermitiano, então os estados | 0) 
não são ortogonais. Desse modo, vamos então introduzir um 


outro conjunto de estados | y;) tal que: 


(w; 6) = dj » com 2a | 0) (ui! = 


Assim, a expressão (3.1.33) torna-se: 


aT o; F(o)=-Im> (y, |D(wo-H,/h-H/ hn |6)= 


=-Imb (0-0, +iw) (y,|Dló)= 


(0-0, —iw,) 
Dld.). 1.34 
A io E TWW, |9,) (3.1.34) 
Fazendo: 


My! DIG)=A+iB, 


46 


teremos: 


V. 
“F(9)=5———Re(y,|Dl|6,) - 
ATP F(o) Dr e(y,|D|9,) 
= Dip: (3.1.35) 


2 2 
(0-0) +w; 


Vê-se pela expressão (3.1.35) que o perfil de linhas 
superpostas (“overlapping lines”) não é composto somente de 
linhas Lorentzianas. Há, também, alguns termos assimétricos. 

O trabalho de Baranger (1958) nos dá apenas a forma 
analítica da linha. Ele chega aos resultados de Anderson (1949) 
e aos da Teoria da Integral de Fourier efetuando algumas 
simplificações em seu resultado formal. Vejamos então como 
fazer isso. 


Tomemos uma representação em que o operador Ho é 
diagonal, isto é: 


Ho ly) = E ly. 


No caso em que podemos considerar H 
perturbativamente e levar em conta somente o termo em 
primeira ordem, obtemos os resultados de Anderson (1949) para 
o caso de não-superposição: 


(ya E Ly =(yil HI) dia. 
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Esta consideração impõe que (y;l H | w;) seja pequeno 
comparado com o espaçamento entre os níveis atômicos. 
Portanto, valendo apenas para o caso de linha isolada, teremos 
para a expressão (3.1.35), uma vez que D é um operador 
Hermitiano e portanto Im (y;l D | wy;) = O, a seguinte forma: 


TTp; Ho) = 5 [E (ID Iv). (3.136) 


Sendo: 


(0-1w)=E+(ulI HIy)=E+Regyl HIy)+ 
+iIm6yi HI yj, 
teremos: 
o=E+RetyilH Iy)p=E+S > 
Si=RetyilH | yj, (3.1.37a) 
w=-Imtyl EH Iwy. (3.1.37b) 


Portanto, a forma do espectro dada pela expressão 
(3.1.36) torna-se: 


Fw) = 


ID Iv), (3.138) 
PRATO E Wi vo), « 
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em que Si é o deslocamento e wi; a largura da linha 
isolada e Lorentziana. Esse resultado foi obtido por Anderson 
(1949). 


Na teoria das colisões podemos sempre escrever: 
(ul S ly) = exp(-idi), 
em que 0; é o deslocamento de fase (“phase shift”) da função 
de onda e o; < 1 caracteriza o tipo de colisão. Então, usando-se 


as expressões (3.1.24a), (3.1.37a,b), a expressão acima e 
considerando a ortogonalidade de (y.|, 


btt) = (ut = 0-S)avty) = 
=(i[ (wii S ny)av)-i fav= 
= [1] a, exp(rióy) dv | — À fav = 
=| | 04 (cos 6,-i sen 6) dv] — i fav. 


w=-ImGyl HH ly) = 7! (1 - a; cos 4) dv, 
(3.1.39a) 


S=Re(ylH Iy)= 23 o, sen 4, dv. (3.1.139b) 
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Se todas as colisões são elásticas, isto é: o; = 1, teremos: 
S;=]sen Odv e w=](I-cos O) dv, 
(3.1.40a,b) 


que é o resultado da Teoria da Integral de Fourier [Lindholm 
(1942), Foley (1946), etc...]. 


Por outro lado, se todas as colisões são inelásticas, isto é: 
oi = 0, ou se di é muito grande, o que implica ser (yil S hyi) = 
oi exp(-idi) = 0, o deslocamento Si é nulo; este é o resultado 
obtido por Anderson (1949). A largura wi é igual à fregiiência 
de colisão, que é o resultado obtido por Lorentz (1906). 

A aproximação de impacto se aplica quando o 
alargamento for muito menor do que a energia cinética média 
das moléculas A mesma relação é válida para o deslocamento 
[Baranger (1958), Cattani (1968a)]. 


Di Giacomo (1964) calculou a forma de linhas espectrais 
partindo da hipótese de que o coeficiente de absorção ” da 
radiação incidente formado por partículas com momento de 
dipolo elétrico d (ou magnético) é dado por: 


7 =(87/cE) (EM 1 dr), 


em que f...) indica uma média no tempo, 
E(t)= E, cos(wt — k-Wéo campo elétrico (magnético) incidente 
eM (t) é o momento de dipolo elétrico (magnético) do sistema 
formado pelas moléculas absorvedoras (emissoras) por unidade de 
volume. Ele é dado por M (t) =n0 « d(t) >= nºTr(p(md) , em 
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que nº é o número de moléculas emissoras (absorvedoras) por 
unidade de volume e o(t)é a matriz densidade do sistema. A 
evolução do estado do sistema perturbado pelo campo incidente é 
determinado usando a equação de evolução temporal da matriz 
densidade dada por i(h/27)090/dt=[H,p(t)], em que Hé o 
Hamiltoniano total do sistema formado por moléculas emissoras 
(1) e perturbadoras (2). 

Desse modo, efetuando um cálculo relativamente 
extenso, Di Giacomo (1964) mostrou que a largura AV e o 
deslocamento S de uma linha isolada e Lorentziana são dados 


por: 


lf = 1 = 
Av=—Re(H.) e S=—Im(H,), 
27 (E) 27 (E) 


(3.1.41a,b) 
em que: 


= n2p?(B/27m a 
H = (6, ) Air bo 2, Cy x 


: Tr [exp(-8nº)] M; Mr 


x Tr? [exp(-8n)) | dg exp(-Ba”/ 2m) CREA la) 
(3.1.410) 


nl) é o número de moléculas perturbadoras por unidade de 


volume, h() é o operador de energia interna da molécula 
perturbadora, e: 
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o ANTA 
A, = Ste) (3.1.41d) 


[IADE +D]" 


Cum, =C(5:13,: M[0Ms) . C(S1J,: M0M;), (3.141€) 


com Tr(2) indicando o traço que envolve os graus de liberdade 
interno das moléculas perturbadoras, m a massa reduzida das 
moléculas emissoras e perturbadoras, q o momentum linear 
relativo, 


Lym, = 200 I,M, 


ll 4 


Rjo,1,M;) dum, 200, ],M,|Rº 


UI 4 


QJM,)ô 


MM, 


47" (a, ],M,|RÓ(E,-H,) 


GSM, MOI MR jo] ;M,), 
(3.1.41f) 


e R vem da matriz espalhamento S , isto é: S= 1 — 2miR 6 (E-H5). 


3.2 Alargamento e Deslocamento Devido às Colisões 
Eletrônicas 


A expressão (3.1.36) dá a forma da linha no caso em 
que as colisões são tratadas segundo a aproximação de impacto. 
Ela se aplica às linhas isoladas com apenas um estado 
perturbado. Griem, Kolb e Shen (1959) e Griem, Baranger, Kolb 
e Oertel (GBKO) (1962) mostraram que esse é o caso de todas 
as linhas do Hélio neutro, com exceção da linha 2P — 4ºD de 


2 


o 
4471 A [Berg, Ali, Lincke e Griem (1962)]. Quando as 
partículas perturbadoras são elétrons, o cálculo da média Av. se 
torna mais simples, pois os elétrons, nesse contexto, não têm 
graus de liberdade internos e não são consideradas mudanças de 
spin. Assim, o número de colisões com velocidade relativa entre 
ve v + dv eo parâmetro de impacto entre p e p+dp 
podem ser obtidos da Teoria Cinética dos Gases [Ch'en e Takeo 
(1957)]. 


dv =Nvf(v)dv.27pdp 


em que f (v) é a distribuição de velocidade de Maxwell- 
Boltzmann, e N é a densidade de elétrons. 

Nesse caso, as expressões para a semilargura (Av) e para 
o deslocamento (S) (em unidades de frequência linear) de uma 
linha isolada e Lorentziana, em estados de energia não 
degenerados | y >, tornam-se: 


Av = 5 Re (E), (3.2.1) 
S = 5 Im (H), (320) 


com: 


H=Nº [vi(mdav|2mpapf<w I-SIw >), 
(3.2.3) 


em que ( ) indica uma média sobre todas as possíveis direções 


do elétron incidente. A expressão de H para o caso em que os 
dois estados são perturbados pode ser vista em Baranger (1958) 
e Cattani (1972, 1973). 
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A matriz S é geralmente calculada [Griem, Kolb e Shen 
(1959), Lewis (1961), GBKO (1962), Griem (1962a)] usando a 
interação dipolar em segunda ordem de Born e trajetória 
retilínea para os elétrons. Assim: 


<yS -1ly,>= + fat<yr VCD Iyç>— 


+ ET di 2 de, expf (gt + 0nto)] X 


x <yri V(t) lyn>.<WynlV(t) Iyr>, (3.2.4) 
com: 
Rd a (3.2.4b) 
É + v? t A 


em que rf é o vetor posição do elétron atômico; pevo 
parâmetro de impacto e a velocidade, respectivamente, do 
elétron colidente. 

Sendo o estado | y; > de paridade bem definida, o 
primeiro termo da expressão (3.2.4a) é nulo. No cálculo do 
segundo termo, a média em todas as direções das componentes 
de v ede p e termos do tipo psp, também é nula, bem como 
termos da forma vxvy não contribuem à integral, pois eles estão 
multiplicados por funções ímpares de t. Para uma distribuição 


isotrópica de elétrons, os termos fp? e [vi (1 = X,y,Z) 


1 


podem ser substituídos por p2/3 e v2/3, respectivamente. Com 
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essas simplificações, e chamando zm = wmp/v e x = vt/p, 
teremos: 


2 2 
(<yIS —1hy, 3) > evilrshy,>.<Wyilrohy, >) x 


x[ dx; [ dx, expli FA (x) ; (2.2.5) 
“o oo (1+x7) , (1+x5) 


Fazendo a integração em x [GBKO (1962)] e levando o 
resultado em (3.2.3), resultará: 


2 
a=o4m(h) d <Vilry, > 
es ain 

x [a(z5) + ib(z5)]. (3.2.6) 
em que: 
E sad TAC) 
a, = — = raio de Bohr; a(z) =| dz"; 
e 


GA 
V4 


b()=[ dr ; (o) =2[Kiz)+Ki(z)] ; 


ate) 


Z? 


B(x)= ZE P | dz! , 
0 


Z 
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sendo Ko (|z|) e Ki (|z|) funções modificadas de segunda 
espécie de Bessel [Watson (1966)]. 


Para considerar o efeito das colisões fortes, usaremos o 
modelo de Anderson (1949). Assim: 


H e His + Mortes , 
em que: 
Ho: ÉS t2(2) np. po o fed) vs 
ii 3 lm v gds a, 
X La(a)+ib(zp")), (3.2.7) 
Hines = —N [v dy f(V) 7 Pain) (3.2.8) 


com: 


representando o parâmetro de impacto mínimo ou raio de 
colisão forte. 


A determinação de Pmin é uma questão muito delicada, 
pois esse “raio” é definido de tal modo que, para p < Pmin, O 
tratamento perturbativo deixa de ser válido. As aproximações 
usadas por GBKO (1962) e Griem (1962b), em geral [Cattani 
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(1973)], não são muito rigorosas e provavelmente constituem as 
principais fontes de erro nos cálculos semiclássicos. 
No limite de altas temperaturas (altas velocidades dos 


elétrons), GBKO (1962) mostraram que o Ha definido pela 
equação (3.2.7) torna-se: 


E ax (hn n(d 
EE RA = Rs (2) RA > » | <Wiliçhy, > P x 


x (ag? [enfizg? ei , (3.2.8) 


em que os sinais + são os de Om e Pmin é calculado por 
intermédio da seguinte relação: 


II 
WIN 
3 
< 

o | 
ER 
B 
Sn 2 
18) 
DN 
fe 
O 
b 
XxX 


E <yy S(Pmin) — 1Iyy> | = 


2 
XE Lt <wilr > E = 3( h Jat=r. (3.2.9) 


MVP min 


em que a é o número quântico principal do estado ly; >. 


No caso de baixas temperaturas (baixas velocidades 
eletrônicas), as colisões são adiabáticas e o tratamento 


perturbativo é válido para qualquer parâmetro de impacto p. 
Nesta situação, H é dado por: 


SL 


NR 


BU nt? 


2 323 lm 


Om 


sos 
x z z ema (14103). (3.2.10) 


4Z 


No caso intermediário, isto é, para temperaturas nem 
muito altas nem muito baixas (que ocorre na maioria dos casos 
analisados), uma das possíveis equações que serve para 
determinar Pmin , pois ela satisfaz as condições limites de altas e 
baixas temperaturas, é dada por: 


2 
o h -2 
| <yil SidPmin) — 1 hy> = 4 do] (ao) x 


A [z > | < Vilts UA > E Ae) + 
ny? 
+ Di<wiçh > PB(zp") ] = 


3/2 


=(3/9)", (32.11) 


sendo agora H dado por: 
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2 
H = — N | dy f(v) [Eros Ea (2) da 


xDD i<yshy, > À falei) +i biz e) . GDA 


Uma outra alternativa (que é menos usada) para a 
determinação de Pmin , que não escreveremos aqui, pode ser vista 
no artigo do GBKO (1962) [vide expressão (3.17), pág. 182]. 

No caso de altas densidades e de níveis pouco afastados, 
a expressão (3.2.12) tenderá a superestimar a largura e o 
deslocamento porque não foi considerado o efeito de screening 
de Debye. GBKO (1962) sugerem como calcular 
aproximadamente apenas a largura da linha usando um resultado 
obtido no cálculo das linhas do hidrogênio (H). Assim: 


A, =—7N [vt dy [Plot Ea “(oo a) 
(3.3.13) 


em que Pmin é calculado pela expressão (3.2.9), pp é o raio de 
Debye e t, é o raio vetor do elétron atômico. 


A expressão (3.2.13) é calculada dividindo-se as colisões 
consideradas [Griem, Kolb e Shen (1959)] em três tipos: fortes 


(p < Pmin ), fracas (Pmin < p <pp)e “screened” (p > pp). 


O limite mínimo de N para a aplicabilidade de (3.2.13) é 
definido por: 
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pp=V/O - (3.2.14) 


Usando-se as expressões que definem o raio de Debye e 
a velocidade média (V), tem-se: 


Nin E Mm 05 (27º e). (3.2.15) 


min 


3.3 Alargamento e Deslocamento de Linhas Isoladas 
Devido a Ions 


A Teoria Estatística ou Estática desenvolvida principal- 
mente por Holtsmark (1919), além de Margenau (1932; 
1933a,b;1935; 1951), Kulp (1932; 1933), Kuhn (1934) e Kuhn e 
London (1934), é frequentemente usada para tratar a perturbação 
provocada pelos íons. Porém, como próximo do centro da linha 
essa teoria falha [Griem, Kolb e Shen (1959)], utiliza-se então 
uma aproximação “dinâmica” na qual são feitas duas hipóteses: 
(1) as perturbações dos íons não se superpõem no tempo, sendo, 
portanto, escalarmente aditivas, e (2) as colisões dos íons com 
os átomos são adiabáticas. Como usualmente as velocidades dos 
íons são baixas, a condição de adiabaticidade está satisfeita; por 
outro lado, a aditividade das perturbações dos íons é válida 
porque estas são pequenas, o que nos permite considerar as 
colisões como sendo binárias (aproximação de impacto). 

Deste modo, o perfil Pions() devido aos íons (omitindo 
o fator numérico e o elemento de matriz de dipolo) é dado por: 
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Bons(0) = Re ] dt exp(iAot) CON (3.3.1a) 
0 


em que: 


O = ep psi e(O = exp - > [s9,(t) a (3.3.1b) 
FI 


Flo 


sendo Ni a densidade de íons e P;j(t) = Am;(t”) a mudança de fase 
provocada pelo j-ésimo íon. Usando a teoria de perturbações 
independente do tempo, o Aw, é dado por (fazendo uma média 


sobre os números quânticos magnéticos): 


2 
-l( h l 2 — 280 
amoalçÃ) E Eae tem > é =, 


(3.3.2) 


que nada mais é do que o efeito Stark Quadrático. 
Admitindo trajetórias retilíneas para os íons, isto é: 
2 
r; =p+v(t-t) a 


teremos [Smimov (1969)]: 


P(D = [not dt' = 21C Í [ei (e-0]] de = 
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ro (3.3.3) 


A média indicada na expressão (3.3.1a) é calculada 
exatamente como no parágrafo anterior. Admitindo que os íons 
deslocam-se com uma velocidade média v , virá: 


josf 3 no = lexp[-i O 


H 


—co o 


= exj2rva ja [ pdpfexpl-iP, 1-9 = exp| (7.0) | 


(3.3.4) 


Portanto, a forma da linha será: 


Pons(0) = Re Í exp[ i AQt + g(LV.C) | af (3.3.5) 
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3.4 Alargamento e Deslocamento de Linhas Isoladas 
Devido a Colisões de Elétrons e Ions 


Considerando as perturbações provocadas pelos elétrons 
e íons como independentes, os resultados dos itens 3.2 e 3.3 
permitem escrever, para linhas isoladas e com apenas um estado 
perturbado, o seguinte perfil: 


P(w) = Re | dtexp[ i (o-d)t — wt + g(t7,C) | dt, 
(3.4.1) 
em que w e d são, respectivamente, a semilargura e o 
deslocamento da linha devido a colisões eletrônicas [esta 
notação é a de GBKO (1962)]. 


Vamos agora introduzir as seguintes variáveis: 


r=P,;/Pn:0=wt;T=Wi;0C=WPn/V; 


4 
WPm (o) nº 


2nC) (4 RAR o já 

|2nC| (4x 2 

a=| ] [8 IN; E [E] > I<yylro y, > x T | : 
(3.4.2a-g) 


em que pm é a distância média entre os íons. 
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Em função dessas novas variáveis, a expressão (3.4.1) 
tomará a forma: 


j(x,0,0) = Re | dtexp[ixt— 1 + g(,0,0)], (3.4.3) 


em que [usar (3.3.3) e (3.3.4)]: 


oo 


ame f 43 O a(lt-8 
s(TO,O) = E [rár [ao [exp [ ia E (te fre )+ 


isa Ô EE (t-9dro ôro E 
S (ro) + (7 — 3” (ro) + a] É Ri 


Griem (1962b) mostrou que AVtotal € Stotal conseguidos 
através do perfil reduzido j(x,0,,0) dado pela expressão (3.4.3) 
são essencialmente os mesmos obtidos com uma aproximação 
estática que leva em conta a correlação entre os íons e o 
screening de Debye produzido por elétrons [Mozer e Baranger 
(1960)]. A largura e o deslocamento totais são dados 
aproximadamente por: 


AViotal E W+1.75 0 (1- 0.75) w (3.4.5) 


Sotal=d + 200(1-0.75)w, (3.4.6) 
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em que : 
R = pa/pp= 68. 16. Nº (Ka T/ 2, (347) 


e q é dado pela expressão (3.4.2g). O parâmetro & para diversas 
transições do hélio (He) está tabelado no artigo de GBKO 
(1962) em função de Ni; e T. As equações (3.4.5) e (3.4.6) são 
válidas para wu <0,5eR< 0,8. 


Ainda é uma questão aberta se certas discrepâncias entre 
perfis de linhas teóricas e experimentais (especialmente no 
centro de linhas de Balmer de componentes interditas de hélio 
neutro) são devidas ao tratamento estático dos íons. Segre e 
Voslamber (1974), aplicando a teoria unificada de Voslamber 
(1972), não encontraram nenhuma diferença apreciável entre o 
tratamento estático e o dinâmico [GBKO (1962) e Griem 
(1962b)]. 
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CAPÍTULO 4 


TEORIA QUÂNTICA DO ALARGAMENTO 
E DO DESLOCAMENTO DE LINHAS 
ESPECTRAIS ATÔMICAS 
PROVOCADOS POR COLISÕES 
DE ELÉTRONS RÁPIDOS 
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4.1 Introdução 


Vimos no Capítulo 3 como se calcula o alargamento e o 
deslocamento de linhas espectrais atômicas, devido a colisões 
eletrônicas, quando as trajetórias dos elétrons perturbadores são 
consideradas classicamente. 

Veremos agora como calcular, com um formalismo 
puramente quântico, o alargamento e o deslocamento de uma 
linha atômica provocados por colisões eletrônicas. Assim, 
segundo Bassalo e Cattani (1972, 1974) e Bassalo (1973), a semi- 
largura Ave e o deslocamento Se de uma linha isolada 
Lorentziana são dados, no caso de estados de energia não 
degenerados (hy >), por: 


Av, O TO (4.1.1) 
2m 
1 ui 
S,=-— Im(H ,,.), (4.1.2) 
2m 


em que os índices I e F se referem aos estados inicial e final, 
respectivamente, e: 


sendo N a densidade de elétrons perturbadores, B = 1/(KpT), 


Kg é a constante de Boltzmann, T a temperatura absoluta, m a 
massa reduzida do elétron e da molécula emissora, q o 


momentum linear relativo, 


Ee 


Lr=2mi(yil R yo) -—2mi (ye Rº Ipe) - 


-47(yil R SE E9) yo (yrl Rº Ivo), 
(4.1.4) 


em que R vem da matriz de espalhamento S = 1 -2 mi (E — Ho). 

Em muitos casos, os estados | yr) são escritos como 
combinação de estados de energia degenerados | «J M):;Jéo 
momento angular total, M sua componente z e o são outros 
números quânticos. 

Quando os estados | « J] M) são apenas degenerados em 
M, Lir pode ser escrita, para transições | oJ) > | arJe) 
[Cattani (1970)], como: 


Ly = Ar) bo Cm f2ui < IM jIRIo)M, > ÔMMET 
My MF 
— 2mi < O) EMpIR“IoTEM, > ômpMET 


— 47º < 0 J;MjIRSE,-Ho)la);Mp> x 


x < 4);MEIR'IoSEM, > |, (4.1.5) 


em que: 


Ap = ITEM, + D(23, + DP? (4.1.6a) 


13 


Cmmp = COF; MOM5) . (13; Mj0M)) > (4.1.6b) 


em que Ar e Cm foram definidos no capítulo anterior, 
expressões (3.1.41d) e (3.1.41e), respectivamente. 


Trataremos, neste Capítulo, apenas do caso degenerado 
podendo o não-degenerado ser tratado de maneira similar. 


4.2 Cálculo da Largura e do Deslocamento 


O alargamento e o deslocamento de linhas atômicas 
produzidos por colisões eletrônicas são geralmente calculados 
[Griem, Kolb e Shen (1959), Lewis (1961) e Griem (1962a,b)] 
usando o método semiclássico do parâmetro de impacto em 
segunda ordem de Born ou usando o tratamento quântico em 
primeira ordem de Born [Kivel, Bloom e Margenau (1955) e 
Kivel (1955)]. 

Em 1970, Vainshtein e Vinogradov mostraram que 
melhores resultados para as secções de choque são obtidos 
quando usamos o tratamento quântico em vez do método do 
parâmetro de impacto, ambos em segunda ordem de Born. 

As limitações introduzidas pelo método do parâmetro de 
impacto no cálculo do alargamento e do deslocamento de linhas 
espectrais foram cuidadosamente analisadas no trabalho de Di 
Giacomo (1966). 

Moiseiwitsch e O'Brien (1970) demonstraram que as 
secções de choque elásticas obtidas usando o tratamento 
quântico em segunda ordem de Born, para energias do elétron 
diferentes de zero, são superiores às obtidas com a de primeira 
ordem de Born, considerando-se ou não o efeito de “exchange”. 

Assim, previsões mais precisas das formas de linhas 
espectrais poderiam, em princípio, ser obtidas expandindo-se a 
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matriz R na expressão (4.1.4) ou (4.15) até a segunda ordem de 
Born. Entretanto, é muito difícil dizer “a priori” quão boas 
seriam essas previsões, pois as limitações da aproximação de 
Born, em segunda ordem, não eram bem conhecidas [Holt e 
Moiseiwitsch (1968), Holt (1969)]. As teorias de formas de 
linhas, citadas acima, que existiam até a época em que 
desenvolvemos os nossos cálculos, usavam a primeira ordem de 
Born em um formalismo quântico, ou a segunda ordem de Born 
usando um formalismo semiclássico em que calculavam 
aproximadamente as colisões fortes por intermédio do método 
do parâmetro de impacto crítico [Anderson (1949), Griem, 
Kolb e Shen (1959), Griem, Baranger, Kolb e Oertel (GBKO) 
(1962) e Griem (1962a,b)]. Em nosso formalismo quântico não 
precisamos estimar os parâmetros de impacto críticos. As 
colisões fortes e fracas são obtidas de modo natural dentro de 
um tratamento completamente quântico de colisões. 

Desse modo, expandindo a matriz R na equação (4.1.5) 
até a segunda ordem de Born e substituindo este resultado na 


equação (4.1.3), Hp torna-se, desprezando o efeito de 
“exchange” entre os elétrons incidente e atômicos: 


o 3/2 
Hp =N (5) Ap Cymp X 


My MF 
2 
x [d'g exp e) (a | (2mi (0,],/M,| V I0,)Môymmp — 


— 21 (04) ;Mgl V 10);M5) ômmy + 


; 1 
+ 271 (ot V E,-Ho+ig V JM) ôMpMpT 


o 


o 


: 1 
— 211 (01Ma V EH EÉ V 0,3,M,) ôMMET 


di (1)My/l ki E,-H) lot] ,M) (Ot) AM Dé lot Mh)) ) É 
(4.2.1) 


em que V é o potencial de interação entre o elétron e o átomo, 
ionizado ou neutro. 

Bassalo e Cattani (1972) e Bassalo (1973) usaram para o 
potencial V a seguinte expressão: 


2 2 
ms ZE os cer. (4.2.2) 


em que R, e Y, são os vetores posição do elétron incidente e do 


a-ésimo elétron atômico, respectivamente, e Z é a carga 
nuclear. 

Entretanto, devido ao efeito dos íons presentes em um 
plasma, levaremos em conta o “screening” de Debye-Hiickel 
[Griem, Kolb e Shen (1959), GBKO (1962) e Griem (1962a,b)]. 
Assim , o potencial V dado pela expressão (4.2.2) foi substituído 
em nossos cálculos [Bassalo e Cattani (1974)] pelo potencial V” 
representado por: 


V'= Ze exp(-R/lp)/R — Se exp(+ R - El/Opk (IR — EO, 


(4.2.3a) 


em que /p é o raio de Debye-Hiickel dado por [Jancel e Kahan 
(1966)]: 
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lp = (KT Ne (1 + a”, (4.2.3b) 


em que KB, T e N foram definidos anteriormente, e z é a 
carga iônica dos íons no plasma. 

No caso em que o plasma contém somente uma espécie 
de íon de carga elétrica ze, a neutralidade do plasma é expressa 
pela relação Niz = N, em que Ni é a densidade iônica. 

Para a obtenção da largura e do deslocamento, 
seguiremos o trabalho de Bassalo e Cattani (1972) e Bassalo 
(1973), usando (4.2.3a,b) em vez de (4.2.2), na expressão 
(4.2.1). Desse modo, introduzindo os projetores: 


jagig>.<ãl e SlaM>.<aM, >, 


a expressão (4.2.1) tomará a forma: 


E (By a Ba 
Ho =Nh [5] App à A Cm o fa q exp [- se x 
x 2mi (04];M4l (GV 13) lo,3/My) dupmy o 


— 27 (04,3;Myl (6] V'|ã) 1o,3;M$) ômpmp + 


+2i dum, X | dl (0,3:M (giv'1ã) 0,3 M,) P x 


E 


— 27i ômime > Jog (00,3 ;My(dl v' da3M, x 


2.— 
Ra pa +ina|44 JA, - 
q — q + 2mA,, 2m 


= 47 [dg (2)My (IV IG) loJM,) x 


2-2 
x (07) My l(g| Vº ED x AS a | 


(4.2.4) 


em que Ask = Es,” Exu; (K = [ou F) é a diferença de 


energia entre os estados 


nn 


JS ) e [O J 6) » lembrando que foi 
usada a seguinte expressão [Davydov (1968)]: 


lim e = ix 809 + P(1/X). 


e>0 


Calculemos agora as contribuições dos termos de 
primeira e de segunda ordem. 
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4.2.1 Contribuições dos Termos de Primeira Ordem 


Lembrando que 1d) = exp (q-R/n)/h”, o termo de 
primeira ordem da equação (4.2.4) torna-se: 
N 


2nip (Wi — Vo), (4.2.1.1a) 


em que: 


Vi=[ PR VIR) = | d'R Po » Cymy (MM, IV” cam) 
I 


(42.1.1b) 


Ve=[dºR VER) =[dºR Pe > Cueme (O) Me VI est) 
Mr 


(42.1.10) 


com V” definido pela equação (4.2.3a). Vi e Vr são as 
interações eletrostáticas do elétron incidente, espalhado em todo 
o espaço, com o átomo no estado I ou F, respectivamente. 


Os potenciais de interação Vi(R) e Vr(R) são calculados 
exatamente apenas para átomos com um elétron. Para átomos 
com mais de um elétron eles só podem ser calculados 
aproximadamente. Moiseiwitsch e Smith (1968) dão referências 
de vários trabalhos, em que cálculos exatos e aproximados para 
V«(R) (K = I ou F) são realizados para átomos neutros ou 
ionizados. 
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Como Vr ou Vr são reais, os termos em primeira ordem 
em V” contribuem apenas para o deslocamento S. Portanto, 
usando as equações (4.1.2) e (4.2.3a,b) teremos: 


SO = A (Vi — Ve), (4.2.1.2) 


em que o índice (1) significa “primeira ordem”. 
4.2.2 Contribuições dos Termos de Segunda Ordem 


Tomemos inicialmente o primeiro dos termos de segunda 
ordem da equação (4.2.4). Introduzindo os projetores 


[a qiã> <ql e Sla)M,>.<a)M,, 
teremos: 
2( B ne js DE 
2m Nh [5] Ddr A Cum, JJ dé à dº q x 
Ba” na 2 
X exp —5m | < ÔwJM;j<qlV'Ig |>0,3,M,>| x 


2 —2 
q -q'+2mA,, 2m 


(4.2.2.1) 


Calculemos inicialmente: 
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<qIVIG> = RE [a R exp(i sã. R/h).[Ze? exp(-R/t5)/R — 


(4.2.2.2) 


em que Aq=G-q. Integrando em R (Apêndice C4-D, 
teremos: 


<qIV'Iãl> = É foroira Po Eemtaa-nm) 


(4.2.2.3) 


Tomando o eixo dos z na direção do vetor momentum 
transferido Aq , virá: 


4 
e 


xa) Mi<g | VIq>m])M,>P=——— >> x 
en (h/to)) + aq" | 
x IZ 3 — <0€];M,| 5 exp(iAg . 2,/h) lo,],M, >. 
(4.2.2.4) 


Definindo o fator de forma (“form factor”) 
Fim(Ag) por: 
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Fo (Ag) = Ar > Cum! Zô — 


M1 


—- <oJ;M;| > exp(iAg 2,/h) Lo,),M, >? . 


(4.2.2.5) 


Quando My = Mr = 0, o fator de forma F2 (Aq) assume 
uma forma mais simples (Apêndice C4-ID): 


2 
Fo (ag)-[2 dm - (04],/0] >) exp(iAq. 2,/h) L0,3,0> 


(4.2.2.6) 
Usando (4.2.2.5), a equação (4.2.2.1) torna-se: 


mi Nt(s o fade so Sa F(Ad) 


) loves) + aa] 


E) 
q—q +2mA,, 2m 


A integração em q e q dará (Apêndice C4-I]): 


(4.2.2.7) 


dias do 
211 Ne Es pa s tn (PB; A) —1 E w Win (B; E) 5 


(4.2.2.8a) 
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em que: 


oo 2 
Wi (BAm) = exp(2Y1) | A dx H (6) ãl x 


o (a +x 


x exp (va/0 =x). (4.2.2.8b) 


Sufita) = | A 6, (88) x). Dto) 


º (120) 
(4.2.2.80) 
com: 
1/2 
sd 5) po E) q: a =D 
y=(Yn/x)-x ; (4.2.2.9a-d) 


y 
D(y) = D(x,tn) = exp(=y”) x Erfi(y) = exp(-y”) x ] exp(t?) at 
(4229) 


em que D(y)= Integral de Dawson [Abramowitz e Segun 
(1969), Tricomi (1964)]. 


Como o segundo termo dos termos de segunda ordem na 
equação (4.2.4) é formalmente igual ao primeiro, ele poderá ser 
escrito na forma: 
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2m 


PRE 
— 2mi Ne E) a Sia(Bra pi E Wia(BiAm 
(4.2.2.10) 


Sendo elástico o terceiro termo da equação (4.2.4), ele 
simplesmente tomará a forma: 


Net (p')” 
E Es Es Wir(B;4=0) , (4.2.2.11) 
em que 
o 1/2 
Wir(B; 4=0) = | e so (82) E i exp(-x?) ? 


(4.2.2.12) 


com o fator de forma definido por: 


My MF 


1/2 
8 
K (82) d = Ap) D, Cum ÔMiME x 
x 2 — <omJ;M,| 3 exp(iAq z, /h) Mp x 


x E — <0J;Mel 3 exp(iAq z, /h) laJiMs + (42213) 
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Quando M; = Mr = O, Fr(Aq) assume uma forma mais 
simples, dada por (Apêndice C4-ID: 


Fe(Aq)= |Z — <o,J,01 3 exp(iAq z, /h) oã0> 


x |Z — <oyJ,0I 5 exp(iAg z, /hi) 1,05] = 


= | Fi(Aq) . Fela) JÊ. 


Portanto, a contribuição total de todos os termos de 
segunda ordem para o alargamento e o deslocamento, usando as 
equações (4.1.1), (4.1.2), (4.2.4), (422.89), (422.10) e 
(4.2.2.11), será: 


Av = Net (8) | [Wa (B; 4,0) + WB; Ao)]-2W (84-0) 
e 27 a nt In Fn > —Pn TEN ? 


2m 


(4.2.2.15) 


': 3 1/2 
se ERAS Ne E) > ER (BA) o Sen (B:Ap5) | 2 


2m 


em que o índice (2) significa “segunda ordem”. 


Como apenas os termos de segunda ordem contribuem 
para a largura de linha, a semilargura total Av. da linha é Av? 


e 2 


dada pela expressão (4.2.2.15). 
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Para obtermos o deslocamento total da linha, somamos 
Se. dado pela expressão (4.2.1.2), que dá as contribuições dos 


termos de primeira ordem, com so obtido pela expressão 


(4.2.2.17), que considera as contribuições dos termos de 
segunda ordem. 

Um aspecto que precisa ser salientado é que as 
expressões deduzidas neste item para obter as expressões gerais 
para a largura e o deslocamento reproduzem as expressões 
equivalentes obtidas por Bassalo e Cattani (1972) e Bassalo 
(1973), quando se considera apenas o potencial Coulombiano 
(4.2.2), ou seja, quando pp, >o (a = 0). Neste caso, as 


expressões (4.2.2.8b-c) transformam-se, respectivamente, em: 


oo y2 
Win(Bi4m) = exp(2Y1n) ] dx su (65) l d x exp (1/82 =) , 


x 
(4.2.2.17a) 
a 8 1/2 
Sm (B: Am) = | su (82) g D(x,y) 2 
sa 
e a expressão (4.2.2.12), em: 
pe 8 1/2 
Wir(B;4=0) = | E so 2) + exp(-x?) 
(4.2.2.18) 


Com o formalismo quântico referido acima, calculamos 
larguras e deslocamentos de várias linhas do íon He* (em A: 
3188, 3389, 3964, 4121, 4437, 4713, 5016, 5048, 5876) e 
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comparamos as nossas previsões com resultados experimentais 
[Bassalo (1975), Bassalo, Yamamoto e Cattani (1975), Bassalo e 
Cattani (1975), Bassalo e Cattani (1976a,b), Bassalo e Cattani 
(1977), Bassalo e Cattani (1978)]. No Apêndice C4-IV 
mostramos em detalhes somente os cálculos efetuados para 
estudar o caso da transição Psp (4713 ÂÀ) [Bassalo 
(1975), Bassalo, Yamamoto e Cattani (1975)). Os cálculos para 
outras transições são mostrados nos demais artigos. 

Comparando as nossas previsões com resultados 
experimentais para as linhas de He*, verificamos que havia bom 
acordo entre teoria e experiência para os deslocamentos S, mas 
as semilarguras AV eram cerca de duas vezes maiores do que os 
valores medidos. Por intermédio de uma análise cuidadosa dos 
nossos cálculos, verificamos que isso era causado 
essencialmente por duas razões: (1) a aproximação de segunda 
ordem de Born não era rigorosamente válida para estimar as 
colisões eletrônicas com os íons do He*:; (2) o uso das funções 
hidrogenóides para calcular os elementos de matriz das 
transições eletrônicas amplificava os erros introduzidos pela 
não-validade da aproximação de Born de segunda ordem. 
Diante desses fatos, tomamos a decisão de desenvolvermos dois 
novos modelos para calcularmos as formas de linhas. 

Um dos modelos, que denominamos Aproximação 
Dipolar, é mostrado na próxima Seção, 4.3. Ele é obtido 
diretamente do nosso formalismo quântico calculando os 
elementos de matriz das transições eletrônicas na aproximação 
dipolar. 

O outro modelo, mais elaborado e mais preciso, que 
utiliza um formalismo semiclássico, conhecido como Teoria 
Semiclássica  Convergente, será apresentado no 
Capítulo 5. 
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4.3 Aproximação Dipolar 


O cálculo do fator de forma indicado pelas expressões 
(4.2.2.13,14) é uma tarefa laboriosa e difícil, conforme 
verificamos nos diversos trabalhos que realizamos (relacionados 
no item anterior), nos quais calculamos a largura e o 
deslocamento de várias linhas do hélio neutro em um plasma 
quente. Além do mais, para calcular a secção de choque de uma 
colisão elétron-átomo na aproximação de Born em segunda 
ordem, é necessário realizar uma soma infinita [Holt e 
Moiseiwitsch (1968)] sobre todos os estados, discretos e 
contínuos, do átomo-alvo. Como essa soma é extremamente 
difícil de ser realizada, são usados métodos simplificados que 
complicam as interpretações das secções de choque obtidas 
[Holt e Moiseiwitsch (1968) e Holt (1969)]. No cálculo da 
forma da linha, não é necessário realizar essa soma sobre todos 
os estados do átomo emissor. Como podemos ver das expressões 
(4.2.2.8b,c) e (4.2.2.15,16), basta considerar somente os 
estados n que obedecem à condição Ex — En = Akn 2 — KBT. 
A energia transferida Akn é positiva se 2k (1 ou F) é um estado 
de energia maior que o estado n. Quando a transição virtual é do 
estado mais baixo de energia k ao estado mais alto n então Ak,n 
é negativo. Assim, o número de estados virtuais n que deve ser 
considerado torna-se menor à medida que k se aproxima do 
estado fundamental. Além de todas essas dificuldades, as nossas 
previsões para AV e S não estavam em bom acordo com os 
resultados experimentais, conforme registramos anteriormente. 


Assim, de acordo com estudos efetuados por Griem 
(1968, 1974), para muitos íons e átomos neutros, é muito mais 
simples e suficientemente preciso considerar o potencial de 
interação elétron-átomo (ou molécula) como uma interação 
dipolar. Desse modo, substituindo o potencial de interação V, 


88 


dado pela expressão (4.2.2), pelo potencial dipolar 
V=ed-R/R'(t), em que d=er é o momento de dipolo 
elétrico da molécula, na expressão (4.2.4), as expressões 
(4.2.2.15,16) passarão a ter as seguintes formas [Cattani (1983), 
Bassalo e Cattani (1993)]: 


SUR WX)+RZ, WXn)] 


n 


(4.3.1) 


n 


47º e! hº 2m Nº 
s,=- en | DR SG) Rã, SK) | 
(4.3.2) 


em que R: (j=I,F) são os elementos de matriz dipolar [Griem 
(1968)], e: 


Xa=EME,-E,), E=3KkgT/2, E -E A, 
(4.3.3a-c) 


WO) =(303/27º).K 3/4] X)l-expl3/(4|XP1, 34]|XP=27,, 
(4.3.3d,0) 


sl? Pp 
SC0=(27/7') [ar DG/8Xt-nIt, (4.3.3) 
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sendo K (z) a função modificada de Bessel [Watson (1966)] 
e D(z) a função de Dawson [Abramowitz e Segun (1969)]. 

Comparando as nossas predições quanto-mecânicas, dadas 
pelas expressões (4.3.1,2), com as correspondentes fórmulas 
clássicas semi-empíricas obtidas por Griem (1968), dadas, 
respectivamente, pelas expressões (24) e (35) (de seu artigo), 
vemos que em vez das funções semi-empíricas (“Gaunt factors”) 
B(X) e BalX), consideradas por Griem, elas são 
substituídas, respectivamente, por W(X)e S(X). 

É importante observar que nosso tratamento somente 
pode ser aplicado quando a segunda ordem de Born é satisfeita 
e quando as distorções das ondas planas |d) produzidas pelo 
campo do íon tornam-se desprezíveis. Isto implica que as 
expressões (4.3.1,2) são válidas somente para X >> 1. Na 
Figura 4.3.1, são mostrados as funções W(X)e S(X) para 
X| >5e os “Gaunt factors” g (X) e g(X). Como podemos 
ver nessa Figura 4.3.1, para X| >> 1 há um bom acordo entre 


os tratamentos semi-empírico e quanto-mecânico, mostrando 
que os cálculos semi-empíricos nessas condições são fidedignos. 


2.0 


WX (—) ToolX) (000) 


10 Ss (X) t---) Top 1X) (000) 
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Figura 4.3.1: As funções quanto-mecânicas W(X) e S(X), para |X| e 


e os “Gaunt factors” semi-empíricos 8, (X) e G,(X). 


Ainda dessa Figura 4.3.1, verificamos que para |X | J1 


os resultados clássicos semi-empíricos podem ser obtidos pela 
teoria do “cutoff” desenvolvida por Griem (1974), na 
aproximação dipolar, considerando o parâmetro de impacto 
crítico como p |. =h/mv. 

As fórmulas de Griem (1968, 1974) e versões 
ligeiramente modificadas delas têm sido aplicadas para estimar 
as formas de linhas de átomos neutros e íons em plasmas 
[Dimitrijevic e Konjevic (1981), Puric e Lakicevic (1981), Hey 
e Berger (1981), Dimitrijevic e Sahal-Brechot (1983), 
Dimitrijevic (1985), Konjevic e Pittman (1987), Djenize, 
Sreckovic, Platisa, Labat, Konjevic e Puric (1990)]. Dos vários 
tratamentos que têm sido desenvolvidos para calcular os 
deslocamentos Stark (“Stark shifts”), o tratamento semi- 
empírico de Griem (1968, 1974) para linhas de íons fornece a 
mais transparente demonstração da conexão (simplificada) entre 
as secções de choque de colisão eletrônica e os parâmetros do 
alargamento de linhas [Wiese e Konjevic (1992)]. A 
aproximação semi-empírica não fomece, algumas vezes, 
resultados muito precisos, mas ela contém os elementos 
essenciais das condições físicas do sistema átomo-plasma. Ela, a 
aproximação semi-empírica, apresenta fórmulas simples que 
podem ser facilmente calculadas, dando resultados razoáveis 
[Wiese e Konjevic (1992)]. 
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Apêndice C4-I 
CÁLCULO DO ELEMENTO MATRIZ <qlv'lã > 


Tomemos a expressão (4.2.2.2): 


<qivlãi> = 5 [d'R exp(i Ag. R/h) x 


x [ze exp(-R/(p)/R-D e? exp(-IR — E 1/05) /R — alas , 
em que: 


A = fe. fd R exp(i Aq. R/hi) exp(-R/(p)/R, 


B= E. DS fáR exp(isg. R/hi) exp(AR — El/0p)/IR -tl= 


= =. > exp(iag.7,/h).Jd' (R — 7) x 


x exp(ig. (R — %) /h). exp(AR = 1/05)/R — 1. 


Usando coordenadas esféricas [Gradshteyn e Ryzhik 
(1965)], virá: 
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Ze? f rr 
= Le | | IR dR exp(iAq R cos0/h) exp(-R/(5) x 


Z 


x sen9 dO dy = € E [ec + 2] 


O cálculo de B será análogo ao de A, bastando para isso 


fazer: 
R-t=i, 
então: 
2 
pet Do Ag) 

Ma) rag] ; 

Assim 
2 
<gIV'Iq> = € 1 a: > exp(iAg . 0) 
+ Aq | a 
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Apêndice C4-II 


CÁLCULO DE Fe (Ag) e Fo(Ag) 
Vamos partir das expressões (4.2.2.5) e (4.2.2.13): 


En (Ag) = At bo Cum x 
My 
xIZ dj —<4]My; > exp(iAg. z, /h) 10,3,M, >P 


Re(Ag) = Ag >, Cum * 
My 


x|Z — <aJMy;| 5 exp(iAq z,/h) o49M>| % 


x|Z — <a JM, > exp(iAq z, /h) GM : 


com: 


ASAE, ++) ]O, 


Cum = CF, MFOM$) . C(H:M,OM)) . 


No caso particular em que My = Mp = 0, teremos: 
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Coo = C(JrlJy; OOO) x C(Jj1JF; OOO). 


Usando as propriedades dos Coeficientes de Clebsch- 
Gordan [Rose (1967)], virá: 


Coo == "E . [(23,+1) (23,+1)] Px (153,00) x 
x 0(1,1/:1000) > Cop = 1/AIF. 


Portanto: 


Fe(Ag) =1Z dk, — <0G]KOID exp(iAg. z,/h) Lo,],0>P., 


e. 


FE(Ag) = E — <a,J,01 > exp(iAqg. z,/1) | 10>) x 


x E — <oJpOI > exp(iAg . z,/h)| 01,05 
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Apêndice C4-II 


OBTENÇÃO DAS FUNÇÕES Wk(A:Bkn) € S xn(A:Bkn) 


Vamos partir da expressão (4.2.2.7): 


: B lá et 3=493= q) 
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= ff n Eca H, (Ag) 
J=||dqd q-ex E RR RS 
) qa q 2m [to Es Aq” | 


xP|5"SE |. 
=" EMA, 


a) Função Wkn(A:Bkn) 
Calculemos inicialmente a integral I: 


Fo (Ag) > AS +A. | 
[to + ag?) | 


I=n]] d'q d'q- exp [ a . 


Façamos a seguinte mudança de variável [Cattani 
(1968)]: 


então: 


dqd'g= Id QUA = 5 E QU'Ag. 
Portanto: 


=4 [la(cos 045) d(045) fafeos 85) dfcos 05) x 
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x If a(ag) aa GF ag ex [E (O + 9) | 


x ng +] x 
2m : 


em que 0. Ong» % e do são os ângulos polares dos 


Ag 
2m +] Ê 


vetores Aq e Q, respectivamente. Chamando de O o ângulo 
entre Q e Aq e tomando o eixo dos z na direção do vetor 


momento transferido Aq , teremos: 


Et . F(Ag) o? poe EM nBefas ANE. 
et q ore] 


x | d(cos 9) exp - Ê (QAg cos o] ô Em cosB + Aa) E 
Fn(Ag) 
[Gu 5)! + aqº| 


exp |- É (ag) |x | 400. op(- BO), 
am (sab [x | im) 


Qmin 


Pa CR q 
=21m en É e | | Aq d(Ag) x 
0 


em que Qmin é calculado da expressão (cos O)min = [ ae | 


Integrando em Q , teremos: 
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'm? BA, r E Ag 
t-Smpm”. ex( Epi) ag a(s) — Fita 
$ (ht) + Aq” | 
2 
Bo BA 
exp|-s— (Ag) — 
-& 2(Aq)' 
Fazendo: 
1/2 
8 BA, , 
[5] o a 
[nto rag |= Et) 2 =(h/t,)(BI8m), 
teremos: 
1/2 
r dx 8m Y 
I=mm: 27.) E F aa — fin 2 | 
m-exp(27.,) [rey (5) : eu É e] 
Portanto: 


DE ou cd 
2miN(5h) El = 2mi Ne* x “E D Win (B; 415). 


21m 21 | 2m 


em que: 
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oo 1/2 
Wo Bin) = exp(2n,) | LES, (62) À E 


0 (a?+ x) 
x exp (10/82) 


b) Função Skn(A;Bkn) 


Calculemos a integral J: 


SEA E Ba H, (Ag) 2m | 
J= ciegen(-5 fa pub of in 
JJ [057 + Aq] q'—-q +2mA,, 


Seguindo as mesmas transformações do caso anterior, 
teremos: 


ie Ha od 


x] a (cos 6) exp [- É (Qua 058) Po dama | 


A integral em cos O = z torna-se: 


+ 
2 
P | dz exp - + (0 sq) QAgz FORTE E 
—1 1n 
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2mA. +QA 
Pa fd a 
OAg l 2mA q OAg 


(1 og) (| OM) 
+ ETR + 2ma, | | ga) 


in 


em que foi usada a representação em série da função: 


exp [da (Q447)] =1+ 5 a [ É. (Qua9 


e a seguinte mudança de variável: 


Portanto: 
J=2mmr exp (2 yin) (+45) 
com: 


F, (Ag) 


, = JJ Q dQ Ag d(Ag) x ; 
[neo + sa? ] 


x 


cer dio tuo) 
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1,=[JQ dQ Aq d(Ag) + E | exp im [ag á 


o (1 l as E ma SÊ 
A a (Cria) ( e Ss] aU E ea | 


Para realizarmos a integral J1, integraremos por parte, 
dividindo o intervalo de integração O > co em duas partes: O 
> Qo, em que Aq.Q —-2.m. Ajn<0, e QG 500, com 
Ag.Q — 2.m .A;n 2 O. Então: 


- — 4m eo Bal.  Fnlda) 
É B Fada P| 8m (aa [1/05 + aq”) ê 


Og 
Ea E QAg+2mA, 


oo 


16m'A, FA 
= ES (aa d(AQ) - expl- Ra (ag) | — HolMg) o x 
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Integrando em Q [Cattani (1968)], virá: 


Tre Dr [Ag d(Ag) - exp E É (aq) -«exp(—b)- 


GRC 


Taro sap om [tio 


em que: 


Agora, integrando em y [Gradshteyn e Ryzhik (1965)], 
teremos: 


Fm (Ag) 


- 8m f B (aq) 
dai - im (MM | | [to] + sg'] 


«exp (— b) x 
0 


Usando a definição e propriedade da função de Kummer 
[Slater (1960)], resultará: 
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E E 
à lados À do o (to) + aq] 


Tb - D(1,3,2:-b) . 


Chamando: 


alii 


e usando a relação entre a função b e a função Erfi (uy) 
[Tricomi (1954)], virá: 


= [des 1 e E = o A (83) +] eme] | 


Calcularemos, agora, a integral J> 
binomial, obteremos: 


2 (Dos If, OM) QAg ) | 
Dt ha |! timaç) (o Imã) | 


2q+n+ 2q+n+H 2q+H 
-25 E gd Bra) [ Qua [E R(f). 
nq=0 À n! (2g +n+ D)! 2 2mA, 


Derivando a expressão acima em relação a PB, teremos: 


. Usando a série 
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25 (enfrro ne BA E B QAq aaa 

E q=0 n! (2q+1)! f Bl âm 
E “Ban B QAq 

E exp > senh ea : 

Integrando a expressão acima com a condição de R(0) = 0, 
virá: 

a vao: An? QAgz 

R=-2 | E exp [- gº sem (Sia : 


Integrando em Q [Roberts e Kaufman (1966)], 
obteremos 


[o-do-en[-2.0º | sem Ca) 
8m 4 


m 


2 
— V21m (Agz) 
= p32 Ag z exp Bêm | 
Portanto: 
amo B A. 
= 2 V21m in 
Jr= = Bi? pa exp | >. | x 


a — ho) Bial 
[Agdiho) lee a »| gm (Aa) 3] 
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12 
Colocando: Ag = (6) x e z=PBt, teremos: 


oo 1/2 
=| ads o (50) 8 (9) a 


1 
x] dt exp (- A E a RL E 
0 
Agora, realizemos a integral em t: 


Í dt exp (- A ris 
0 


Usando a definição da função Erfily) [Jahnke, Emde e 
Lôsch (1960), Tricomi (1954)], qual seja: 


y 
Erfi(y) = [ exp(t”) dt, 


a igualdade anterior ficará: 
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Portanto: 


J --JT[ 40 E (82) ".» “exp a bg 
E o (24x) PIB x? 


x (em - E to) -— Erfi [te 
XxX XxX 


De posse dos valores J1 e J2, o valor de J será: 


J=2m7 exp (2 Y%;n) =] FE - exp - a — o pe 
0 


Chamando: 


y=x—-YWh/X, 


e usando as definições da função Erfi(y) e da Integral de 
Dawson [Abramowitz e Segun (1969)]: 


NT 2 
D(y) = exp(-y ) Í exp(t ) dt = D(x:Y). 
0 
J tomará a seguinte forma: 


J = -2mn?? Í sa E, (62) D(x:Y). 


0 (a? + xº 


Apêndice C4-IV 


CÁLCULO DA LARGURA E DO DESLOCAMENTO 
PARA A LINHA 4ºs — 2*p DO HÉLIO IONIZADO 
(He) EM UM PLASMA 


Usando o formalismo quântico que desenvolvemos no 
Capítulo 4, obtivemos as larguras e deslocamentos de várias 
linhas do He* (em Á: 3188, 3389, 3964, 4121, 4437, 4713, 
5016, 5048, 5876), por intermédio de cálculos muito extensos e 
laboriosos. Para ressaltar essa afirmativa, vamos apresentar 
neste Apêndice somente os principais resultados do cálculo 
da largura (Av) e do deslocamento (S) de apenas uma linha, a 


4ºs — 2ºp (4713 A) do hélio (He*) em um plasma. Os detalhes 
desses cálculos podem ser vistos em Bassalo (1975) e Bassalo, 
Yamamoto e Cattani (1975). 


1 Termo de Primeira Ordem 


O termo de primeira ordem do deslocamento S/ 


representado pela expressão (4.2.1.2) é zero porque Vi = Vr = 
O, como pode ser verificado facilmente levando em conta a 
simetria da distribuição de carga dos elétrons no átomo no He*. 
Este resultado vale para todas as linhas do He* analisadas 
por nós [Bassalo (1975))]. 


2 Termo de Segunda Ordem 


Em nosso cálculo, o maior problema foi obter os fatores 
de forma Fkn(Ag) e Fir(Ag) definidos pelas expressões 
(4.2.2.5,6) e (4.2.2.14), respectivamente, porque o número de 
funções de onda do hélio que são encontradas publicadas é 
irrisório. Por isso, numa primeira aproximação, aqueles fatores 
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foram calculados usando as funções de onda hidrogenóides 
encontradas em Pauling e Wilson (1935). Destarte, verificamos 
que somente as transições virtuais 4s > 4s,4s > 4p, 45 > 4de 
4s — 4f contribuem para a largura e o deslocamento eletrônicos. 
As demais transições são muito menores e podem ser 
desprezadas. 

Numa segunda aproximação, os fatores de forma foram 
recalculados usando as funções de onda (1s)(4s), (1s)(4p), 
(1s)(4d) e (1s)(4f) obtidas numericamente com auxílio do 
método de Hartree-Fock-Slater com “exchange” de Kohn-Sham 
[Leite e Ferreira (1971)]. Verificamos que as diferenças entre os 
resultados precedentes e os novos são desprezíveis. Ou seja, em 
nosso caso o cálculo dos fatores de forma [Fass (Ag) (! = s, p, d 
e f)] com auxílio das funções de onda hidrogenóides é plenamente 
satisfatório (ver seus valores mais adiante). A carga nuclear 
efetiva Zerr para um elétron com o número quântico principal n 
= 4 foi tomada igual a um [Bethe e Salpeter (1957)]. 

Deste modo Ave e Se, dados pelas expressões 
(4.2.2.15,16), ficam, respectivamente, definindo a variável 


y = (B/8m)2.2a,/h: 


em que 9=2a,//p., ao é O raio de Bohr, E=(2mB)“a,/h e 
n=(B/2m)” “h/4a,. As diferenças de energia Ass4y (medidas 


em cm» entre os estados 4s e 4/(/ =Pp,de f)são iguais 
as 092108 cem, STA dO” = e SiS 4 10º cem; 
respectivamente [Bethe e Salpeter (1957)]. 


: (0) (0) O (0) 
3 Fatores de Forma: EF, E, Sha OR 


No cálculo dos fatores de forma, as integrais serão 
realizadas usando-se as tabelas de Gradshteyn e Ryzhik (1965), 
com a seguinte transformação de coordenadas: 


— 2aç 
—h 


Aq . 


3.1 Cálculo de F2 


48>4s 


Calculemos inicialmente 


RR = [ exp(iAgcose/h) Vaso Vaso dº F = 
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E E da — 36r 
= = Es 7 fr dr exp[ a E a a, 


| | exp(iAgcos8/h) d (cos0) = 


1 hq r Aq 
= : r dr ef] ; sen( 
1152(2a,)” Aq Zã 


=h+b+ cases + Lc, 


com ea : 
1152 (2a,) Aq 


I =A ] rdr ef] , 2424sen| Sar > 


o 


n=(y+ DS; 
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L=A rdr expf- 5 | | 2a(- 3 Jen( Sr) > 
L = (3/)(y7 - 3(yº + D; 
L=A rdr expf 5 =). 24, Ea sen( ddr) = 
L= 67 My” + DR 


E A 
-—+ Aa 
1 = A rdr expf 2a =). 2. | 5 sen( f ) = 


L=(1/9(y* — 109º +5)(yº + DS; 


I,=A rdr expf- =] (- ae), 24. sen( Sr] — 
= (3/7 - 3º +): 


lE=A ] rdr expf- =| [- de [- de sen( dr) EN 


= 277º = 1)(y” + : 
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E 2 
L=A | rdr expf =] ; [- e l Ea ; sen( Sr) > 
-5 
L = —9(y* = 1097 +5)(y? o 1) 5 
SRT EE oo O PR PR a 
I; = À ] rdr ef =] Se) | a sen( f ) > 


I = (15/2)(3yº — 109º +3)(y? + nº; 


(o) 


na RE op pr cu DR a 
L = A | rdr expf 2a J | De 24 sen( f ) E 
—4 
L = 6 -1y +1): 
ja 2 
Lo = A | rdr expf- E ] : E : [- dé) é sen( Sr] > 


ho = 9 (y! — 109245) (2 + 1); 


f r 121? 121? Ag 
ea food) [85] [85] 


I = 30 (3yº — 109743) (yº + nº; 


2 3 A 
Il, = A rdr expf- a =): | : - EE 5| l sen( ddr) => 


Lo = (15/2) (9º — 21yº+35yº-7) (yº + NE 


3 
1, =A rdr expf 2a =| . - CE 5| 24: sen( Sr] =: 


—5 
Ls = (1/2) (yº = 107 +5) (7º +1) : 
' ; 7 
teca free) [ag] (38) cof) = 


La = (15/2) (3º — 109º +3) (yº + nº; 


E 3 2 A 
no a ex *) - Za 5| | - E | sen( Sr) = 


hs = (15/2) (9º — 21y* +35? = 7) (yº + NE 


E 3 3 A 
I6=A | rdr expf- =] - E 5| - TE s| sen( Sr] E 


ho =35 (9º +19! 77 +41) (yº + . 
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Sendo: 
Eca Fa |! o Laisa , 
teremos: 
. (4y*-43yº +31) (19yº-306y! +355y? -40) 
Ads | A, AÃO POD HP AG DST 


ao qa) 


a 
(15yº -335yº +455y* +175y 70) 


(7+1) 


ou: 


po (y) =| 1-47"-46y"+166y'—215yº+109yº —19yº + 
4s—4s () — O 


8 
(y7+1) 
3.2 Cálculo de Eai 
Sendo: 
O o E 
Ras E | fus ap q 


com: 
Esso | exp(ixq r cos0/h) Waso Wap dr, 


e usando a mesma técnica do item 3.1, teremos: 
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-16 2 
Ep sap (0=45(9"+1) (2y"-18yº +43] -37yi et dy" y) 


3.3 Cálculo de F,. ,4a 


Sendo: 


2 


Es sad = [sai > 
com: 
Lis4d E [exp(isg r cos0/h) Vaso Wado ds 


e usando a mesma técnica do item 3.1, teremos: 


E sua)= (199 +1) (647! +s68yº 99299 + 5367! 807)" 


3.4 Cálculo de EF. ,4r 


Sendo: 


2 


2 


Es s4t = [Fes4e 
com: 
Lisa = [ exp(isq r cos0/h) Vaso Waro PT, 
e usando a mesma técnica do item 3.1, teremos: 


ES, (0)=(1/320)(y2+1) “(-1603º+1470y7 —2240yº +560y').. 
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4 Largura e Deslocamento Totais 


Com as expressões para AV, e S, vistas acima, 
podemos calcular a semilargura e o deslocamento produzidos 
por colisões eletrônicas. Para calcular o efeito dos íons usaremos 
a aproximação de GBKO (1962) e Griem (1962b). Segundo 
esses autores, a semilargura total Av e o deslocamento total 
S são dados pelas expressões: 


Av=Ave +1.75a(1-0.75R) Ave, 


S=Se+2.00(1-0.75 R) Ave, 


conforme vimos no Capítulo 3 [expressões (3.4.5,6)]. Os 
parâmetros « e R estão definidos e tabelados nos artigos de 
GBKO (1962) e Griem (1962b). 

Usando as expressões acima, verificamos que, nestas 
condições experimentais, em que somente íons de He* estão 
presentes, 85% do alargamento e do deslocamento são 
provocados por colisões eletrônicas. As transições virtuais 4s 
S 4s,4s 5 4p, 45 > 4d e 45 > 4fcontribuem para Av ou 
S com cerca de 8%, 90%, 1.7% e 0.3%, respectivamente. 
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CAPÍTULO 5 


FORMALISMO SEMICLÁSSICO 
CONVERGENTE 
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5.1 Introdução 


O alargamento e o deslocamento de linhas espectrais 
emitidas (ou absorvidas) por átomos neutros e íons em plasmas 
são causados por numerosos processos, conforme vimos no 
Capítulo 1. 

Nossos estudos limitam-se a plasmas localmente 
homogêneos e em equilíbrio termodinâmico. Neste Capítulo 
desenvolveremos uma teoria que denominamos de Formalismo 
Semiclássico Convergente que pode ser aplicada para átomos 
neutros e íons. Assumiremos que as condições de temperatura e 
pressão do plasma são tais que alargamentos e deslocamentos 
das linhas são provocados principalmente por colisões 
(eletrônicas e iônicas) e pelo efeito Doppler, sendo esse último 
processo cerca de uma dezena de vezes menos importante que o 
primeiro. 

A forma da linha devido ao efeito Doppler é descrita por 
uma gaussiana [Vainshtein e Sobel“man (1959)] dada por: 


Fp = (12 qp)" exp [- (0- 05) 2/05], 
(5.1.1) 


em que a largura Doppler da linha é op = wo v/c, sendo q% e 
v, respectivamente, a frequência da linha e a velocidade média 
das partículas emissoras. 

O efeito provocado pelas colisões pode ser calculado 
quântica ou semiclassicamente, resultando em uma forma de 
linha F(w) que contém todas as contribuições dos diferentes 
perturbadores constituintes do plasma. 

A forma de linha observada I(w) é, então, calculada 
como a convolução das curvas obtidas para cada um dos efeitos 
individualmente, isto é: 
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K0)=[" (0) FO-0) do. (5.1.2) 


Muitos resultados experimentais foram obtidos em 
diversos tipos de plasmas e simultaneamente várias teorias 
foram desenvolvidas para explicar a contribuição das colisões na 
forma das linhas observadas. 

Para densidades eletrônicas relativamente baixas as 
linhas observadas resultam de transições que obedecem às regras 
usuais que se verificam em átomos não submetidos a pertur- 
bações externas. Nessas condições as linhas observadas são 
isoladas e denominadas permitidas. Quando a densidade de íons 
do plasma se torna relativamente alta, a perturbação iônica pode 
induzir transições que, nas condições anteriores, não eram 
permitidas. Estas novas transições são denominadas proibidas 
[Vainsthein e Sobelíman (1959), Griem (1974), Barnard, 
Cooper e Smith (1974)]. 

Assim, ao lado de uma linha permitida pode surgir uma 
proibida, ocorrendo a superposição das duas linhas. Contudo, o 
estudo da forma de linha neste último caso requer uma 
abordagem diferente daquela do primeiro caso. 

Griem, Baranger, Kolb e Oertel (GBKO) (1959), Griem 
(1974) e Barnard, Cooper e Smith (1974), por exemplo, desen- 
volveram formalismos para tratar dos dois casos. 

Neste Capítulo 5 utilizaremos uma Teoria Semiclássica 
Convergente apenas para estudar linhas isoladas [Cattani 
(1972), Walder (1982)]. 

Nos próximos itens mostramos o formalismo geral 
para obtenção da forma de linhas espectrais em plasmas e as 
hipóteses feitas nesse formalismo, que dão origem à Teoria 
Semiclássica Convergente, objeto principal de estudo deste 
Capítulo. 


129 


5.2 Teoria Geral de Forma de Linha Espectral 
Isolada 


De acordo com que estudamos no Capítulo 3 (item 3.1), 
a forma de linha espectral isolada F(qw), originada pelas 
colisões dos perturbadores com o emissor no plasma, é 
calculada por meio da transformada de Fourier da função de 
correlação temporal d(t), ou seja: 


A 1fS io 
F(w)=Re by e o] (524) 


Designando os estados inicial e final por Wj e Wr, 
respectivamente, da linha isolada I > F, a função de correlação 
d(t) é dada por: 


H0=p,|<WOIdIy(O> - <y(DIdly >] , 


Av 
(5:22) 
em que pr é a matriz densidade para o estado inicial, do 
operador momento de dipolo elétrico que define a interação da 
onda electromagnética emitida com o emissor, t o intervalo de 
tempo da correlação e [ Jay indica que deve ser feita a média 
sobre os estados internos dos perturbadores, média sobre os 
estados translacionais dos emissores e perturbadores e média 
sobre as orientações relativas entre os sistemas emissores e 
perturbadores. 
A evolução temporal do sistema é descrito pelo operador 
unitário T(t), que transforma a função de onda de acordo com a 
relação: 
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T(O Iy(0) >= Iy(t) >. (5.2.3) 


Se introduzirmos o operador evolução na função de 
correlação temporal, podemos escrevê-la como: 


OD=pil <Weldlyç >| <ye ITC ly > + <p TON] 
(5.2.4) 


Consideremos o operador Hamiltoniano do sistema 
átomos neutros (emissores), íons e elétrons dado por: 


H=H+ VD + Vi(d, (5.2.5) 


em que Ho é o Hamiltoniano do sistema na ausência de 
perturbadores, e Ve(t) e Vi(t) são, respectivamente, as interações 
dependentes do tempo do emissor com todos os elétrons e íons 
do plasma. 

A interação entre os átomos neutros não é levada em 
conta porque é desprezível comparada com Ve(t) e Vi(t). Por sua 
vez, o operador evolução temporal T(t) obedece à equação 
diferencial: 


na TO =[H, + Vo + VOTO. (5.2.6) 


Agora, vamos definir o operador T(t), que leva em conta 
as perturbações eletrônica e iônica, como um produto de 
operadores Uc (t) e Ui(t), ou seja, 


T(O) = Uc(t) Ui(b). (5.2.7) 
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Substituindo a expressão (5.2.7) na expressão (5.2.6), 
obteremos as seguintes equações diferenciais para Uc(t) e Ui(t): 


ih q VO=[H+VOVO, (528) 


nau, D=UMDVIDU(DUD. (5.2.8b) 


Como os elétrons se movem com velocidade muito 
maiores do que a dos íons, estes podem ser considerados 
parados durante cada colisão elétron-emissor. Devido a isso 
consideramos em nossos cálculos os íons estáticos [GBKO 
(1959), Griem (1974)], isto é, consideraremos Vi(t) = Vi. 
Portanto, o operador Ui(t) fica escrito de forma mais simples: 


U()=exp 


(8, ev | (5.2.9) 


O efeito provocado pelos íons na dinâmica dos elétrons 
está contido no pequeno deslocamento AE do nível que eles 
produzem. Porém, como Te AE << 1 (em que Tt; é o tempo de 
colisão elétron-emissor) e AE é pequeno comparado com o 
espaçamento entre os níveis, assumiremos que Ho +Vi; = Ho. 

Assim, a expressão (5.2.8b) se simplifica: 


AZUL = VD U.(D, (5.2.10a) 


em que: 
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V; (O=e/ 2H] V.(t) exp [tia . (5.2.10b) 


Esta aproximação permite que os efeitos de colisões 
eletrônicas e iônicas sejam desacopladas e estudadas 
separadamente, como passamos a fazê-lo. 

Desse modo, a função de correlação assume a forma: 


di (D=p ey ldhy, >P| ey IUZ (OU ; Nye><W UU CO | 
Av 
(5.2.11) 


Introduzindo os projetores: 


y m > W m | e >hy, ><W | Há 


a expressão (5.2.11) ficará: 


O (0=paeyeldyç > E | WE UZO ><a UM Or 


x <y U(Oy,> - <y IU(Oly, dl (5.2.12) 


Av 


em que a média indicada (Av) pode ser feita separadamente para 
os elétrons e íons, como indicamos a seguir: 
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elétrons 


dir (D = pil<yçldhy >P [WU Oyo> <p, U. (ON | » 


Av 


íons 


x ya UM Olyr> “<td TU C(Oy, >| 


Av 


(5.2.13) 


Para calcular as médias indicadas acima precisamos 
definir como são as perturbações geradas pelos elétrons e íons. 

Como dissemos, os íons são considerados estáticos 
durante as colisões eletrônicas. Assim, podemos afirmar [GBKO 
(1959)] que as colisões eletrônicas ocorrem na presença de um 
campo elétrico estático F criado instantaneamente pelos íons, 
dado pela distribuição de Holtsmark W(F), sendo esta 
aproximação chamada quase estática. 

As funções de onda W, do sistema são autofunções do 
operador Ho, com autovalores Er, isto é: 


Holy, > = Ecly, >. (5.2.14) 


Admitindo que o campo elétrico estático F pro- 
voca no emissor um efeito Stark de segunda ordem, 
temos 
LES anP 


U(Dly,> = e” 


Ny > + Lhy><yiVilyo JE, Edo 


jm 


(5.2.15) 
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em que: v=d.F e F=5E(b. Portanto, teremos: 
j 


Epa? Xt 


<yIU;(Dly,> = e dh (5.2.16) 


Desse modo, o termo correspondente à perturbação 
iônica fica: 


[Vl UMONr>  <yiU, (Oy, aa e 


íons 


ne Aar-apF2 
=, e a sa (5217 


Av 


em que W=(Ep-E?)/h e, consegientemente, a função de 


correlação se reduz a: 


ia 


dl) = prxyeldiy,>P [VEIUÍONy:> “<p UU Oy, 


> 
Av 


| K > E íons 
eee” o (52.18) 


Av 


A contribuição das colisões eletrônicas na função de 
correlação pode ser calculada a partir da expressão (5.2.10a) 
para o operador Uc(t). Dessa equação temos que Uc(t) é dado 
por: 
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U.(D=3 exp [-À [vio ti] (5.2.19) 


sendo 3 o operador ordenação temporal ou operador de 
Dyson. A perturbação criada pelas colisões entre os elétrons e o 
emissor está contida no operador V;(t), e as colisões devem ser 
consideradas na ordem cronológica em que ocorrem. 

O problema da ordenação temporal das colisões é 
extremamente complicado, pois, devido ao longo alcance das 
forças eletromagnéticas, o emissor está interagindo, num mesmo 
instante t, com muitos perturbadores e, além disso, a 
perturbação criada por uma colisão num instante t + At encontra 
o emissor perturbado pela colisão ocorrida no instante t, o que 
gera um efeito cumulativo das colisões. Além disso, devemos 
considerar as interações entre os perturbadores. 

Para resolver esse problema, Baranger (1958a,b) 
introduz duas aproximações: a primeira considera os 
perturbadores estatisticamente independentes, isto é, não existe 
correlação entre eles; e a segunda supõe válida a aproximação 
de impacto. 

A aproximação de impacto permite resolver o problema 
da ordenação temporal das colisões e constitui a hipótese 
fundamental da Teoria Geral de Forma de Linha Espectral. 


Aproximação de impacto e Condições de Validade 
Fundamentalmente, a aproximação de impacto consiste 


em supor que duas ou mais colisões fortes não ocorrem 
simultaneamente e que o intervalo de tempo entre elas é muito 
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maior que o tempo de colisão 7... Esta hipótese leva à conclusão 
de que a interação média entre o emissor e perturbadores é fraca. 

A aproximação de impacto é válida quando estiverem 
satisfeitas as condições acima, o que quer dizer que o 
alargamento We e o deslocamento Se provocados pelas colisões 
com um dado perturbador, geralmente elétrons, devem ser 
menores que o inverso do tempo de colisão. Como we e Se são 
da mesma ordem de grandeza, podemos escrever que: 


Q<T te E eo'=%m cc 1, (5.2.20a-c) 
v 
em que pm é a distância média entre os emissores e v a 
velocidade do perturbador. 

Baranger (1958a,b) mostra que, como consegiiência des- 
sas aproximações, a interação dependente do tempo entre os 
elétrons e o emissor pode ser substituída por uma interação 
efetiva constante representada pelo operador não Hermitiano 
Hrr, calculado a partir dos elementos de matriz Sn e Srr da 
matriz de espalhamento S dos estados inicial e final. Ou seja: 


[<p IUÍONy> : <yvU.(Oy, >|, = el!t)Hrrt = e Me +iSet 
(5.2.21) 
em que: 
Hr= [1-SnSfe |, »  Sy= <ViSy,>, 


(k = I ou F) são os elementos diagonais da matriz S de 
espalhamento entre o elétron e o emissor. A semilargura Qe e O 
deslocamento Se, devido às colisões eletrônicas, são dados por: 


O. = — Im Hr 5 Se =R. Hr (5.2.22,23) 
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Escrevendo explicitamente a média [1-Sy Si], para 
V 


um perturbador qualquer, temos: 
Hr=—i27N [Pdbb [“dvv Fw) Sp, - <nl(l-SpShe)ln>, 


(5.2.24) 


em que N é a densidade de perturbadores, b o parâmetro de 
impacto da colisão, v a velocidade relativa entre emissor e 
perturbador, F(v) é a distribuição de velocidades de Maxwell- 
Boltzmann, |n> os estados internos dos perturbadores e pn a 
probabilidade de encontrar um perturbador no estado |n>. 

Uma análise da expressão (5.2.24) mostra que existe 
uma correlação entre os estados final e inicial dos perturbadores, 
pois: 


<nISy Sppln> = X<nlSgln'> - <n1Sgeln>. 
a 


Quando não há mudança do estado do perturbador du- 
rante a colisão, o que corresponde a considerar apenas n = n',a 
matriz de espalhamento S pode ser escrita na forma 
exponencial, definindo o que denominamos de Formalismo 
Semiclássico Convergente, como mostraremos no item 
seguinte (5.3). 

No caso em que há mudança de estado dos 
perturbadores, todos os termos do tipo n nº devem ser levados 
em conta e, consegiientemente, a matriz de espalhamento S não 
pode ser calculada na forma exponencial, sendo necessário 
utilizar o tratamento de Anderson (1949). 
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Como os perturbadores no nosso caso são os elétrons e 
portanto não possuem estados internos, escrevemos Hr como: 


Hr=—127N [Pdbb [Pdvv F(y) - [1-SySie). (5.2.25) 


Agora, podemos voltar à expressão (5.2.18) para a 
função de correlação Orr(t) e calcular o termo correspondente às 
colisões eletrônicas. Então, substituindo a expressão (5.2.21) na 
expressão (5.2.18) e fazendo a média sobre a contribuição 
iônica, resultará [Walder (1982)]: 


: : 21 Ds: 
= + uia ssestrcar : 
IF(O = ye ldiy, >. [“w(F) dF-[“dt' x e ê 

0 0 
(5.2.26) 


Finalmente, podemos calcular a forma da linha, que leva 
em conta os efeitos dos íons e elétrons, resultando em: 


1 Aura oe+Se+2E ca) 
F(w)=-R, [dt e”! ["w(F) dFe : 
To 0 


(220) 


em que foi usada a notação de GBKO (1959), quando 
escrevemos O -—Op = ER sendo C4 o parâmetro Stark 
e 


quadrático. 
Portanto, integrando no tempo a expressão acima, virá: 
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w(F) dF 


F(w=Se[- (5.2.28) 


ú cÊ+|w-w-S 


2 
21 

e 2 C, E] 
e 

Novamente introduzindo a notação adotada por GBKO 


(1959), obtemos a forma reduzida da linha j(x,0) em unidades 
adimensionais: 


joso)= " (5.2.29) 
0 1i(x-0p) 
em que: 
4 
= WD WeDS. ia 21C , + 3 É 
= CD, É = OD, po í Pm= 471N ? 
(5.2.30a-c) 
e: 

wH(B) = wH(F/F), F=261eN?. (5.2.30d,e) 


O alargamento qe e o deslocamento Se eletrônicos dados 
pelas expressões (5.2.22,23) dependem do operador Hyr que, 
por sua vez, é calculado a partir dos elementos da matriz Sr e 
Syr » estudada no próximo item. 

Outras equações para o alargamento e deslocamento 


devido às colisões iônicas são discutidas e mostradas no item 
5.6. 
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5.3 Cálculo da Matriz de Espalhamento S na Forma 
Convergente 


Neste item mostramos como calcular o elemento da 
matriz <nlSln> da matriz de espalhamento S e como esta pode 
ser escrita na forma convergente. 

Como o Hamiltoniano do emissor é dado por 
H=Ho+V(t), podemos obter a evolução temporal do sistema 
utilizando o operador unitário T(t) na representação de 
Heisenberg: 


Te =; (5.3.1a) 
em que: 
“Ho, “H 
lise a ve tee". (5.3.1b,0) 


Por sua vez, o operador de evolução temporal T(t) 
obedece à equação diferencial: 


.H 
=-"04 


Bo: 
nSTO e PN e PTD: (5.3.2) 


com a condição inicial T(-co) = I, sendo I a matriz unitária. 

Como sabemos [Davydov (1968)], os elementos de 
matriz Tmn(t) estão relacionados com a matriz de espalhamento 
S(t) por intermédio da expressão: 


t 
E f de Vim (t) 


TaD=S (De = (5.3.3) 
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Por outro lado, os elementos de matriz do operador 
Tron(t) satisfazem a equação diferencial indicada abaixo: 


0 q 
ns Tan(D = <mlih5-T(Oln > = 


-Ho “Ho , 
= Z<mle ? V(De 1? Ik><kIlT(Dln> 
k 


ou seja: 


SD Tool) = — > ValO Tot) eiomkt, (5.3.4a) 


em que: 
Vmk(t) =<mIV(Olk>, Tal) =<KIT(Dln>, 
Omk = (Em— Eo)/h.  (5.3.4b-d) 
A partir da expressão (5.3.3), podemos verificar que: 
al alo (5.3.5) 
Derivando em relação ao tempo t a expressão (5.3.3), 


virá: 


e 
E Ton(D) = SS) exp E! dt Vamít + 


E Vam(D San (D a d dt Van(O | (5.3.6) 
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Substituindo (5.3.6) em (5.3.4a), vemos que um 
elemento de matriz qualquer Smn tem a seguinte forma: 


DS an(O = VanlO Sal) E E VOS ka (O EXPO it Ma (0) 
k 
(5.3.7a) 


em que: 
Mn(D=5 [de am (t)— Va]. (5.3.7b) 


Como estamos interessados em obter somente termos 
diagonais de Smn, isto é, Sm e S*rr, vamos analisar o elemento 
de matriz Snn(t) que obedece à equação diferencial indicada a 
seguir: 


Is(ov-ivos( + DATONHO expfio, +]. 


ot” E h 
(5.3.8) 
Porém: 


E Va n(D expo, t+HNG(D] = Van(D San (D exp [io t+ (D] + 
k 


42 VaDSsa(D exp fio, t+ina(O]. 


em que o símbolo >indica a soma sobre todos os possíveis 
k 


valores de k exceto k=n. Assim, a expressão (5.3.8) fica dada 
por: 
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9 MEG 
om ==a z VÁ (DS a (D (5.3.9a) 
em que: 
Va(D = Va(D expfio,t + Ma(D]. (5.3.9b) 


Para resolver a expressão (5.3.9a) precisamos conhecer o 
elemento de matriz não diagonal Sen(t) da matriz de 
espalhamento que é obtido por meio da expressão (5.3.7a): 


a) pes do 1 * = 1 1 
BE SiaD = — dj Vin Sal — E Vil Smo(D (5.310) 


Integrando a expressão acima, virá [Yamamoto (1977)]: 


Sa (D=— 5 [ dit) Salt) — Ef at Sant) 


mzk —oo 


(3119 


Substituindo a expressão acima na expressão (5.3.9a), 
temos a seguinte equação diferencial a ser obedecida pelo 
elemento diagonal Snn(t): 


SSD = — DD [AVG Vinlt) San(t) — 
k 


—oo 


E mo Epa VÃO Vea Sao MS) 


k mk —o 
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Repetindo os cálculos dos elementos não diagonais de 
Snm(t) por meio de expressões similares à expressão (5.3.10) e 
substituindo na expressão acima, obteremos ESTO como uma 


série de potência de (ih)! que contém apenas o elemento 
diagonal Snn(t). Assim, [Murphy e Boggs (1966), Yamamoto 
(1977)]: 


b) LET USARÁ 
IrSmm(1) Ro E dE ViO Vinlt) Salt) + 


; t gi] 
+55 DU Tdi | dt VaO Vit) Vento) San(to) + 


k mk —o 


(5.3.13) 


Os termos de ordem superior na expressão acima podem 
ser obtidos tomando-se como modelo os dois termos dados 
nessa mesma expressão. 


Embora E Sm(t) seja expresso em termos de Sn(t) e não 


contenha quaisquer elementos não diagonais, a expressão não 
pode ser calculada facilmente, desde que os elementos diagonais 
aparecem dentro do sinal de integração. Entretanto, podemos 
usar a relação: 


Sm(1D=S(D + [ad S nt), (5.3.14) 


de forma que, substituindo-se a expressão (5.3.13) na expressão 
acima, resultará: 
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1 ft ' t ' ' * ' : 
Salt) = Sal + o Jdt [de Vo) Vin (to) Sia ty) Fe. 
k 


(5.3.15) 


Substituindo (5.3.15) e o termo equivalente Snn(t2) na 
expressão (5.3.12), teremos: 


o) = 1 ' ; ai E 
E na AV (O Vin(to) X 
x [Sn(D + SS! fi defó dt, VA) Vea (LD Sant) +... | + 
nn h2 : H “so 1 ink knt1 nn 1 


; t y 
E is s fat, ] dO Vit) Vin(to) x 
k es, 


mzk —o 


x EO + a [def de VA) Vit) Salto) + | 
k 


Aqui, novamente, colocamos Snn(t1) em termos de 
Snn(t) segundo a expressão (5.3.15), e assim sucessivamente, de 


forma que podemos colocar Snn(t) em evidência. Desse modo, 
teremos: 


ô Spa 
DES (D) = - opa ! ALVO Vin(t) + 


: t 1 o ; E; 
+ e E fd | dt VO Vim(t) Ven(to) +... | San (D > 


k mk so 


—oo 


(5.3.16) 
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em que explicitamos apenas os dois primeiros termos. 
Integrando a expressão (5.3.16) (lembrar que 


9S,(t) - dS,. (t) Ro 
dt dt 


Í seio ; dt - 22 fas VARO) Vin(t) + 
2. 


[a o db VÃO Vinlt) Vin(to) +... 


—oo 


Chamando: 
oo t 
At) = + no. á dt l dt VÍ Vito — 


EL Bs í dt í dt, 5 AVIÃO Vinlt) Vento) +... 


mzk —co — —oo 


teremos: 
San(+9)=S,(-00) nd, (5.3.17) 


em quen= Iou F. 

A matriz de espalhamento escrita na forma 
exponencial indicada pela expressão (5.3.17) é a essência do que 
chamamos de formalismo convergente. Como veremos a 
seguir, é possível descrever as interações fortes e fracas de modo 
contínuo em função do parâmetro de impacto b sem que 
apareçam divergências nas probabilidades de transição quando 
b > O (região de colisões “fortes”. 
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Na expressão de An(t) podemos considerar apenas o 
primeiro termo da expansão, pois a maior contribuição para a 
probabilidade de transição é dada quando V(t) é grande o 
suficiente para provocar transições entre os níveis de energia do 
emissor. 

Como a dependência da probabilidade de transição em 
V(t) é exponencial, quando V(t) é grande os termos de ordem 
superior são pequenos comparados com o primeiro termo e, 
portanto, podem ser desprezados. Desse modo, resulta [Murphy 
e Boggs (1966)]: 


+oo il 
Ab) A fd favíiO Ve(tp. 
E = —oo 


Levando em conta ainda que nnx(t) — nun(t) = O, a expressão 
acima ficará: 


oo t 
Ab o | de fat V oO Vento) ent, 
k -o —oo 


A integração em t; indicada na expressão acima pode 
ser feita considerando t — ti = T. Assim: 


oo t 
At E | de f dr Va Vin lt) ein, (5.3.18) 
—oo 0 


Como os operadores Vnk(t) e Ven(t-t) são Hermitianos e 
o potencial de interação pode ser escrito como a transformada 
de Fourier de uma função G(w) que descreve a interação no 
espaço de energia, virá: 
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1% : 
V(t-D=Vi(t-o), Go)=— | dt Ve”, 
at D=ViMt—) (09) Rd (1 
(5.3.19,20) 


e portanto An(t) fica: 
; gi Es ã 
Ao (28) D | dei [dt <nl | doG(o) esmik> x 
k 0 —oo —oo 


x«<n | | do (0) eiSCk>. 


—oo 


A integral acima tem como resultado [Murphy e Boggs 
(1966)]: 


A(O (22) 5 [ do [não -oHP(— )Jtenio (09) lk>? 
n da h E nk D,—O , 


(5.3.21) 
em que P indica a parte principal. 


Introduzindo as funções Ik e dk , a expressão acima 
será escrita na forma: 


A(D=5 Eidos (5.3.22a) 
em que: 


2 
T; 22) D<niG (o, >, (5.3.22b) 
k 
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+00 2 
“2, <n1G (O) kl 
O = P ) Vo (5.3.220) 


Em vista das expressões acima, a matriz de espalhamento 
Skkx(+oo) será escrita na forma: 


af es 
Sado E (5.3.23) 
em que tomamos, no instante inicial, Skx(+o) = 1. 


As expressões  (5.3.22b,c) que representam, 
respectivamente, Ik e O , mostram que essas funções 
dependem dos elementos de matriz das transições virtuais do 
nível considerado na presença da perturbação eletrônica. 


5.4 Cálculo do Alargamento e do Deslocamento 
Eletrônico 


Agora estamos em condições de calcular o 
alargamento e o deslocamento eletrônicos de uma linha 
isolada I —s F dados, respectivamente, pelas expressões 
(5.2.22,23), que dependem do produto dos elementos Sr e 
Sgr da matriz de espalhamento S . 

Esses elementos de matriz podem ser escritos usando os 
resultados do item anterior (5.3), na seguinte forma: 


ro; 
ad 
2 


<yIShyy> = Sy=e (5.4.1) 
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r + 5pidp 
<yp IS Iyp> = Syp=e * é (5.4.2) 
e, por conseguinte: 
+ = —STypidyp 
<vyiISiy ><yçIS Iyç = e E (5.4.3a) 
em que: 
Pr =I+ e, Or = O -—Or. (5.4.3b,c) 
Como estamos interessados no cálculo de alargamentos e 
deslocamentos de linhas emitidas por átomos neutros do plasma, 


sendo os perturbadores elétrons, podemos escrever o operador 
Hrr dado pela expressão (5.2.24) como: 


US E 
He =-i21N fab b fav v F(y) e qe É 
(0) 0 
(5.4.4) 


A semilargura We e o deslocamento Se eletrônicos, 
medidos em hertz, são dados por: 


E É Es 
,=N Jdbb [dy v F(v) - cos(dr) € ? | (545) 


0 0 


S.=-N | db b ] dv v F(v) sento Eu (5.4.6) 
0 0 
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Para calcular as funções [ir e dir, precisamos calcular a 
transformada de Fourier G(w) do potencial de interação V(t) 
entre o elétron e o átomo emissor, o que faremos a seguir. 


5.4.1 Cálculo de G(qw) 
Para calcular G(mw), definido pela expressão (5.3.21), 


vamos considerar o potencial V(t), entre um elétron e átomo 
neutro, escrito como uma expansão multipolar, ou seja: 


V(y=-e> o E Ta (5.4.1.1) 


em que r é a posição do elétron em relação ao centro de massa 
do átomo, db o multipolo de ordem / definido como [Cattani 
(1968), Bassalo (2006)]: 


47 AM vm La! 2 
Es [pr Yº (0,p) dv”, (5.4.1.2) 


em que r”, 60º, q” são as coordenadas das partículas 
constituintes do átomo medidas no centro de massa deste, e 
Y”(0,9) é o complexo conjugado do harmônico esférico 
[Jackson (1992)]. 

Desse modo, a função G(w) fica escrita na seguinte 
forma: 


G(0)= Sa e e Y (69) q bi (5.4.1.3) 


tH 
r 
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Levando em conta somente as interações de dipolo e 
quadrupolo, teremos: 


G(0) = Gaipolo(00) + Gquadrupolo(0D) + .... (5.4.1.4) 


Consideremos a trajetória do elétron de velocidade v 
como sendo uma reta contida no plano (ZY), conforme mostra a 
Figura 5.4.1. 


e v(t) 


Figura 5.4.1. 


A Figura 5.4.1 indica que a distância r(t) entre o elétron 
e o átomo é dada por: 


0=D+y)= pe BL (5.4.1.5) 
y 
em que b é o parâmetro de impacto da colisão, e: 
vt 
(pv)? 
(5.4.1.6,7) 


cos0= send= 


EN 
(b? pi o ? 
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(a) Interação dipolar 


O termo / = 1 na expressão (5.4.1.1), correspondente à 


interação dipolar, será dado por [lembrar que m = —/,...,0,...,4 
e Y"=(-1 ( X2) em que (+) significa complexo 
conjugado]: 


Y(0,7/2 
Vaipolo(D = Ef é » e dT T” 


EA [T cu co E (T+ T7!) El 


2 


Portanto: 


+oo 0 1Ot 
in eds Tibe a! VE 
pe 


que resulta, efetuando as integrais (vide Apêndice C5-D: 
2ew ob 2 - ob 
Oraipoto(OD =— eo e (e) Zerterço e] 


(5.4.1.9) 


ou mudando de variável q) para Z = wb/v: 


Gapool(Z)=Gi (ZTE + G5(Z) o TOR OR 
(5.4.1.102) 
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com: 


G(D=- p2KD. GD=- SE ZKD 
(5.4.1.10b,c) 


(b) Interação quadrupolar 


O termo / = 2 na expressão (5.4.1.1), correspondente à 
interação quadrupolar, será dado por [lembrar que m = 
Ca APCP À PRP “A o DÁ =(-Dº(yr ), em que (+) significa 
complexo conjugado): 


47 Rm x 
Visa unia (1) = e Se e Y (0,712) 1, = 


=-e E - [5 cos'o-1] -1 - dorm RA) — 


3 sen”0 E 
«E 3 +17) (5.4.1.11) 


Assim, Gquadrupolo(0D) fica: 


—+oo . 
BERRO RL E 
G quadrupolo (0)= rea 2 (bt)? H 


= 1Ot a E tt?e iQt 
ENE by (ri+r7)[ dite Ao vfrivT?) ate 


uv 4 EU as 


(5.4.1.12) 
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Efetuando as integrais (vide Apêndice C5-I) e mudando 
a variável q para a variável Z = wb/v, virá: 


Gquadrupoo(D=GZT; + GHZ) (T+T5) + GHZ) (TI+TD)., 


(5.4.1.13a) 
em que: 
UD=- — € 2 = 
GD= Es [É RAD - 2640], 
q Tue. e 2 2 
GHZ)= E 3 2 K(Z), 


Rs: 2 [52 
GD-— <|z Ko(D + ZK(D). (54.1.13b-d) 


5.4.2 Cálculo das funções Mr e Orr 

As funções Tx e dx (com K = I,F), como já vimos no 
item 5.3, dependem dos elementos de matriz InIG(m)lk>Ê, de 
acordo com as expressões (5.3.23b,c). Assim: 


I<nlG(O)Ik>Ê = I<nlGaipoto(0) + Gquadrupoto(0D) + ... Ik>I? = 
p q p 


= <nIG,(0) + G(DIk>-<kKIG;(0) + Gi(9)n> = 


<nIG(Olk>:<kIG;(0)In> + <nIGg(m)lk>:<kIGi(0)In> + 


+ 


<nIG(Dlk>-<kIGi(m)In> + <nIG (Olk>-<kIG/(0In>. 
(5.4.2.1) 
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Os estados |n> e |k> escritos explicitamente em 
função do momento angular total j e de sua projeção m são 
dados por: 


In>=Ing jumo>, lk>=In juemç>. 


Fazendo uma média sobre as projeções m do estado 
|k>, as funções I'k e 4x tomam os seguintes aspectos: 


2 
o 2a) l : 
Ex =| hi 2jç+1 2a nan Sn jam, >, 


(5.4.2.2) 


pes plam Kengo 1600) Im jpm, > 
E 


2Ju tl SE EE Oo O 


(5.4.2.3) 


Os elementos de matriz de G(w) que aparecem nas 
expressões acima são calculados a seguir, utilizando as 
expressões (5.4.1.4), (5.4.10a-c) e (5.4.1.13a-d) que são dadas 


em função do operador momento multipolar T/”. 
Vamos, agora, deduzir uma equação auxiliar no cálculo 


desses elementos de matriz, considerando os operadores T)” e 


T;" e calculando a seguinte somatória: 


nN+ ' 
pa ada! TE Ini < mam (O) Ja >. 
mamk 


(5.4.2.4) 
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Lembrando que fipjs CD" To” e usando o 


teorema de Wigner-Eckart [Rose (1967), Bassalo e Cattani 
(2005)], a expressão acima ficará: 


> ED" <ndanT In jm, ><n 


Mom 


. —-m' . NE 
aaa! T,: NJaN> rr 


e CD” <jtmmijm><jm—mljm,> x 
MaMk 
x <A JoL ln, Je ><, JT lin 


oJo? T 


ES ente 


Mamk 


<Jj/mynlm,><j!-mp—mlj,—m,> = 


21H 
ZjoH 


x <A Jo T ln) ><m Tn Ja> = 


— [Di 
2j 


DDS <itmnlm>< MM > X 


MoMk 
x<n Tn << dT lo, ju> 


De posse desse resultado, a expressão (5.4.2.4) será 
escrita na forma: 


N+ 
ria 
DX <a iT Inçjem,><n jam] (17) Inçjodny> = 


MoMk 


=A(GeJo!) <a DT nje><nD T nada > 0,7 Omm'- 
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Usando a definição das delta de Kronecker 6, e ó 


(valores nulos para /4/'e mm e unitários quando são 
iguais), a expressão acima ficará [Rose (1967)]: 


N+ 
5 <a oa em >< jam (To In 


Momk 


la = 


= 5 nda In jam >P= 5 <jmemlim ><a jadT An je 
mom MoiMk 
RA 
= 5 Ijtmyemliam> kn jdT nd = E Km ddT nd. 


MoMk 


Quando aplicamos essas relações na expressão (5.4.2.1), 
verificamos que: 


. . . + . 
>» <a) In jam, >-< jm IGG) Ini > =0, 
Dik MM, 


Mk Memo 
Portanto sobram apenas os termos G(0)Gj;(0) e 
G(DG,(0) na expressão (5.4.2.1), que passamos então a 


calculá-los, usando os valores de Ga(w) e Gy(0) dados pelas 
expressões (5.4.1.10a) e (5.4.1.13a). 


Para o termo dipolar, teremos: 
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> kndoamiSA(Dn dm P = 5 kn SUDT + 


Dkjkmkma, Nkjkmkma, 


+ Gigi! +17 Onde > s 


2 
E ninja 5 [GD] encima ciano + 


Dejo, Demo, 


2 
AGIDO [iam dan + nim iam) = 


: 2 
z Era [e SADEn indo, 
DkJk 
(54.2.5a) 
em que: 
A(D=2"|KÍ(Z) + KID) |. (5.4.2.5b) 


Por sua vez, o termo quadrupolar tomará a seguinte 
forma: 


> kndanS (Dna? = 5 kn ijançISUZDT+ 
Mem Meia 


+ 642) [D+T; | + OZ [D+TZ [In jam > = 
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ue (7 0h: 2, GUD2 
= 5 SAD kona Ini, > + GIZ) x 
Dk eMkMa, 


g 1 A 2 ) —1 . pj 
x [nino nim, >É + kndJomT, In Jim > | F 


+ GU nada Im en AH nada! TS Im je > | 


Pig 
415 


E 2 A (DknçaTolny >P, (5.4.2.6a) 
nv Nkjk 


em que: 
AD=L [KU7+4K ID +3K7D)]. (5.4.2.6b) 


Em vista dos resultados obtidos acima [expressões 
(5.4.2.5a,b,6a,b)], a soma dos elementos de matriz G(w) será 
dada por: 


2Ju + 
DD lenda Sn jam, = a ( À :)Z A(Z)x 


Mk Memo TVD Pnnik 


2Jj + E 
x I<njoTln, je > + da [+] SA ncia | 


Nkjk 
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Substituindo a expressão acima nas expressões 
(5.3.23b,c), obtemos I'k diretamente. Contudo, o valor de dk 


será conhecido depois do cálculo da integral indicada abaixo: 


pf dep! Poda SO nie, Eua 
Ok — WO 


—oo 


— Pjutl E, 5 o. 
= Ha (e eng lT Indo «e SEo dE 


“dz A(Z) 
Za Z 


q 2 


2 
E E) I<n, ju lT, IndetP[ 


A parte principal dessas integrais vale [GBKO (1959), 
Klarsfeld (1970), Deutsh e Klarsfeld (1973)]: 


A RA 1 =12[Ko(DI(Z)-K(DI(Z]=B(Z) , 


TdZ AZ) 
= a Z,-Z = Z[B(Z)+b(Z)-1K (21 (D)]=B, (2) 


em que: 


b(Z) = + 7 - 1 ZK(DI(D . 
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De posse das expressões acima, as expressões finais para 
as funções 'k e O serão dadas por:: 


2 
Te= a), E nao lin je 2 A(Z) + 
Nkjk 


PN 

e 

Es is E ] DE KnodalHT, Im, je É Ag (2) , (5.4.2.8) 
nkeik 


2 
de =. [e] DD KnçiallT lin jo >P B(Z) + 


nkjk 


5) 2 

2 

+B(Ei ] DX kndalT In j>ÉBD » 
nkjk 


(5.4.2.9) 


emque T=T/ay e G=T SER sendo a, o raio de Bohr. 


0º 


As funções A(Z), B(Z), Ay(Z) e Bu(Z) comportam-se de 
acordo com a figura 5.4.2 [GBKO (1959)]: 
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mm DIPOLO 
——— QUADRUPOLO 


Figura 5.4.2 


A seguir, vamos calcular os alargamentos e os 
deslocamentos para as linhas emitidas pelo átomo de hélio 
neutro. 


5.5 Equações para o Cálculo do Alargamento e do 
Deslocamento de Linhas Espectrais do Nélio 
Neutro 


Quando o emissor é o átomo de hélio neutro, os estados 
inicial I e final F envolvidos na transição são definidos por: 


= a 
IK> = Inç/ m,>, 
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com K=IouF,emque a = 1 parao caso do para-hélio e a = 
3 para o orto-hélio. 


Os estados virtuais | k> envolvidos na determinação do 
alargamento e do deslocamento são definidos por: 


k> = In“ m'>. 


Assim, as funções Ix e dx, com K = [ou F, dadas 
pelas expressões (5.4.2.8,9), tomarão as seguintes formas: 


2 2 a api 
E s(5 a sp tentelmiInte> 


Ko 3 (by) & a? 


22 E a a gi 2 
E a ste Alim TAS AD. Es 


np! dç 


2 É a agi 
—2/ea d<nC AT In“ 2 
dx= E! [e] A oO. B(Z) + 


dao a mpi 
2 ed den AT In > 
+ B(Z), Dez 
15 ES z ay at2) E 


em que Z=(b/v) O 


nf ndp E 


Podemos mostrar que a contribuição de quadrupolo nas 
funções I'x e Ox é desprezível para o átomo de hélio neutro, 
calculando a razão entre os elementos reduzidos de matriz da 
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contribuição dipolar e quadrupolar, a partir das expressões 
(5.5.1,2). Com efeito: 


Co, E: > ken“am in" o'>P A(Z)/a5 

np 

(as a D kn“T ln" >A (Z) /ag 
ny 


kd 


> ienamIn? 0'>P B(Z)/a? 


(dx), E 5b” n'dy! 
(00), as D kn “aTlIn 0>PB(Z)/a0 
no 


Os elementos reduzidos de matriz de dipolo e 
quadrupolo são dados pelas expressões indicadas abaixo [Oertel 
e Shomo (1968)], quando as transições virtuais ocorrem 
mantendo o número quântico principal nº, ou seja, nºº = na = 
n: 


an “(nº mibr )—— para !'=/—1 


ã a 
knº(ULln0'>?= 
do 


pie 
ER 


2-(4+1) a. para (=/+] 


T 
kn/-2|n?/'>2= 
4 
E) 
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Efbê- fo , paral=(-1 


Doe pi-()s para(= (41. 


Das expressões acima podemos observar que quanto 
maior o valor de / menores são os elementos reduzidos de 
matriz de dipolo e quadrupolo. Portanto, substituindo os 
elementos reduzidos de matriz nas equações acima, passamos a 
analisar apenas a relação (T'k)a/([k)a. Levando em conta que, 
em média (vide Figura 5.4.2) A(Z) - Ad(Z), obtemos: 


(To, > 9 b? 
(De), 25 ão 


x 


: (ia tn (un 


af s fio np MA af ns 27 (420) 
n [(n Yu) )| e a [é | [nen] Ro 


Uma estimativa pode ser feita, considerando os casos 
limites para /,istoé, (=0 e /=n- 1. Desse modo, teremos: 


Numericamente, observamos que essa razão é da ordem 
de - 102 b2/aj e como b/ao - 1 na região de forte interação, 


resulta que: 


Do, = 02 
CTA 


e, portanto, podemos desprezar a contribuição de quadrupolo. 
No caso do deslocamento da linha, de modo semelhante 
podemos mostrar que a contribuição quadrupolar é desprezível 


comparada com a dipolar. 


Do que vimos acima, para o hélio neutro, podemos 


escrever que: 


2 2 
ea a “Ip tã 01 
Rea [e] Den atT in e> A(Z), (5.5.3) 


nã 


2 ea i * 2 
dx= 512] Dent M/lin'te>P B(D À (5.5.4) 
ny 
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Desse modo, o alargamento e o deslocamento da linha 
do hélio neutro são obtidos, então, substituindo-se as expressões 
de IK e dx, dadas pelas expressões (5.5.3,4) nas expressões 
(5.4.5,6). 


5.6 Alargamento e Deslocamento Provocados por 
Elétrons e Ions 


O efeito dos íons na aproximação quase estática foi 
levado em conta na função de correlação temporal obtida pela 
equação (5.2.26), resultando na forma de linha dada pela 
equação (5.2.28). 

Griem (1962, 1968, 1974a,b) mostrou que, quando os 
íons são tratados na aproximação quase estática, como foi feito 
no item 5.2, o alargamento total w e o deslocamento total S da 
linha, assim obtida, são muito bem representados pelas 
expressões abaixo, como a soma dos alargamentos Qe e Q e 
deslocamentos Se e Si eletrônicos e iônicos, ou seja: 


0=0,+0 =[1+1.75A (1-0.75R) | 00, (5.6.1) 
S=S,.+S,=S, + 2.04 (1-0.75R) 0,, (5.6.2) 
em que: 
Cc E? 3/4 
des [Eai Do R=Pm, (5634) 
e O, Pp 


sendo C4, Fo e Pp, definidos no item 5.2 e pp é o raio de 


Debye. Os índices e e i representam as contribuições geradas 
pelos elétrons e íons, respectivamente. 
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As expressões (5.6.1) e (5.6.2) são válidas para valores 
de R e de A, obedecendo aos limites: R < 0.8 e 0.05 < A < 
0.5. A limitação no parâmetro A é feita inferiormente, pois para 
A < 0.05 a interação quadrupolar não é desprezível, isto é, essa 
condição implica que a intensidade máxima Fo do campo 
estático dos íons é pequena e portanto, apesar de o efeito de 
quadrupolar ser pequeno, ele começa a ficar da ordem de 
grandeza do efeito dipolar, não podendo ser desprezado. O 
limite superior de A = 0.5 ocorre quando a intensidade máxima 
Fo do campo eletrostático dos íons é grande, aparecendo o efeito 
Stark linear. Nestas condições, a linha não é mais isolada, 
aparecendo a superposição das linhas permitida e proibida. 

Nas expressões (5.6.1) e (5.6.2) o alargamento e o 
deslocamento eletrônicos podem ser calculados pelas teorias 
existentes. Em particular, usaremos a nossa Teoria 
Semiclássica Convergente desenvolvida no item 5.4. 

A utilização das expressões referidas acima para o 
cálculo de alargamento e deslocamento iônicos depende da 
condição de validade da aplicação da aproximação quase 
estática. Nessa aproximação deve estar satisfeita a condição de 
que o tempo de colisão iônica é muito maior que o tempo de 
colisão eletrônica, o que pode ser verificada por meio do 
parâmetro 6 >> 1, sendo esse parâmetro definido como: 


— DPm 
OR Te (5.6.5) 


ion 
Assim, enquanto o elétron colide com o átomo (emissor), 
o íon está praticamente parado e, consequentemente, utiliza-se a 
distribuição de Holtsmark, para o campo elétrico E no interior 
do plasma, conforme fizemos no item 5.2. 
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Uma questão importante a ser observada é que a 
introdução da contribuição iônica, segundo a aproximação quase 
estática, provoca uma assimetria na forma da linha espectral, no 
mesmo sentido em que ocorre o deslocamento desta, quando é 
feita a convolução do efeito das colisões iônicas com a 
Lorentziana que descreve a forma da linha gerada apenas pelas 
colisões eletrônicas. 

Quando a condição de validade 6 >> 1 de utilização da 
aproximação quase estática não for preenchida, as expressões 
(5.2.29), (5.6.1) e (5.6.2), que representam, respectivamente, o 
perfil, o alargamento e o deslocamento da linha espectral, devem 
ser substituídas por expressões que levem em conta o efeito da 
dinâmica dos íons. 

Para o < 1, GBKO (1959) e Barnard, Cooper e Smith 
(1974) propõem expressões para o perfil da linha e seus 
respectivos alargamentos e deslocamentos que levem em conta o 
movimento dos íons. 

GBKO (1959), por exemplo, obtém um perfil reduzido 
para o caso em que a velocidade do íon (Vion) é grande, dado 
por: 


i(x,0,6) = | R far exp/ixt=2 > (5 r(Bos ot 
io t É 4142 3 l 


(5.6.6) 


em que x, O, e O foram definidos pelas expressões (5.2.30a,b) e 
(5.6.5). 

Desse modo, o alargamento e o deslocamento totais da 
linha nessa aproximação são obtidos por meio da expressão do 
perfil reduzido, como [Griem (1962)]: 


o=[1+1300"]m,, (5.6.7) 
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S=S,.+2360%º 0" 0,, (5.6.8) 


Por sua vez, o perfil reduzido obtido por Barnard, 
Cooper e Smith (1974), usando a aproximação adiabática, é 
dado por: 


1+ 
i(x,0,0) = À. RE pá (5.6.92) 
T (xd + (l+w) 
em que: 
“0 89 Bull 
w= 7 =000w (k), (5.6.9b) 
e +/oto o D! (k) — 304 Pa, (5.6.9€) 
O, Po 
k=1.03 06º 02º x, (5.6.9d) 


e Wik) e DI(k) são funções tabeladas [Barnard, Cooper e 
Smith (1974)]. 


Uma análise da expressão (5.6.9a) mostra que, nos 
limites Xu > +00 e x/y — —o, o perfil comporta-se 
segundo as expressões: 
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e no centro da linha x > O o perfil j(x, o, 6) se aproxima de 
uma Lorentziana cuja largura e deslocamento são dados por: 


(5.6.10) 


e 2 


o = É + 1.36 08º RA 0) 


S=S. + 2.56 08º 618 — 3048 Ba o. (5.6.11) 


0 


As duas aproximações expostas para o caso em que 6 < 1 
levam a expressões idênticas para o alargamento, diferindo 


apenas no deslocamento da linha pelo termo -3 0º? (p. /po)o,, 


que corresponde à introdução do “screening” de Debye. 

Uma outra forma de levar em conta o efeito dos íons é 
tratá-los na aproximação de impacto, válida quando estiverem 
satisfeitas as condições. 


o<1 e im is (5.6.12a,b) 


5.7 Aplicações 


Algumas aplicações das expressões deduzidas neste 
Capítulo podem ser vistas em trabalhos realizados pelos autores 
e colaboradores. Por exemplo, em 1980, Bassalo, Cattani e 
Walder apresentaram o cálculo da forma (semilargura e 
deslocamento) de várias linhas espectrais isoladas Lorentzianas 
(em Á: 7281, 7065, 6678, 5876, 5048, 5016, 4713, 4121, 3965, 
3889, 3188, 2945) do hélio neutro (He [) em um plasma quente. 
Enquanto a contribuição eletrônica (com e sem “screening” de 
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Debye) para aquela forma foi calculada por intermédio do 
formalismo convergente, desenvolvido por Cattani (1972), e 
Yamamoto e Cattani (1978), e apresentado neste Capítulo, a 
contribuição iônica foi obtida usando-se o formalismo 
semiclássico de GBKO (1962). Além do mais, os elementos de 
matriz reduzidos foram calculados usando-se os valores 
“oscillator strenght” tabelados por Wiese, Smith e Glennon 
(1966) e a diferença de energia h la, dados por Moore 


nã! 
(1949). Os resultados teóricos são comparados com resultados 
experimentais, conforme se pode ver nas tabelas apresentadas no 
Apêndice C5-II. 

O alargamento Stark da componente permitida da linha 
3965Ã do hélio neutro foi calculado por Araújo, Bassalo, 
Cattani e Walder (1981). Como no artigo anterior (1980), neste 
de 1981, usamos o formalismo convergente para obter a 
contribuição eletrônica, e o formalismo semiclássico de GBKO 
(1962) para a contribuição iônica. Nossos resultados teóricos 
foram comparados com os resultados teóricos e experimentais 
de Pilon e Barnard (1980) e com os resultados experimentais de 
Diatta (1977). Ainda usando esse mesmo formalismo, Bassalo, 
Cattani e Walder (1981) analisaram o perfil Stark e as larguras 
da linha SOI6Á do hélio neutro. Nossos cálculos teóricos 
mostraram um bom acordo ao serem comparados com os 
resultados experimentais de Chiang, Murphy, Chen e Griem 
(1977). 

O formalismo convergente foi também usado por 
Bassalo, Cattani e Walder (1982) no cálculo das larguras (0,) e 


deslocamentos (d,) eletrônicos de 42 linhas espectrais isoladas 
Lorentzianas (em Á: 584, 537, 522, 20581, 5016, 3965, 7281, 


5048, 4438, 6678, 4922, 4388, 15084, 11013, 9603, 19089, 
12968, 11045, 18697, 10830, 3889, 3188, 2945, 2829, 7065, 
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4713, 4121, 5876, 4472, 4026, 12528, 9464, 8362, 21120, 
12846, 10668, 17002, 11969, 19543, 12985, 10997, 18686) do 
hélio neutro (He 1) para uma densidade eletrônica de 10'º /cm” 
e temperaturas T = 5.000, 10.000, 20.000 e 40.000 K. 
Além disso, nesses trabalhos foram também calculados os 
parâmetros A e R definidos por Griem (1974a,b). Nossos 
resultados estão apresentados no Apêndice C5-III. 

É oportuno registrar que os nossos resultados para as 
larguras e deslocamentos das linhas citadas acima são usados até 
hoje [Pelaez, Gonzalez, Rodriguez, Aparício e Mar (2006), 
Omar, Gunter, Wierling e Ropke (2006), Bartecka, Wujec, 
Halenka e Musielok (2004), Milosavljevic e Djenize (2003, 
2002a,b, 2001), Pandey, Rao, Lambert, Jeffery e Asplund 
(2001)] por pesquisadores do mundo inteiro que trabalham em 
diagnóstico de plasmas utilizando as transições do hélio como 
referência. 


Apêndice C5-I 
CÁLCULO DAS INTEGRAIS DA FUNÇÃO G(o) 
(mv) eG (00), dadas 


pelas expressões (5.4.1.10a-c) e (5.4.1.13a-d), envolve integrais 
do tipo: 


O cálculo das funções G 


dipolo quadrupolo 


DO q dte'” — te dt(cosat+i senot) 
L, o IN (b? pad no = [m (»? e poe ? (All) 
moço dit" e? 
e com m=1€2. (A1.2) 


Ex (b? 4 vp A ” 
1) Integral L, 


Tomemos a expressão (Al.l). Como a função 
envolvendo sen wté ímpar, a sua integral no intervalo 
(— oo + co) é nula e, portanto, resta apenas a integral envolvendo 


cos wt, que é uma função par. Portanto, a expressão (Al.1) 
tomará a forma: 


n+1/2. 


E E dt e'” : ir dt cosat 
Ed es n+1/2 ti 0 (b? +vr) 


(b? + VP) 


a: ie dt cosat 


de (A13) 
0 (22+4) 1/2 
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onde Z, =b/v. 


Por outro lado, sabemos que [Gradshteyn e Ryzhik 
(1965)]: 


D(n+1/2) (22) (+= dtcosot T(n+1/2) (Qn)! 
Pupo Jo (240) Ta?) 2º n! 


0 


K (0Z,) = 


Levando as expressões acima na expressão (Al.3), e 
chamando Z = wb/v, teremos: 


92141 o n! 


vi Qm! Cb 


 K,(Ob/v) = 


qu+l n! 
[= CARO: (AL4) 
v(2n)!b” 


Agora, estudemos os valores particulares de n. 
la)n=1 


Neste caso, a expressão (A 1.4) tomará a forma: 


(AÍ.5) 
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lb)n=2 


Neste caso, a expressão (A 1.4) tomará a forma: 


ate ARO AD (A1.6) 


2) Integral J” 


Tomemos a expressão (A1.2): 


ias dt t” ei? 
fi (5º + RR com m=le2. 
v 


JP = 


Agora, estudemos os valores particulares de m. 


2a) m=1 
n Sir diteo. «Ta E diet 1d 
nO Jo (pve)? E i do + (pve)? E i do. 
(A1.7) 


Usando a definição de Z (Z=wblv) e a expressão 


(A1.4), a expressão acima ficará: 


A e dEte o . e do Ao dio 02... 
a (prvi) i do i dZ do 


Dl d 


Evo le KZ] (A1.8) 


Considerando que [Gradshteyn e Ryzhik (1965)]: 
ag K,(Z)|=nZ" K,(2)+2" K/(Z) 
dZ n n n 2 


onde a linha (”) indica a derivada em relação ao argumento (Z), 
E: 


2K(D)=K, (2)+K, (2) € > K(Z)=-K (2)+K (Z) 


K(D)=- = KlZ) +K, (2) = K,(Z)|=-Zº K, (2) 


(A1.9) 


Substituindo a expressão (Al.9) na expressão (Al.8), 
teremos: 


o iot n+1 | 
T, =] aa ENE Ia E — Zº Ku(Z). 
“(brvr) ivó b” Qn)! 


(A1.10) 
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Agora, estudemos os valores particulares de n. 
2a.1)n=1 
Neste caso, a expressão (A1.10) tomará a forma: 


pet Ca) (ALII) 


"(bavryCo iv 
2a.2)n=2 
Neste caso, a expressão (A 1.10) tomará a forma: 


RN dt te'” o 2 ã 
da rr RE RR o 


2b) m=2 


Neste caso a expressão (A 1.2), assumirá a forma: 


» qto dtte? 1d o dte' d'I 
J, =[ 2 SA O 2 IN E 99H O a" 
pb] rdo o (p+vp) do 
Usando a definição de Z (Z=wb/v) e as expressões 
(A1.4,9), a expressão acima ficará: 
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P si ge o cod [) dz. 
e fat (be) do dZl do) do 
--bdfdi, dy) bl dd, bl dadi, 
v di dz do vº dZ dz v? dZ dz 
28 fears 
EaD mo  molaoZ  K(Z)] E 
vº (Qn)!b dZ A dZ 
te dit? e! Rol dorso, 
1 = IB nHID (3 22 [Z K(Z)]. 
(+ vr) v* Qn)!b dz 
(A1.13) 


Considerando que [Gradshteyn e Ryzhik (1965)]: 


d n — n n-1 
o [7º KulZ) |= 2º Ko(D) +72" Ka(D), 


a expressão (A1.13) ficará: 
dt t? ei? qu+l n! 


2 — qe: — n n-1 
J, = Í. (b? +v'P js E v (2n)! bp? IZ K,(Z) +Z K,1(Z2)] 


(A1.14) 
Agora, estudemos os valores particulares de n. 


2b.1))n=1 


Considerando que K (Z)=K,(Z) [Gradshteyn e 
Ryzhik (1965)], a expressão (A1l.14) tomará a forma: 
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=2 PERDA KMD. (ALIS 
V 


3 


2b.2)n=2 


Neste caso, a expressão (A1.14) tomará a forma: 


2 E dt t? ei? E 3) ; 
n=[. rapa str É Ko) + Za 


(A1.16) 


Deixamos ao leitor a demonstração das expressões 
(5.4.1.10a-c) e (5.4.1.13a-d), usando os resultados obtidos neste 
Apêndice. 
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Apêndice C5-II 


CÁLCULOS COM A TEORIA CONVERGENTE DO 
ALARGAMENTO E DO DESLOCAMENTO DE LINHAS 
ESPECTRAIS ATÔMICAS (CONTRIBUIÇÕES 
ELETRÔNICAS E IÔNICAS) DO HÉLIO NEUTRO 
EM UM PLASMA 


Neste Apêndice vamos apresentar os resultados da forma 
de linha espectral (largura e deslocamento) de algumas linhas do 
hélio neutro em um plasma quente. Os cálculos da contribuição 
eletrônica (com e sem “screening” de Debye) foram realizados 
por intermédio do formalismo convergente, desenvolvido por 
Cattani (1972), e Yamamoto e Cattani (1978), enquanto a 
contribuição iônica foi obtida usando-se o formalismo 
semiclássico de GBKO (1962). 


TABELA I. Comparação das semilarguras medidas e calculadas 
para algumas linhas do hélio neutro 


(A) 
o 
Semilarguras (A) 
Comprimento Calculada Medida 
de onda Densidade Sem Debye Com Debye 
ú eletrônica Temperatura R A “screening” “screening” 
(ema)  (0%%em)  (10K) 

7281 .03 20.9 0.29 0.076 0.80 0.80 0.88º 
7065 .03 20.9 0.29 0.062 0.46 0.46 047 
6678 .03 20.9 0.29 0.190 0.78 0.78 0.98º 

5876 16.0 49.0 0.30 0.162 4.50 4.50 5.5? 

13.0 43.0 0.31 0.152 3.68 3.68 49º 

41 20.0 0.37 0.104 1.22 1.22 
2.0 18.0 0.35 0.088 0.58 0.58 

.03 20.9 0.29 0.076 0.28 0.28 0.39º 

5048 3.2 30.0 0.29 0.173 Seo 5.da 4.6º 

2.0 18.0 0.35 0.157 3.10 3.10 3.4º 
.03 20.9 0.29 0.132 1.58 1.58 1.68º 

5016 3.2 30.0 0.29 0.257 2.48 2.44 1.9 

17.0 24.0 0.43 0.409 14.90 13.36 13.0" 

41 20.0 0.37 0.262 3.36 3.24 3.84 

2.0 18.0 0.35 0.213 1.60 1.58 184 
.03 20.9 0.29 0.184 0.78 0.78 0.86" 

4713 3.2 30.0 0.29 0.138 2.98 2.98 29º 
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13.0 20.0 0.45 0.202 12.06 12.06 14.0º 
41 20.0 0.37 0.151 3.66 3.66 4.6 
2.0 18.0 0.35 0.128 1.72 1.72 2.054 
03 20.9 0.29 0.107 0.90 0.90 0.95* 
4121 3.2 30.0 0.29 0.189 8.12 8.12 6.2º 
2.9 18.0 0.37 0.190 6.98 6.96 7.0! 
03 20.9 0.29 0.145 2.42 2.42 2.22? 
3965 4 18.0 0.33 0.395 3.34 3.14 ai 
03 20.9 0.29 0.353 2.38 2.28 2.30 
3889 3.2 30.0 0.29 0.103 0.70 0.70 0.73º 
15.0 26.0 0.41 0.151 3.46 3.46 4.5º 
41 20.0 0.37 0.109 0.92 0.92 1.25º 
2.0 18.0 0.35 0.093 0.42 0.42 0.554 
03 20.9 0.29 0.078 0.22 0.22 0.24 
3188 15.0 29.0 0.38 0.289 11.74 11.32 13.4? 
41 20.0 0.37 0.202 3.00 2.98 3.0" 
q 18.0 0.34 0.161 1.20 1.20 1.35º 
03 20.9 0.29 0.143 0.72 0.72 0.78 
2945 03 20.9 0.29 0.220 1.98 1.98 2.0º 


* Kelleher (1977). 

* Berg, Ali, Lincke, Griem (1962) 
“Wulff (1958) 

“ Diatta (1977) 


TABELA I. Comparação dos deslocamentos medidos e calculados 
para algumas linhas do hélio neutro 


o 
Deslocamentos em (A) 


Comprimento Calculado Medido 
de onda Densidade Sem Debye Com Debye 
o. eletrônica Temperatura R A “shielding” “shielding” 
(emA) (0tcm*) (IO K) 
7281 03 20.9 0.29 0.076 +0.37 +0.36 +0.473 * 
7065 .03 20.9 0.29 0.062 +0.23 +0.23 +0.302 * 
6678 03 20.9 0.29 0.190 +0.22 +0.19 +0.355 * 
5876 16.0 49.0 0.30 0.162 -1.37 -1.22 +0.7º 
13.0 43.0 0.31 0.152 -1.13 -1.02 0.0º 
41 20.0 0.37 0.104 -0.41 -0.38 -0.309 
2.0 18.0 0.35 0.088 -0.20 -0.19 -0.10º 
03 20.9 0.29 0.076 -0.097 -0.097 -0.074º 
5048 3.2 30.0 0.29 0.173 +2.19 +2.05 +2.0º 
2.0 18.0 0.35 0.157 +1.43 +1.35 +0.90º 
03 20.9 0.29 0.132 +0.71 +0.68 +0.891º 
5016 3.2 30.0 0.29 0.257 -0.83 -0.71 -0.55º 
17.0 24.0 0.43 0.409 -5.14 -3.91 -4.8? 
41 20.0 0.37 0.262 -1.14 -0.94 -0.854 
2.0 18.0 0.35 0.213 -0.54 -0.46 -0.40º 
03 20.9 0.29 0.184 -0.25 -0.23 -0.283º 
4713 3.2 30.0 0.29 0.138 +1.34 +1.28 +1.4º 
13.0 20.0 0.45 0.202 +5.98 +5.37 +6.0º 
41 20.0 0.37 0.151 +1.81 +1.71 +2.3º 
2.0 18.0 0.35 0.128 +0.87 +0.83 + 1.0 
03 20.9 0.29 0.107 +0.44 +0.43 + 0.538º 


4121 3.2 30.0 0.29 0.189 +3.61 +3.31 +2.8º 


2.9 18.0 

03 20.9 

3965 4 18.0 
.03 20.9 

3889 32 30.0 
15.0 26.0 

41 20.0 

2.0 18.0 

03 20.9 

3188 15.0 29.0 
41 20.0 

BIA 18.0 

03 20.9 

2945 03 20.9 


* Kelleher (1977). 

* Berg, Ali, Lincke, Griem (1962) 
“Wulff (1958) 

“ Diatta (1977) 


0.37 
0.29 
0.33 
0.29 
0.29 
0.41 
0.37 
0.35 
0.29 
0.38 
0.37 
0.34 
0.29 
0.29 


0.190 
0.195 
0.395 
0.353 
0.103 
0.151 
0.109 
0.093 
0.078 
0.289 
0.202 
0.161 
0.143 
0.220 


[oro] 


+3.10 
+1.09 
-0.95 
-0.70 
+0.13 
+0.66 
+0.19 
+0.99 
+0.044 
+2.36 
+0.68 
+0.28 
+0.17 
+0.43 
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+1.124º 
-0.739 
0.771? 
+ 0.3? 
+12» 
+0.254 
+0.10º 
+0.080" 
+4.1º 
+0.94 
+0.30º 
+0.258º 


+0.683º 
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Apêndice C5-III 


CÁLCULOS COM A TEORIA CONVERGENTE DO 
ALARGAMENTO E DO DESLOCAMENTO CAUSADOS 
POR COLISÕES ELETRÔNICAS DE LINHAS 
ESPECTRAIS ATÔMICAS DO HÉLIO NEUTRO EM UM 
PLASMA E DOS PARÂMETROS IÔNICOS A E R 


Neste Apêndice são apresentados os resultados dos 
cálculos usando a Teoria Convergente dos alargamentos e 
deslocamentos causados por colisões eletrônicas de 42 linhas do 
hélio neutro em quatro diferentes temperaturas. São também 
apresentados os cálculos dos parâmetros iônicos À e R. 


Tabela 1. Larguras (w.), deslocamentos (d.) e parâmetro Stark A para o He(l) com 
densidade eletrônica N = 10!%/em” e temperaturas de 5000, 10.000, 20.000 e 40.000 K. 


o o o A 
L(A) Me (A) de (A) 
584 0.1 108E-03 -0.1371E-03 0.0126 
0.1290E-03 -0.1328E-03 0.0112 
0.1472E-03 -0.1202E-03 0.0102 
0.1709E-03 -0.1019E-03 0.0091 
537 0.4066E-02 -0.2887E-02 0.1604 
0.3863E-02 -0.2383E-02 0.1667 
0.3548E-02 -0.1904E-02 0.1777 
0.3169E-02 -0.1474E-02 0.1934 
522 0.1624E-01 -0.1089E-01 0.2965 
0.1519E-01 -0.8826E-02 0.3118 
0.1373E-01 -0.6924E-02 0.3363 
0.1205E-01 -0.5269E-02 0.3710 
20581 0.3066E+00 -0.4131E+00 0.0428 
0.3724E+00 -0.4121E+00 0.0370 
0.4524E+00 -0.3812E+00 0.0320 
0.5427E+00 -0.3258E+00 0.0279 
5016 0.3562E+00 -0.2587E+00 0.1610 
0.3394E+00 -0.2168E+00 0. 1670 
0.3138E+00 -0.1770E+00 01771 
0.2832E+00 -0.1406E+00 0.1913 
3965 0.9373E+00 -0.6303E+00 0.2958 
0.8769E+00 -0.5122E+00 0.3110 
0.7932E+00 -0.4035E+00 0.3353 
0.6968E+00 -0.3090E+00 0.3695 
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7281 


5048 


4438 


6678 


4922 


4388 


15084 


11013 


9603 


19089 


12968 


11045 


18697 


5000 
10000 
20000 


0.2754E+00 
0.3178E+00 
0.3519E+00 
0.3698E+00 
0.5347E+00 
0.6064E+00 
0.6528E+00 
0.6622E+00 
0.1265E+01 
0.1406E+01 
0.1488E+01 
0.1481E+01 
0.3954E+00 
0.3576E+00 
0.3128E+00 
0.2683E+00 
0.2254E+01 
0.1923E+01 
0.1598E+01 
0.1300E+01 
0.5443E+01 
0.4620E+01 
0.3810E+01 
0.3069E+01 
0.1375E + 02 
0.1300E + 02 
O.1196E + 02 
0.1073E + 02 
0.2121E+02 
0.1970E+02 
0.1773E+02 
0.1552E+02 
0.3777E+02 
0.3471E+02 
0.3086E+02 
0.2666E+02 
0.3630E+02 
0.3156E+02 
0.2675E+02 
0.2217E+02 
0.4826E+02 
0.4122E+02 
0.3423E+02 
0.2777E+02 
0.8128E+02 
0.6844E+02 
0.5611E+02 
0.4502E+02 
0.1873E+02 
0.1555E+02 
0.1266E+02 


0.3595E+00 
0.3546E+00 
0.3278E+00 
0.2838E+00 
0.6515E+00 
0.6235E+00 
0.5578E+00 
0.4684E+00 
0.1460E+01 
0.1371E+01 
0.1203E+01 
0.9927E+00 
0.2332E+00 
0.1801E+00 
0.1362E+00 
0.1014E+00 
0.8776E+00 
0.6398E+00 
0.4526E+00 
0.3120E+00 
0.1915E+01 
0.1350E+01 
0.9213E+00 
0.6133E+00 
-0.9741E+01 
-0.8162E+01 
-0.6658E+01 
-0.5272E+01 
-0.1340E+02 
-0.1091E+02 
-0.8620E+01 
-0.6613E+01 
-0.2138E+02 
-0.1707E+02 
-0.1321E+02 
-0.9903E+01 
0.1588E+02 
0.1200E+02 
0.8817E+01 
0.6311E+01 
0.1774E+02 
0.1273E+02 
0.8873E+01 
0.6044E+01 
0.2582E+02 
0.1763E+02 
0.1178E+02 
0.7753E+01 
-0.5803E+01 
-0.4125E+01 
-0.2834E+01 


0.0905 
0.0813 
0.0753 
0.0725 
0.1517 
0.1380 
0.1306 
0. 1292 
0.2171 
0.1995 
0.1912 
0.1918 
0.1539 
0.1659 
0.1835 
0.2059 
0.6947 
0.7827 
0.8991 
1.0497 
1.2667 
1.4325 
1.6552 
1.9468 
0.2973 
0.3101 
0.3301 
0.3581 
0.4571 
0.4830 
0.5227 
0.5775 
0.6497 
0.6922 
0.7560 
0.8436 
0.6707 
0.7449 
0.8434 
0.9709 
1.2556 
1.4132 
1.6246 
1.9006 
1.9498 
2.2182 
2.5744 
3.0369 
0.7524 
0.8651 
1.0094 


10830 


3889 


3188 


2945 


2829 


7065 


4713 


4121 


5876 


4472 


4026 


12528 


9464 


40000 


5000 
10000 
20000 
40000 

5000 
10000 
20000 
40000 

5000 
10000 
20000 
40000 

5000 
10000 
20000 
40000 

5000 
10000 
20000 
40000 

5000 
10000 
20000 
40000 

5000 
10000 


0.1013E+02 
0.4046E-01 
0.5128E-01 
0.6718E-01 
0.8896E-01 
0.8537E-01 
0.9425E-01 
0.9863E-01 
0.9836E-01 
0.2672E+00 
0.2891E+00 
0.2959E+00 
0.2871E+00 
0.6823E+00 
0.7357E+00 
0.7458E+00 
0.7139E+00 
0.1563E+01 
0.1649E+01 
0.1639E+01 
0.1542E+01 
0.1532E+00 
0.1798E+00 
0.2055E+00 
0.2253E+00 
0.2923E+00 
0.3391E+00 
0.3780E+00 
0.3989E+00 
0.7728E+00 
0.8905E+00 
0.9807E+00 
0.1018E+01 
0.1343E+00 
0.1377E+00 
0.1338E+00 
0.1253E+00 
0.1360E+01 
0.1188E+01 
0.1015E+01 
0.8495E+00 
0.3605E+01 
0.3129E+01 
0.2649E+01 
0.2187E+01 
0.4110E+01 
0.4495E+01 
0.4657E+01 
0.4585E+01 
0.7176E+01 
0.7716E+01 


-0.1882E+01 
-0.5627E-01 
-0.5718E-01 
-0.5449E-01 
-0.4849E-01 
0.5968E-01 
0.4610E-01 
0.3230E-01 
0.2075E-01 
0.1983E+00 
0.1620E+00 
0.1261E+00 
0.9512E-01 
0.4261E+00 
0.3268E+00 
0.2355E+00 
0.1630E+00 
0.9212E+00 
0.6955E+00 
0.4868E+00 
0.3433E+00 
0.2170E+00 
0.2222E+00 
0.2144E+00 
0.1935E+00 
0.3904E+00 
0.3888E+00 
0.3630E+00 
0.3170E+00 
0.9994E+00 
0.9817E+00 
0.9020E+00 
0.7755E+00 
-0.1157E+00 
-0.9678E-01 
-0.76 15E-01 
-0.5715E-01 
0.7900E-01 
. 0.6972E-02 
-0.2754E-01 
-0.4055E-01 
0.25 18E+00 
0.4626E-01 
-0.6152E-01 
-0.1074E+00 
0.2645E+01 
0.1980E+01 
0.1333E+01 
0.8161E+00 
0.4400E+01 
0.3319E+01 


1.1928 
0.0321 
0.0269 
0.0219 
0.0178 
0.0848 
0.0788 
0.0761 
0.0763 
0.1511 
0.1424 
0.1400 
0.1431 
0.2314 
0.2187 
0.2 164 
0.2236 
0.3157 
0.3032 
0.3046 
0.3188 
0.0758 
0.0672 
0.0608 
0.0567 
0.1284 
0.1149 
0.1059 
0.1017 
0.1726 
0.1552 
0.1444 
0.1404 
0.0722 
0.0709 
0.0725 
0.0761 
0.6144 
0.6802 
0.7653 
0.8746 

1.0564 

1.1749 

1.3312 

1.5370 
0.1468 
0.1373 
0.1337 
0.1353 
0.2274 
0.2153 
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21120 


12846 


10668 


17002 


11969 


19543 


12985 


10997 


18686 
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20000 
40000 
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20000 
40000 
5000 
100º%0 
20000 
40000 
5000 
10000 
20000 
40000 
5000 
10000 
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40000 
5000 
10000 
20000 
40000 
5000 
10000 
20000 
40000 
5000 
10000 
20000 
40000 
5000 
10000 
20000 
40000 


0.7818E+01 
0.7501E+01 
0.1365E+02 
0.1441E+02 
0.1433E+02 
0.1350E+02 
0.6076E+01 
0.7482E+01 
0.8654E+01 
0.9331E+01 
0.7423E+01 
0.8674E+01 
0.9679E+01 
0.1016E+02 
0.1243E+02 
0.1422E+02 
0.1549E+02 
0.1587E+02 
0.2007E+02 
0.1778E+02 
0.1548E+02 
0.1323E+02 
0.3195E+02 
0.2781E+02 
0.2366E+02 
0.1965E+02 
0.1034E+02 
0.1126E+02 
0.1157E+02 
0.1123E+02 


0.1357E+02 


0.1460E+02 
0.1480E+02 
0.1418E+02 


0.2363E+02 


0.2496E+02 
0.2483E+02 
0.2340E+02 


0.1102E+02 


0.9227E+01 
0.7622E+01 


0.6219E+01 


0.2324E+01 
0.1539E+01 
0.8007E+01 
0.6021E+01 


0.4266E+01 
0.2904E+01 
0.6114E+01 
0.6236E+01 
0.6011E+01 
0.5416E+01 
0.9304E+01 
0.9117E+01 
0.8383E+01 
0.7229E+01 
0.1528E+02 
0.1472E+02 
0.1325E+02 
0.1665E+02 
0.4487E+00 
-0.5750E+00 
-0.9799E+00 
-0.1040E+01 
0.1987E+01 
0.1455E+00 
-0.7953E+00 
-0.1160E+01 


0.7777E+01 


0.6317E+01 
0.4818E+01 
0.3518E+01 


0.8659E+01 


0.6698E+01 
0.4883E+01 
0.3430E+01 


0.1403E+02 


0.1064E+02 
0.7652E+01 
0.5329E+01 


-0.1924E+01 


-0.1120E+01 
-0.6090E+00 


-0.295 1E+00 


0.2132 
0.2199 
0.3151 
0.3026 
0.3038 
0.3177 
0.1009 
0.0864 
0.0774 
0.0732 
0.1682 
0.1496 
0. 1378 
0.1329 
0.2454 
0.2217 
0.2080 
0.2042 
0.6020 
0.6593 
0.7313 
0.8228 
1.0529 
1.1683 
13191 
1.5161 
0.1546 
0.1451 
0.1422 
0. 1454 
0.2295 
0.2172 
0.2151 
0.2221 
0.3162 
0.3035 
0.3046 
0.3186 
0.8650 
0.9884 
1.1408 
1.3289 
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CAPÍTULO 6 


INTERAÇÕES FRACAS E 
ESTABILIDADE QUIRAL! 


Trabalho dos autores publicado no Journal of Quantitative Spectroscopy and 
Radiative Transfer 102, p. 441-449 (2006) e no livro Sense of Beauty in 


Physics: A Volume in Honor of Adriano Di Giacomo, pgs. 197-206 
(Edizione PLUS — Pisa University Press, 2006). 
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6.1 Introdução 


Como é bem conhecido [Barron (1982), Mason (1982)], 
a atividade ótica de um material oticamente ativo muda com o 
tempo. Uma amostra, contendo predominantemente um 
estereoisômero, pode tornar-se uma mistura de iguais 
quantidades de cada isômero. Este processo de relaxação, que é 
chamado de racemização, ocorre espontaneamente ou é devido à 
interação de uma molécula ativa com o meio ambiente. Muitas 
teorias têm sido propostas para descrever a racemização [Vardi 
(2000)]. 

A atividade ótica ocorre [Barron (1982), Mason (1982)] 
quando a molécula tem duas distintas configurações, esquerda 


|E) e direita |D), as quais são degeneradas para o operador 
paridade P=P(x) , isto é, P()/E)=|D) e P(x)|D)=|E). O 
isomerismo esquerda-direita pode ser visto em termos de um 
poço duplo de potencial e os estados |E) e |D) podem ser 
descritos como configurações moleculares que estão 
concentradas à esquerda e à direita do poço de potencial. Os dois 
enantiômeros de uma molécula quiral são descritos por 
superposições de auto-estados de paridade par e ímpar de um 
poço duplo localizado em torno dos mínimos do potencial, 
x=a e x=-a. A coordenada x é inserida no operador 
paridade P(x) e conecta os dois mínimos do potencial. Ela pode 
representar a posição de um átomo, a rotação de um grupo em 
torno de uma ligação, alguma outra coordenada, ou coordenadas 
coletivas de uma molécula. 

Vamos definir por H, a Hamiltoniana de cada lado do 


poço duplo e da barreira de potencial V (x) separando os dois 
mínimos do poço duplo. Nesta situação, E) ê |D) são auto- 


estados de H,, isto é, (EJH,|E)=(D|H,|D)= E, e há uma 
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pequena penetração destes estados interior da barreira V, (x) de 
modo que (EV,|D) = (DV, |E) =). 

Assumiremos que o poço duplo de potencial tem a forma 
de potenciais harmônicos superpostos [Merzbacher (1970)]. 


Indicando por w a frequência fundamental de cada oscilador 
harmônico e por 4 a massa reduzida das partículas vibrando 


entre x=a e x=-a, os estados vibracionais fundamentais 
|D(x)) dos osciladores harmônicos esquerdo e direito são 
escritor, respectivamente, como [Merzbacher (1970)]: 


9,09) = (ua! mi) exp quo! 2h + a), (6.1.1a) 


[9 09)= (ua! mi)“ exp quo! 2h a], (6.1.1b) 


Os estados de configuração esquerdo e direito da 
molécula ativa serão escritos na aproximação de Born- 
Oppenheimer (aproximação adiabática) como 
|E) = vs; )D (x) e |D) = yo Pp (0) , em que |W) descreve 
todos os graus internos de liberdade da molécula ativa exceto x. 


Nestas condições, ó = (EV, (OD) = (D V,OO| E) é dado 


por [Merzbacher (1970)]: 


o =(ho/m7º Juma? h)º expl- uma? fhi). (6.1.2) 


Na década de 1980, experimentos ópticos [Bouchiat e 
Pottier (1986), Noecker, Masterson e Wieman (1988)] 
mostraram casos nos quais a simetria especular (paridade) em 
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átomos estáveis é quebrada durante absorção de luz. Estes 
resultados apóiam a teoria da unificação das forças eletromagnética 
e fraca. A descoberta da violação da paridade em processos 
atômicos foi o resultado de muitos anos de esforço experimental. 
Depois do aparecimento das teorias unificadas no começo dos anos 
1970, muitos experimentos foram projetados para testar novas 
teorias, para escolher entre elas, e para medir as constantes 
fundamentais envolvidas [Bouchiat e Pottier (1986)]. 

Se os efeitos de interação fraca estão presentes, a 
paridade é violada e os lados esquerdo e direito do potencial 
duplo não são mais  simétricos. Deste modo, 


(EH,|E)=E,=E,-e e (D|H,|D)= E, =E,+e, em que 


2 


2e é a diferença de energia entre as configurações direita e 
esquerda devido à interação de violação da paridade. De acordo 
com cálculos realizados por Di Giacomo, Paffuti e Ristori, em 
1980, €/h é tipicamente da ordem de 10 Hz para transições 
rotacionais e vibracionais e da ordem de 10 Hz para transições 
magnéticas nucleares [Harris e Stodolski (1978), Mason e 
Tranter (1984, 1985), Tranter (1985a,b), MacDermott, Tranter e 
Trainor (1992)]. 


Na Secção 6.2 mostraremos como calcular, dentro do 
contexto da equação de Schrôdinger, a racemização quando a 
molécula ativa é submetida a um potencial perturbativo genérico 
U(t). 


6.2 Molécula Ativa Interagindo com o Ambiente 


Em trabalhos precedentes [Cattani (1991, 1992, 1993, 
1994, 1995a,b, 1996), Cattani e Bassalo (1997, 1998, 2001)] 
calculamos, usando a equação de Schrôdinger, a racemização 
quando as moléculas ativas, embebidas em um gás, liquido ou 
sólido, são submetidas a um campo externo genérico U(t). Em 
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nosso modelo assumimos que a racemização é produzida 
essencialmente pelas transições entre dois estados vibracionais: 


|E) e |D). Desta maneira, a função estado [W(1)) da molécula 


ativa é representada por 


PP (1)) = a; (D|E)+a,(D|D), (6.2.1) 


e obedece à equação iho/P(1))/9t=[H, +V,C)+U(O|PCO). 
Assim, a,(t)e ap(t) são governados pelas seguintes equações 
diferenciais: 


da (1) dt= ilha (E, -E+U e) +a((6+U sp), 
(6.2.2a) 


dap()/ dt=«ilhla(0E, +E+U pp) ta (0O(6+U pel, 
(6.2.2b) 


em que os elementos de matriz U,,, com n,k = E,D, são dados 
por U,, = (ju coJk)- 

Como as interações quirais são iguais [Vlot, van 
Miltenburg e Oonk (1997)], definimos u=U,p=Upp. À 


interação heteroguiral será indicada por P=U,, =Upp. Neste 
caso as expressões (6.2.2a,b) são escritas como: 


da (t)/ dt = (ilha (E, -e+u)+a,(N(6+9)], 
(6.2.3a) 
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dap(t)/ dt=ilhla,(O(E, +e+ru)+a (6 + q). 
(6.2.3b) 


Nas próximas Secções 6.3, 6.4 e 6.5, essas expressões 
(6.2.3a,b) gerais serão aplicadas para calcular a racemização e a 
estabilidade óptica para três casos diferentes. 


6.3 Racemização Produzida por Colisões Binárias 
Assumindo « =0 


Nesta secção assumimos que e=0 e que a molécula 
ativa é embebida em um gás diluído ou comprimido ou em um 
líquido em que o potencial U(t)é devido a colisões binárias 
[Cattani e Bassalo (2001), Ch'en e Takeo (1957)] entre as 
moléculas do sistema ativo. Em gases diluídos as colisões têm 
duração muito curta (em torno de 107" seg para um sistema na 
temperatura ambiente) com uma frequência de colisão muito 
alta, para densidades moleculares N =107 /cm”. As colisões 
moleculares induzem transições entre as configurações E e D, 
descritas por qy nas expressões (6.2.3a,b). As transições 
espontâneas entre E e D são descritas por ó. Fazendo e =0 
nas expressões (6.2.3a,b), a, (1) e a, (t), essas equações podem 
ser resolvidas exatamente dando: 


ab) =exp[-i(E,(1)/h+0,40))| x 


x [aexp(-i8,()) +bexp(i8,p(D)]/2, (63.13) 
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ap(1)=exp| -i(E,(1)/h+ 0,:(1)) | x 


x [aexp(-i8,p(1))-bexp(i9,(D)]/2. (6.3.1b) 


em que a e b são constantes a serem determinadas pelas 
condições iniciais e: 


0:(D= J; (ni Vo (2) + U(D |k) dth, (6.3.1c) 

com n,k = E,D. 
Se em t=( a molécula ativa é preparada de modo que 
[E(0)) = E), obtemos, a partir das expressões (6.3.1ab), e 


colocando a,(0)=1 e a,(0)=1, quea=b=1. Portanto, o 
estado da molécula ativa será descrita por: 


[P(1))=exp[-i(E,t/h+ 0,:(1)) | x 


x fcos[6,,(1)]|E)-isen[0,,(D]|D)) . (6.3.2) 


Dessa maneira, a função básica ro=(DjrO) será 
dada por [Cattani (1991,1994,1995a,b)]: 


r0=(DEO) =(sen'(0, (O, (6.3.3a) 
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Bel) = | (EL VC) +U(D]|D) deth+ [ pm) deth, 
(6.3.3b) 


em que (...) significa uma média sobre os efeitos de colisões 
binárias que são aditivas, independentes e randômicas. A 
interação U(t) é escrita como uma soma de interações binárias 
dadas simplesmente por u(t) = 7/ R(t)”, em que 7 é a constante 
de força entre as partículas interagentes, R(t) a distância entre 
elas e tomada como uma função do tempo t, p=4,5,... se a 


interação é dipolo-quadrupolo, quadrupolo-quadrupolo e assim 
por diante, respectivamente. 


No caso de um gás diluído, tratando as colisões binárias 
na aproximação de impacto, temos que [Cattani (1995a,b), 
Ch'en e Takeo (1957)]: 


r(1)=[1-cosQé i)exp(-49]/2, (6.3.4) 


em que A=(y/h) PD N(KkT Im)? Per». Né a densidade de 
moléculas perturbadoras, k é a constante de Boltzmann, T é a 
temperatura absoluta do sistema e m é a massa reduzida das 
partículas colidentes. 

Para um gás comprimido ou um líquido, no qual as 


colisões são quase estáticas, mostramos que [Cattani e Bassalo 
(2001)]: 


r;(1) =| 1-cos(26t/h)exp(-A*1"") |/2, (6.3.5) 


em que 2*=(87/ p)N(y/2h)"" [ xP sen?xdx. 
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Levando em conta as expressões (6.3.4,5), verificamos 
que em gases e líquidos as colisões binárias randômicas entre as 
moléculas ativas e perturbativas sempre produzem a 
racemização da amostra. Comparando r (t) e r,(t) verificamos 
que rn (t), em gases diluídos, decai no tempo com exp(—1t) e em 
gases densos e líquidos r,(t) decai mais suavemente, com 
exp(-1*1*"), pois p=4,5,...e assim sucessivamente. 

Objetivando fazer uma estimativa de n(t) e n(t), 
assumiremos que há somente uma interação dipolo-quadrupolo 
(p=4) entre as moléculas ativa e perturbativa. Neste caso, 
[Cattani (1993)], desde que y=d(E|Q0)|D)=deexp(-uma?/3h), 
À é dado por: 


A=13.0N kTim)'“(0 dih)" exp-2uma” 13h), (6.3.6) 


em que d é o dipolo elétrico da molécula perturbadora e O o 
elemento de matriz de quadrupolo da molécula ativa entre as 


configurações esquerda e direita. Analogamente, 4 * é escrito 
como: 


1*=2.867 N(0 d /2h)* exp(-3uma” | 4h). (6.3.7) 


Então, de acordo com as expressões (6.3.4,5), o sistema 
inevitavelmente racemiza quando e=0 e as colisões entre as 
moléculas ativas e perturbativas são binárias, aditivas, 
independente e randômicas. 

Uma estimativa numérica de n(t) e n(t) será feita 
considerando os seguintes parâmetros moleculares típicos: 
a=10* cm, u=10? g m=10? g d=10" esu,, 
9=10? esu., T=300 Ke N=10"/cm. As fregiiências q 
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serão escritas como wo = Ax10” rad/s. Levando em conta esses 
valores, verificamos que ó/h, A e A*, definidos pelas 
expressões (6.1.2) e (6.3.6,7), são dados por: 


ô/h=1.10x107 4º" exp(-9.524)y”, 
À =5.03x10º exp(-6.354)y” e 
À* = 2.90x 102 exp(-7.144)y”, 


respectivamente, medindo os tempos em anos (y). 

Nas Figuras 6.3.1 e 6.3.2 são mostrados n(t) e n(t), 
respectivamente, como uma função de t, medido em anos, para 
A=4.5. Para um gás diluído, de acordo com a Figura 6.3.2, a 
racemização ocorre para t>1 dia. Conforme a Figura 6.3.1, o 
fator rn (t) para t< 0.002 oscila em torno de 0.5, assumindo 
rapidamente este valor para tempos t > 0.002 y. Como vemos 
na Figura 6.3.2, n(t) oscila por um longo período 
(0<t<10º y) em torno de 0.5 assumindo assintoticamente este 


valor somente para t > 2x10” y . Estes resultados mostram que 
o mecanismo de racemização devido a colisões binárias em 
gases diluídos (aproximação de impacto) é muito mais eficiente 
do que o produzido pelas interações quase estáticas em gases 
comprimidos e líquidos. Note que quando U é estático (ver 
Secção 6.3), não há racemização. 

Em conclusão, quando e=0 e U(t) é criado por 
colisões randômicas binárias, o sistema oticamente ativo 
inevitavelmente racemiza. Isto implica que para tempos muito 
longos a “atividade óptica” ou “rotação óptica” da amostra 


definida por O, (t)=1 — 2r(t) torna-se igual a zero. 
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Fig. 6.3.1. Função (t) para gases diluídos, definida pela expressão (6.3.4), 
em função do tempo £, medido em anos (y). A racemização foi calculada 


considerando a fregiiência de oscilação harmônica ( = 4.5x 10” rad/s. 
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Fig. 6.3.2. Função r,(t) para gases densos e líquidos, definida pela 


expressão (6.3.5), em função do tempo £, medido em anos (y). A 
racemização foi calculada considerando a fregiiência de oscilação 


harmônica O) = 4.5x10º rad/s. 
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6.4 Estabilidade Quiral Gerada por um Ffeito 
Cooperativo entre as Moléculas Interagentes e pelas 
Interações Fracas 


Vamos considerar agora o caso em que £ 0 e a molécula 
ativa é embebida em um gás denso, líquido ou sólido, no qual as 
interações múltiplas dominam as interações binárias e existe um 
efeito cooperativo entre as moléculas interagentes. Devido a este 
comportamento cooperativo assumiremos que cada molécula é 
sujeita a um campo médio resultante dessas interações 
combinadas de todas as outras moléculas do sistema. Esse 
campo médio é entendido como um campo de Hartree 
autoconsistente. 

O potencial de interação cooperativo, que será indicado por 
U(x), será levado em conta no contexto da equação de 


Schrôdinger simultaneamente com a interação fraca. 


Consegiientemente, como Ulx) é estático, as expressões 
(6.2.3a,b) podem ser resolvidas exatamente. Deste modo, se em 
t=0 a molécula ativa é preparada de modo que | (0)) = | E), 
o estado molecular será descrito por [Cattani (1994)]: 


[P(1)) =exp(iEr/h) ([cos” y + sen” gexpliAr ih) || E)- 


— i sen(2%) sen(At/h)|D), (6.4.1) 


1/2 


em que cotg2y=elMó+p) e A=[e" +(6+9)'] 


Sendo sen(2/)=(ó+qp)e” +(6+9)']”, a função 
r(t) será dada por 
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0) =[DIPO) = 0 sen? [e +(6+97]" th), (6.4.2) 


em que O é a amplitude de r(t) dada por 
O=(ó+qp) e” +(6+97)']. 


Assim, quando U é estático, de acordo com a expressão 
(6.4.2), não há racemização e a atividade óptica O,(t) oscila 
com um período dado por T=2mh Me? +(6+ 9) ] em torno 
do valor médio <O, >=1- 0/2. 

Deste modo, para existir uma estabilidade quiral, isto é, 
O, =1, é necessário que O<<l. Então, quando há uma 
interação cooperativa U(x) entre as moléculas da amostra, a 
estabilidade óptica ocorre somente se £ >> (6 +49). 


Vamos considerar o caso particular de gases densos e 
líquidos compostos por moléculas dipolares. Este é um caso 
especial porque o mecanismo de interação cooperativa aparece 
entre as moléculas da amostra e U(x) pode ser facilmente 


calculado [Claverie e Jona-Lasinio (1986)]. Em um meio 
dielétrico [Claverie e Jona-Lasinio (1986), Onsager (1936), 


Bóttcher (1973)] aparece uma campo de reação E, dado por 
E, =Mé-I)d (2é+1)Rº, em que É é a constante dielétrica do 
meio, d é o momento de dipolo médio da molécula ativa e R é 
o raio da cavidade onde d está embebido. Deste modo, o 
potencial de interação entre E, e a molécula ativa é dado por 
U(y)=-d-E , - Uma vez que o elemento de matriz do dipolo da 


molécula ativa entre E) e |D) é zero, a interação heteroquiral 
(EU ()|D)=q desta molécula com E, será calculada levando 


em conta, por exemplo, o momento de quadrupolo O(x) da 
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molécula ativa. Assim, py será dado por p=d (E|OGO|D)! Rº, 
em que (EjO()|D) = 0exp(- uma” /h), sendo 8 o elemento de 


matriz de quadrupolo da molécula ativa entre as configurações 
esquerda e direita [Cattani (1993)]. Desse modo, q é dado por: 


p= (0d /R')expl-uma” !h). (6.4.3) 


Agora, vamos estimar a amplitude de racemização 
O=(ó+q) e” +(6+q)], definida pela expressão (6.4.2), 
levando em conta que ó/h e qy/h são dados, respectiva- 
mente, por 9/h = Qo/1º?) (uwa?/h)'exp(-uma?/h) e p/h= 
=(0d/hR*)exp(—uma? /i), de acordo com as expressões (6.1.2) 
e (6.4.3). Estes dois fatores, usando os parâmetros moleculares 
definidos na Secção 6.3, considerando Rx= 3x10*cm e 
colocando w = Ax10º rad/s, tornam-se: 


ô/h=5.54x102 A“ exp(-9.524) Hz e 
plh=1.51x102 exp(-9.524) Hz. 


Na Figura 6.4.1 as amplitudes O(A), para £/h=10” 


Hz e €e/h=10“Hz, são plotadas como uma função do 
parâmetro 4. 
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Fig. 6.4.1 Amplitude (A) definida pela expressão (6.4.2), em função do 


parâmetro À, definido pela equação O) = AX 10) rad/s. Dois casos limites 
têm sido considerados: £/h = 10º Hz (transições rotacional e vibracional) 


e €/h=10* Hz (transições magnéticas nucleares) 
ç g 


Na Figura 6.4.1 vemos que a amplitude O decresce 
muito rapidamente para fregiuências «w > 4x10" rad/s. Quando 
e/h=10"º Hz (transições rotacionais e vibracionais), O 
localiza-se no intervalo 1.82x10º* <Q <1.67x10?, para 
frequências [47 localiza-se no intervalo 
5x10º <w<6x10Urad/s. Estes resultados mostram que é 
possível conseguir estabilidade Óptica para frequências 


w>5x10ºrad/s. Quando £/h=10*Hz (transições nucleares 
magnéticas), condições de estabilidade são obtidas para 


w>5-5.5x10"rad's. Estes resultados mostram que para 
frequências w > 5.5x10"* rad/s a diferença de energia £ gerada 
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pela força fraca produz estabilidade quiral porque, para estas 
fregiiências, e torna-se muito maior do que ó e q. Nestas 
condições a atividade óptica torna-se igual a um, ou seja, 
O, =1. 

Finalmente, o caso de uma molécula isolada é obtido 
colocando y=0 na expressão (6.4.2). Consegiientemente, a 
amplitude O é dada por O=6" (€” + 6”), mostrando que a 
estabilidade quiral deve ocorrer quando £ >> ó . De acordo com 
a Figura 6.4.1, isto acontece para w>5-5.5x10"rad/s. Se 
€e=0, r(t) é dado por r(t)= sen? (2ót/h), mostrando que há 
uma oscilação E-D com um período T=7mlô. Como 
ô/h=5.54x102 A“? exp(-9.524) Hz, de acordo com a Secção 
6.2, e tomando, por exemplo, A=5, obtemos T = 43 dias. Por 
outro lado, para esse valor de A, se e/h= 10º Hz, teríamos 
estabilidade quiral, pois o parâmetro O é praticamente igual a 
zero, ouseja, 0=1.6x10"2. 

No caso de uma molécula ativa isolada (q = 0), o efeito 
das interações fracas no bloqueio das oscilações E-D, que 
ocorrem quando £ >> ó , pode ser explicado usando a relação de 
incerteza da energia AEAt=h. De fato, como o tempo de 
oscilação espontânea entre as configurações E e D é 
At=T=7mhló, a incerteza na energia é dada por AE =ó. 
Deste modo, se há uma diferença de energia € entre E e D, as 
transições naturais E-D são permitidas somente quando 
AE =60>€. Por outro lado, as transições serão proibidas 
quando £ >> 6. Na presença de uma energia de interação q, 
usando o mesmo raciocínio, as transições E-D são bloqueadas 
quando a condição £ >> 6 + q é obedecida. 
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6.5 Racemização e Estabilidade Quiral quando 
e *z0e U(t) é Gerada por Colisões Randômicas 


Na Secção 6.3 calculamos r(t) assumindo que £=0 e 
que o potencial U(t) dependente do tempo, era gerado por 
colisões binárias randômicas. Na Secção 6.4, determinamos as 
condições para a estabilidade quiral assumindo que £ £ 0 e que 


a molécula ativa é submetida a um potencial estático U(x) 
devido a um efeito cooperativo entre as moléculas interagentes 
do sistema. 

Nesses dois casos, as expressões (6.2.3a,b) podem ser 
resolvidas exatamente. No primeiro caso (Secção 6.3), €=0 e 
U(t) depende do tempo. No segundo caso (Secção 6.4), 20 
e U=Ulx) é independente do tempo. Por outro lado, quando 
c*0 e U=U(t), as expressões (6.2.3a,b) não podem ser 
resolvidas exatamente. Nestas condições, verificamos que r(t), 
no caso de colisões binárias randômicas, é dado por: 


r(b)=(6/A,)[-cos(2A t/h)expl-f(DIV/2, (6.5.1) 
em que A,=(€' +67)", f(t)=AM para gases diluídos e 
f()=A*t'” para gases comprimidos e líquidos, em uma 
aproximação quase estática, de acordo com a Secção 6.3. 

Vamos considerar três casos particulares da expressão 
(6.5.1). (1) Quando £=0 obteremos r,(t) e 1,(t) dados pelas 
expressões (6.3.4,5), respectivamente, que descreve a 
racemização da amostra. (2) Quando £>>ó vemos que 
r(t)=0, que é a condição para a estabilidade quiral. (3) Quando 
£ * 0 e U=0, obtemos r()=(6/A,) sen"(A,t/h) que é o 
caso de uma molécula isolada analisada na Secção 6.4. 

Definindo 1. =0.5(6/A,)”, verificamos, por 
intermédio da expressão (6.5.1), que a atividade óptica varia no 
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intervalo 1-2r,, <O, <1. Também vemos que para tempos 
muito longos, isto é, para f(t)>>1, r(t) tende assintoticamente 
para um valor máximo r, .. Assim, há racemização, isto é, 


(9) E =0, somente quando ô>>€, pois 1, =0.5. Por outro 


mi 


lado, há estabilidade quiral, isto é, O 5 =1, quando £ >> 6 pois 


Fax = O. Vamos estimar os valores de w para esses dois casos 


extremos assumindo que €/h=10º Hz e tomando 
ô/h=5.54x102 A“? exp(-9.524) Hz, dado na Secção 6.3, 
relembrando que A é definido por w = Ax10"rad/s. Assim, 
verificamos que a estabilidade quiral deve ocorrer quando 
A>5.2, isto é, para w>5.2x10ºrad/s e que há racemização 
quando A <4, isto é, para o<4x 10" rad/s. Para freqiiências no 
intervalo 4x10º <w<5.2x10 rad/s r(t) oscila e decresce 
quando o tempo aumenta tendendo assintoticamente para 
Fax QUE se encontra no intervalo O<r,. <1/2. 

Na Figura 6.5.1 a função r(t), definida pela expressão 
(6.5.1), é mostrada como função do tempo no caso de gases 
diluídos para £=10“e A=4.7. De acordo com esta figura, 
r(t) oscila entre O e 0.082 e estabiliza com r,.. = 0.041 para 


t> 10” anos - 3.6 dias, dando uma atividade óptica assintótica 
O, = 0.918. 
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Fig. 6.5.1. Função r(t) definida pela expressão (6.5.1) como função do 
tempo £, medido em anos, para£/h = 10 Hz e A=4/7 


6.6 Sumário e Discussões 


Nas Secções 6.3 e 6.5 calculamos r(t) assumindo que a 
interação U entre as moléculas ativa e perturbativa é devida a 
colisões binárias, aditivas, independentes e randômicas. 
Mostramos que quando £ =0 o sistema óptico inevitavelmente 
racemiza, isto é, O 5 =(0 para tempos muito longos. Por outro 


lado, se 0, há racemização somente quando ó>>€e e há 
estabilidade quiral (O, =1) somente quando £ >> ó. 


Como mostramos na Secção 6.4, quando U é criado por 
um efeito cooperativo no sistema ativo não há racemização 
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independentemente do valor de £ e há estabilidade quiral 


2 


somente quando E>>ó+q. Isto é, quando a diferença de 


energia € é suficientemente grande para bloquear as transições 
E-D induzidas pelo tunelamento natural e pelo potencial 
cooperativo. 

Em nosso modelo, a interação cooperativa U(x) é 
entendida como um campo de Hartree autoconsistente e foi 
levado em conta no contexto da equação de Schrôdinger, 
simultaneamente com a interação fraca. Em nossa aproximação, 
a simetria esquerda-direita é quebrada pela interação fraca, a 
qual é responsável, sob algumas condições, pela estabilidade 
óptica. 

Em um modelo proposto recentemente por Vardi (2000) 
para estudar a estabilidade quiral, o campo autoconsistente tem 
duas componentes: U e U,..» geradas pelas interações 


hom 
homoquiral e heteroquiral, respectivamente. Essas componentes 
foram introduzidas em uma equação de Schrôdinger não-linear 
com o propósito de promover a evolução temporal do sistema 
ativo. Ele mostrou que, quando interações U são 


hom 


energicamente favoráveis em relação às interações U,., há uma 


het > 
quebra espontânea de simetria esquerda-direita que viria 
amplificar a atividade óptica de uma mistura quase racêmica. 

A Mecânica Quântica Não-Linear tem sido usada [Vardi 
(2000), Koschany, Kuffer, Obermaiar e Plessner (1994)] para 
explicar a estabilidade quiral. Esta parece ser uma tentativa 
plausível porque os estados estacionários de uma equação de 
Schrôdinger não-linear [Koschany, Kuffer, Obermaiar e 
Plessner (1994), Smerzi, Fantoni, Giovanazzi e Shenoy (1997)] 
não necessitam ser auto-estados dos operadores relativos ao 
grupo de simetria do potencial. Assim, o termo não-linear 
introduz uma quebra espontânea de simetria [Smerzi, Fantoni, 
Giovanazzi e Shenoy (1997), Davies (1981)] que favorece a 
localização em um dos poços. Contudo, equações de 
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Schrôdinger não-lineares realistas devem ser deduzidas 
levando em consideração os efeitos cooperativos em interações 
de muitos-corpos na amostra [Smerzi, Fantoni, Giovanazzi e 
Shenoy (1997), Davydov (1982), Davies (1981), Yomosa 
(1978)]. Este algoritmo permitiria obter uma fiel equação não- 
linear para estudar a estabilidade óptica. As equações não- 
lineares consideradas por Vardi (2000) e Koschany, Kuffer, 
Obermaiar e Plessner (1994) não foram obtidas dessa maneira. 
Eles propuseram, de um modo arbitrário, equações seguindo 
modelos não-lineares genéricos adotados na literatura [Krenke e 
Campbell (1986)]. Além disso, sabemos que equações não- 
lineares exibem um número grande de soluções ricas e 
complexas dependendo da magnitude dos parâmetros não- 
lineares. Assim, da análise de Vardi (2000) e Koschany, Kuffer, 
Obermaiar e Plessner (1994), é difícil concluir que os efeitos 
não-lineares são ou não realmente mecanismos efetivos 
responsáveis pela estabilidade quiral. 

Finalmente, é importante observar que nossas conclusões 
a respeito da estabilização dos enantiômeros são limitadas às 
moléculas que principalmente racemizam por intermédio de 
apenas uma inversão simples entre dois estados (E e D). 
Entretanto, há muitos outros mecanismos de racemização [Eliel 
e Wilen (1994)]. Em nosso trabalho, esses processos não foram 
considerados. 
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